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Abstract 

The problem of anisotropy and its effects on the statistical theory of high Reynolds- 
number (Re) turbulence (and turbulent transport) is intimately related and inter- 
mingled with the problem of the universality of the (anomalous) scaling exponents 
of structure functions. Both problems had seen tremendous progress in the last 
five years. In this review we present a detailed description of the new tools that 
allow effective data analysis and systematic theoretical studies such as to separate 
isotropic from anisotropic aspects of turbulent statistical fluctuations. Employing 
the invariance of the equations of fluid mechanics to all rotations, we show how to 
decompose the (tensorial) statistical objects in terms of the irreducible representa- 
tion of the SO(d) symmetry group (with d being the dimension, d = 2 or 3). This 
device allows a discussion of the scaling properties of the statistical objects in well 
defined sectors of the symmetry group, each of which is determined by the "angu- 
lar momenta" sector numbers (j, m). For the case of turbulent advection of passive 
scalar or vector fields, this decomposition allows rigorous statements to be made: (i) 
the scaling exponents are universal, (ii) the isotropic scaling exponents are always 
leading, (iii) the anisotropic scaling exponents form a discrete spectrum which is 
strictly increasing as a function of j. This emerging picture offers a complete un- 
derstanding of the decay of anisotropy upon going to smaller and smaller scales. 
Next we explain how to apply the S0(3) decomposition to the statistical Navier- 
Stokes theory. We show how to extract information about the scaling behavior in 
the isotropic sector. Doing so furnishes a systematic way to assess the universality 
of the scaling exponents in this sector, clarifying the anisotropic origin of the many 
measurements that claimed the opposite. A systematic analysis of Direct Numerical 
Simulations (DNS) of the Navier-Stokes equations and of experiments provides a 
strong support to the proposition that also for the non-linear problem there exists 
foliation of the statistical theory into sectors of the symmetry group. The exponents 
appear universal in each sector, and again strictly increasing as a function of j. An 
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approximate calculation of the anisotropic exponents based on a closure theory is 
reviewed. The conflicting experimental measurements on the decay of anisotropy 
are explained and systematized, showing agreement with the theory presented here. 
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1 Introduction 



The statistical theory of fluid turbulence is concerned with correlation func- 
tions of the turbulent velocity vector field u{x, t) where x is the spatial po- 
sition and t the time [1] . Since the velocity field is a vector, multi-point and 
multi-time correlation functions are in general tensor functions of the vector 
positions and the scalar times. Naturally such functions have rather compli- 
cated forms which are difficult to measure and to compute. Consequently, 
almost from its very beginning, the statistical theory of turbulence was dis- 
cussed in the context of an isotropic and homogeneous model. The notion of 
isotropic turbulence was first introduced by G. I. Taylor in 1935 [2]. It refers 
to a turbulent flow, in which the statistical averages of every function of the 
velocity fleld and its derivatives with respect to a particular frame of axes 
is invariant to any rotation in the axes. This is a very effective mathematical 
simplification which, if properly used, can drastically reduce the mathematical 
complexity of the theory. For this reason, it was very soon adopted by others, 
such as T. D. Karman and L. Howarth [3] who derived the Karman-Howarth 
equation (see below), and A. N. Kolmogorov [4,5] who derived the 4/5 law 
(re-derived below). In fact, most of the theoretical work in turbulence in the 
past sixty years was limited to the isotropic model. 

Experimentally, however, we know that isotropy holds only as an approxima- 
tion with a varying degree of justification. In all realistic fiows there always 
exists some anisotropy at all scales; the statistical properties of the velocity 
fleld are effected by the geometry of the boundaries or the driving mechanism, 
which arc never rotationally invariant [6-14]. Therefore, a realistic description 
of turbulence cannot be purely isotropic and must contain some anisotropic 
elements. Yet the problem is that once we take anisotropy into account, we 
face a drastic increase in the complexity of the theory. The number of vari- 
ables that is needed to describe the common statistical quantities, such as 
correlation functions and structure functions of the velocity fleld, increases a 
lot. For example, under shear there is a characteristic length scale, which can 
be constructed from the typical velocity and the typical shear [8,15]. This 
length has to be considered in order to distinguish those scales where the 
turbulent evolution is mainly dominated by the inertial effects of fluid me- 
chanics or by the direct input of energy due to the anisotropic shear [16-19]. 
Similarly, all dimensional estimates acquire a signiflcant degree of ambiguity 
because of the proliferation of different dimensional quantities related to the 
parameters of anisotropy. As a consequence of these inherent difficulties the 
existing anisotropic effects were simply ignored in many of the experimental 
and simulational studies of statistical turbulence. This attitude gave rise to 
ambiguous assessments of important fundamental issues like the universality 
of the scaling exponents in turbulence. 
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The standard justification for ignoring anisotropic effects is that the basic 
phenomenology, since the pioneering works of Kolmogorov [4, 5] , predicts a 
recovery of isotropy at sufficiently small scales of the turbulent flows. Nev- 
ertheless, both recent experimental works and theoretical analysis suggested 
that the actual rate of recovery is much slower than predicted by simple di- 
mensional analysis, pointing out even the possibility that some anisotropic 
correlation function, based on velocity gradients, stays 0(1) for any Re [7,20- 
22]. In order to settle this kind of problems, theoretically or experimentally, it 
is crucial to possess systematic tools to disentangle isotropic from anisotropic 
fluctuations and to distinguish among different kinds of anisotropic fluctua- 
tions. Thus a central challenge in the theory of anisotropy in turbulence is 
the construction of an efficient mathematical language to describe it. Without 
a proper description, the complexities of the formalism can soon obscure the 
physical content of the processes that we wish to study. 

The problem of anisotropy is not disconnected from the other fundamental 
problem which has to do with the nature of universality in turbulence. By 
universality, we mean the tendency of different turbulent systems to show the 
same small-scales statistical behavior when the measurements are done far 
away from the boundaries. Consider, for example, the longitudinal two-point 
structure function 

S^^\r) = (^5u'j{x,r,t)'^ , 5ue{x,r,t) = r ■ [u{x + r,t) - u{x,t)] , (1) 

with r being the unit vector in the direction of r, and (•) stands for an appro- 
priate ensemble average. This function shows essentially the same dependence 
on the separation vector r, whether it is measured in the atmospheric bound- 
ary layer, in a wind tunnel or in a DNS, provided it is measured for sufficiently 
small separations and far from the boundaries. This high-degree of universality 
cannot be expected if anisotropic fluctuations were the dominant contributions 
to the two-point structure functions. Different boundary conditions and differ- 
ent forcing mechanism necessarily introduce different large-scale anisotropics 
in the flow, which would translate to different small-scale anisotropic fluctua- 
tions. Small-scale universality can be achieved only if anisotropic fluctuations 
are sub-leading with respect to the isotropic fluctuations. In the following, we 
also discuss which aspects of the anisotropic fluctuations are universal and 
which arc not. We will see that some aspects of the anisotropic fluctuations 
depend on the boundary conditions while other aspects do not. In fact we will 
show that scaling exponents are expected to be universal whereas amplitudes 
depend on the boundary conditions. 

In the last 5 years, a tremendous progress in the understanding of the two 
aforementioned problems, i.e., finding a mathematical language that prop- 
erly describes anisotropic turbulence and its universal properties has been 
achieved. Not surprisingly, the two problems are closely related, as it often 
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happens in physics - a problem becomes considerably simpler if described in 
the proper mathematical language. The technical core of these recent achieve- 
ments is the S0(3) decomposition [23]. This tool enjoys the advantages of 
being mathematically simple, yet very powerful and systematic. By using it, 
many of the mathematical complexities of dealing with anisotropy in turbu- 
lence and in other hydrodynamic problems arc greatly simplified. The principal 
idea is to represent the main statistical observables, such as structure func- 
tions and correlation functions, in terms of their projections on the different 
(j, m) sectors of the group of rotations. It can be applied to all the statistical 
quantities in turbulence, creating a detailed profile of the effects of anisotropy. 
Additionally, and perhaps more importantly, the SO (3) decomposition reveals 
some new universal properties of fully developed turbulence. It is expected 
that each sector of the SO (3) group has its own universal exponents. In par- 
ticular, it is shown that the exponents associated with the anisotropic sectors 
are larger than the isotropic exponents, in accordance with the isotropization 
of the statistics as smaller and smaller scales are observed. 

As already mentioned, the SO (3) decomposition is useful also to investigate 
isotropic and anisotropic fiuctuations in other hydrodynamic problems. In 
particular we will focus on the case of scalar and vector quantities passively 
advected by a turbulent velocity field. In these cases, one may often elevate the 
phenomenological assumptions made for turbulent anisotropic fluctuations to 
the status of rigorous statements [24]. By using a systematic decomposition 
in different sectors of the SO (3) group one may show that passive scalars, 
advected by stochastic self-similar Gaussian velocity fields, always possess 
isotropic leading small-scale fluctuations. Moreover, one may quantitatively 
distinguish among different kinds of anisotropics, assessing their rate of decay 
by going to smaller and smaller scales. It turns out that the rate of recovery 
of isotropy is typically much slower than expected on the basis of dimensional 
analysis. Moreover, all different anisotropic fluctuations decay in a self-similar 
way but with different rates; the scaling exponents being universal, while pref- 
actors are non-universal [25]. The very same can be rigorously proved for the 
passive advection of vector-like quantities, as for the case of magnetic flelds 
when the feedback on the velocity evolution due to the Lorentz force is ne- 
glected. There, the vector nature of the transported quantity leads to an even 
richer, and more complex, list of possible anisotropic fluctuations [26, 27]. An- 
other important problem which we address in detail is the case of the passive 
advection of a vector-like incompressible quantity, i.e. a passive vector with 
pressure [28]. Although without any counterpart in nature, such a system is 
particularly interesting because it can be seen, for some aspects, as the closest 
linear approximation to the non-linear Navier-Stokes evolution. For example, 
it allows to study in a systematic way some problems connected to the con- 
vergence of integrals involving the pressure term. Similar technical problems 
arise also in the analysis of both isotropic and anisotropic multi-point velocity 
correlations in Navier-Stokes equations. 
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Of course, a significant part of this review will be devoted to applications of 
the theoretical and technical tools to physical experimental data [10, 13, 29- 
33] and numerical data sets [15,34-38]. In order to exploit the entire poten- 
tiahty of the S0(3) decomposition one needs to measure the whole velocity 
field, u{x), in a 3 dimensional volume. This is because in order to disentangle 
different projections on different sectors one needs to integrate the given cor- 
relation function against the proper eigenfunction of the rotation group on the 
3d sphere of radius r. By doing that, the exact projection on each different 
sector of the SO (3) decomposition is under control, with the only practical 
limitations for reaching highly anisotropic sectors being the lack of resolution 
of highly fiuctuating angular properties. At the present stage of experimental 
capabilities the exact decomposition can be carried out explicitly only in data 
sets coming from DNS. Here the velocity field in the whole testing volume is 
available. For experimental data, the best way to exploit the SO (3) decom- 
position is to either select observables with vanishing isotropic components, 
in order to focus directly on anisotropic sectors, or to perform a multi-sector 
analysis, i.e. to fit simultaneously the isotropic and anisotropic components. 



The review is organized as follows. Sect. 2 offers a historical review of isotropic 
turbulence. We present a modern derivation of the exact results pertaining to 
the 3'rd order structure function and the celebrated 4/5 law. We review the 
standard theory for all correlation functions, and discuss the experimental 
difficulties with the isotropic theory. These difficulties included apparent per- 
sistence of anisotropies into the small scales for high Re, apparent location 
dependent scaling exponents etc. In Sect. 3 we review the history of attempts 
to deal with anisotropy. In Section 4, the technical basis of the SO (3) de- 
composition is introduced focusing on the particular statistical problems of 
anisotropic fiuctuations discussed in the previous section. Then, in Sect. 5 we 
switch to study exactly solvable hydrodynamic problems with emphasis on ei- 
ther those aspects peculiar to each different model and to those features in 
common with the non-linear Navier-Stokes case. Among the common aspects 
we cite the possibility to study in these models in full details the foliation of 
the equations of correlation functions in different anisotropic sectors; the uni- 
versality of isotropic and anisotropic exponents; the hierarchical organization 
of exponents -leading to recovery of small-scales isotropy. At the end of this 
section we present closure results for two-point turbulent structure function 
in the anisotropic sectors j = 2,4,6. In Sect. 6, the utility of this language 
is demonstrated by discussing experimental data in atmospheric boundary 
layer and on homogeneous-shear flows. In Sect. 7, we present the analysis of 
anisotropy in DNS of typical strongly anisotropic flows. Two cases are dis- 
cussed in depth: channel flows and random Kolmogorov flows, the latter being 
homogeneous flows stirred at the large scales. Sect. 8 presents a summary and 
conclusions. Technical details are collected in the appendices. 
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2 Historical Review: Isotropic Turbulence 



In the first two sections we present a historical review. We start with the model 
of homogeneous isotropic turbulence, and then turn to previous attempts to 
treat theoretically anisotropy in turbulence. 



2.1 Homogeneous and Isotropic Turbulence 

The Navier-Stokes equations for the velocity field are invariant to all rotations: 

^^If ' + [u{x, t) ■ V] u{x, t) = -Vp{x, t) + uA^u{x, t) , (2) 
ot 

V • u{x,t) = 0, 

with p{x, t) and u being the pressure and kinematic viscosity respectively. 
Since the gradient and Laplacian operators arc both rotationally invariant, 
the rotation symmetry of the equation can be broken only by anisotropic 
forcing terms or anisotropic boundary conditions. Rather naturally then the 
statistical theory of turbulence was mostly developed in the framework of 
isotropic turbulence [2]. The central idea of this approximation is that the 
statistical average of any function of the velocity components in any coordinate 
system is unaltered if this coordinate system is rotated or reflected in any 
manner. The assumption of isotropy was widely adopted. In 1938, Karman and 
Howarth [3] used it to explore the second- and third-order correlation- functions 
of the velocity field. Their use of tensor notation was more elegant and compact 
than that used by Taylor. It enabled them to derive some constraints on these 
correlation functions and express them in terms of a few scalar functions. For 
example, for the second-order correlation- function in homogeneous turbulence 

C'^ir, t) = (u'^ix + r, t)u^{x, t)) , (3) 

they used the representation 

C"^(r, t) = [/(r, t) - g{r, i)]Pr^ + g{r, t)5-^ , (4) 



and then derived a linear differential relation between /(r, t) and g{r, t) using 
the solenoidal condition of daC°'^{r, t) — 0, 

2/(r,i)-2^(r,t) = -r^^. 
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This means that under the assumption of isotropy, and using the solenoidal 
condition, the second-order correlation-function can be written in terms of one 
scalar function instead of nine. Similarly, Karman and Howarth analyzed the 
third-order correlation function by representing it as an isotropic tensor and 
then reducing the number of scalar functions using the solenoidal condition. 
They were also able to connect it to the second-order correlation function in 
decaying turbulence using the Navier-Stokes equations. These computations 
have since found their way into every standard text-book on the statistical 
theory of turbulence. 

The mathematical representation of isotropic turbulence has reached its most 
elegant and powerful form in a paper by H. P. Robertson from 1940 [39]. 
Robertson provided a systematic way to represent isotropic tensors using the 
theory of invariants. For example, to derive the general representation (4) in 
the stationary case using Robertson's method, we consider the scalar function 

C(a,6,r) = C"^(r)a«6^, 

with a and h being two arbitrary vectors. If C"^(r) were an isotropic tensor, 
C(a, 6, r) would preserve its functional form upon an arbitrary (simultaneous) 
rotation of the three vectors r, a, h. Using invariant theory, Robertson deduced 
that C(a, b, r) must be a function of the six possible scalar products (r-r), (r- 
a),... and of the determinant [rab] = enapr^a'^h^ . Additionally, by definition, 
it must be a bilinear function of a and b and therefore must have the following 
form: 

C(a, 6, r) = A{r){r • a)(r • 6) + B{r){a ■ b) + C{r)[rab] , 

where A{r), B[r) and C(r) are arbitrary functions. Finally, recalling that 
C(a, 6, r) is the contraction of C°'^{r) with a and 6, we find that 

If we further demand C"^(r) to be invariant to improper rotations as well 
(i.e., rotations plus reflections), we can drop the skew-symmetric part e'*"^, 
thus retaining a representation which is equivalent to Eq. (4). 

2.2 The 4/5 law in Isotropic Turbulence and its Generalization 

By using the isotropic representation of the third-order correlation function, 
in 1941 Kolmogorov proved the "four-fifth law" well inside the inertial range 
of a fully developed turbulence. This law pertains to the third order moment 
of longitudinal velocity differences, stating that in homogeneous, isotropic and 
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stationary turbulence, in the limit of vanishing kinematic viscosity i/ — > 

{[dui{x,r,t)f) = —er , 

where e is the mean energy flux per unit time and mass, e = u (|Vq,M/3P)- The 
fundamental assumption needed to derive this law is the so-called "dissipation 
anomaly" which means that the dissipation is finite in the hmit z/ — 0. As 
noted in [5], "this is one of the most important results in fully developed 
turbulence because it is both exact and nontrivial. It thus constitutes a kind of 
'boundary condition' on theories of turbulence: such theories, to be acceptable, 
must either satisfy the four-fifth law, or explicitly violate the assumptions 
made in deriving it" . 

To demonstrate how isotropy helps in deriving this result, we present a re- 
derivation in which we will obtain an additional exact relation that appears 
to have the same status as the four-fifth law, pertaining to homogeneous, 
stationary and isotropic turbulence with helicity [40, 41]. Defining the velocity 
v{Xj t) (isv{x,t) = u{x, t) — (u) we consider the simultaneous 3rd order tensor 
correlation function which depends on two space points: 

J°'^T(r) = (v"{x + r, t)v^{x, t)v^{x, t)) . (5) 

We show that in the limit — > 0, under the same assumption leading to the 
fourth-fifth law, this correlation function reads [41] 

ja,/37(^) = -{^(r'SafS + r^K, - ^r"5^,) - ^(^^r^ + e^,,r^)r' , (6) 

where 5^/3 is the Kronecker delta and e^^^ is the fully antisymmetric tensor. 
The quantity h is the mean dissipation of helicity per unit mass and time, 

/i = i/((V"iiO(V"[V X uf)) , 

where repeated indices are summed upon. In the derivation below it assumed 
that h remains constant when i/ — > in the same spirit of the dissipation 
anomaly [42-44]. The first term in Eq. (6) is just the 4/5 law. The new part 
of result (6) can be also displayed in a form that depends on h alone by 
introducing the longitudinal and transverse parts of u: the longitudinal part 
is Ui = r{u ■ r)/r'^ and the transverse part is Ut = u — Ui. In terms of these 
quantities we can present a "two fifteenth law" 

2 

{[5ui{x, r, t)] ■ [ut{r + x,t)x Ut{x, t)]) = —hr'^ . (7) 

io 

We note that this result holds also when we replace u by v everywhere. 
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To derive the result (6) we start from the correlation function J^'^'^^r) which 
is symmetric with respect to exchange of the indices (3 and 7 as is clear from 
the definition. In an isotropic homogeneous medium with helicity (no inversion 
symmetry), the most general form of this object is [41]: 



+b2{r)[ea/3sr^ + e^^sr^y + a3(r)[5a^r^ + + Sp^r'' - 5r"r^rVr'] . 

This general representation is invariant to the choice of orientation of the 
coordinates. Not all the coefficients are independent for incompressible flows. 
Requiring dJ"'^'^{r)/dr" — leads to two relations among the coefficients: 



+ Ih^''^ = + ' + ^)Ni(^) - 4ai(r)] = . 



As we have two conditions relating the three coefficients ai, ai and as only one 
of them is independent. Kolmogorov's derivation related the rate of energy dis- 
sipation to the value of the remaining unknown. Here the coefficient 62 remains 
undetermined by the incompressibility constraint; it will be determined by the 
rate of helicity dissipation. 

Kolmogorov's derivation can be paraphrased in a simple manner. Begin with 

the second order structure function S^'^\r) = {\u{x + r) — u{x)\'^) . Computing 
the rate of change of this (time-independent) function from the Navier-Stokes 
equations (2) we find 

= = -P(^) (r) - 2-e + .V^^(^) (r) , (8) 



where T>^'^\r) stems from the nonhnear term (u- V)u and as a result it consists 
of a correlation function inchiding a velocity derivative. The conservation of 
energy allows the derivative to be taken outside the correlation function: 

P(2)(r) = -^{u''{x,t)u''{x + r,t)[u^{x,t) -u^{x + r,t)]) . (9) 



In terms of the function of Eq. (5) we can write 

p(2)(^) = _^ [j"'/3"(t., t) - r'P'^i-r, t)] . (10) 



Note that Eq. (5) is written in terms of v rather than u, but using the incom- 
pressibility constraint we can easily prove that Eq. (9) can also be identically 
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written in terms of v rather than u. We proceed using Eq. (8) in Eq. (10), 
and find 

= 2-^r^[5a^{r) + 25i(r)] . (11) 



For r in the inertial interval, and for z/ — > 0, we can read from Eq. (8) V^'^\r) — 
— 2e and therefore have the third relation that is needed to solve all the three 
unknown coefficients. A calculation leads to 

ai(r) = -2e/45 , Oi = -e/18 , 03 = . 

The choice of the structure function 5*2 (r) leaves the coefficient h2{r-) unde- 
termined, and another correlation function is needed in order to remedy the 
situation. Since the helicity is u- [V x w], we seek a correlation function which 
is related to the hehcity of eddys of scale of r: 

T(2)(r) = {[u{r + x,t)- u{x, t)] ■ [V x u{x + r, t) - V x u{x, t)]) . 

Using the Navier-Stokes equations to compute the rate of change of this quan- 
tity we find 

= ^^5^ = -G^'\r) -2h- uV^T^^\t) , (12) 



which is the analog of (8), and where 



G(2)(r) = {{u{x, t) ■ [Vr X [u{x + r,t) X [V^ x u{x + r, t)]]])} 

-I- {term r — > —r} . (13) 

The conservation of helicity allows the extraction of two derivatives outside 
the correlation functions. The result can be expressed in terms of the definition 
(5): 

Substituting Eq. (8) we find 

which is the analog of Eq. (11). Using Eq. (12) in the inertial interval in the 
limit u ^ we find the differential equation 

od%2(r) ^ db2(r) , , h 
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The general solution of this equation is b2{r) = —h/30+air~^+a2r~^. Requir- 
ing finite solutions in the limit r — > means that ai = 02 = 0. Accordingly we 
end up with Eq. (6). The moral of this example is that even in isotropic and 
homogeneous systems there exist sub-leading terms which can become domi- 
nant for specially selected objects hke (7). Once anisotropy exists, there are 
many more (in fact infinitely many) sub-leading contributions that need to 
be assessed carefully. Similar results for slightly different correlation functions 
have also been found in [45, 46] 



2.3 Kolmogorov's Theory for 2'nd, 4'th, and Higher Order Structure Func- 
tions 

Unfortunately, the exact result pertaining to the 3'rd order structure function 
is rather unique. The moment we consider 2'nd, 4'th or higher order correla- 
tion functions there is no exact result for the scaling exponents. Kolmogorov's 
1941 theory states that, for Re large enough, small-scales turbulent fluctua- 
tions should recover isotropy and homogeneity (if measured far enough from 
boundaries) and should possess universal scaling properties depending only on 
the mean energy flux, e. For homogeneous and isotropic ensembles one deflnes 

S^-\r) ^ {Sunx,r,t)) = (6-r)t/(-)(f , ^), (14) 

Lo r 

where Lq and rj are the integral length scale and the viscous scale, respec- 
tively. The function f^'^^ is supposed to be well behaved in the hmit of inflnite 
Reynolds numbers for flxed separation, r: ]imx^y-^o f^'"\x,y) — const. In this 
limit, the celebrated K41 scaling prediction for structure functions in the in- 
ertial range, 77 <C r <C Lq, follows: 

^(")(r)~C('^)(er)^^"\ with C^") = -. (15) 

3 

In (15) the constants C^"^ depend only on the large scale properties. Because 
of stationarity, the mean energy flux in Eq. (15) can be equally taken to be 
the mean energy input or the mean energy dissipation. 

Kolmogorov's theory goes beyond the scahng prediction (15). For example, 
any non-vanishing pth order structure functions, including purely transversal 
and mixed longitudinal-transversal velocity increments, must possess the same 
scaling exponents: 

5("'™)(r) = ((5<(r)5Mr(r)) - C("'"^)(er)^, (16) 
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where p = n + m, Sut{r) = 5u{r) — 5u£{r)r, and Sut{r) is one of the compo- 
nents of the two-dimensional transverse velocity difference. Notice that due to 
the assumption of isotropy, only even combinations of transversal increment 
in (16) have a non vanishing average. It is also not difficult to extend the 
K41 reasonings to describe also correlation functions at the viscous scales, i.e. 
observables based on gradients statistics [47, 48] . 

2.4 Experimental Difficulties with the Isotropic Theory 

On the whole, experimental tests of Kolmogorov's theory ran into increasing 
difficulties when the data were analyzed with greater detail. The first system- 
atic attempt to check the isotropic scaling (15) for high Re number turbu- 
lence was [49]. These authors performed a high statistical test of K41 theory 
by going beyond the usual two-point correlations. They measured structure 
functions of higher order, reaching good evidence that there exist anomalous 
deviations from the scaling exponents (15). Their data substantiate a power- 
law behavior with C(n) 7^ n/3. At that time, and for many year later, the 
situation was very controversial both theoretically and experimentally. There 
appeared many claims that the observed deviations were due to sub-leading 
finite-Reynolds effects. One should not underestimate the difficulties of get- 
ting reliable estimates of the scaling exponents. First, one must expect finite 
Reynolds numbers corrections which may strongly reduce the inertial range 
where scaling laws are expected or introduce anisotropic corrections to the 
isotropic K41 predictions. Both effects are usually present in all experiments 
and numerical simulations. Nowadays, state-of-the-art experiments of turbu- 
lence in controlled geometries reach a maximum Re numbers measured on the 
gradient scales, A"^ = (| Vu|)/(|tt|), of Rx ~ 5000 where ^ ^ and U 
is the typical large scale velocity. In atmospheric flows Rx can be as high as 
Rx ~ 20000 but at the expense of high anisotropy. More complex is the situ- 
ation of DNS where the best resolution ever reached up to now is 4096'^ [50], 
corresponding to a i?A ~ 1100. DNS allow a minimization of the anisotropic 
corrections, by implementing periodic boundary conditions and fully isotropic 
forcing, something which is not experimentally feasible. However, also in DNS 
the discrete symmetries induced by the finite lattice-spacing do not allow for 
perfect isotropic statistics. We thus either have high-Reynolds-numbers exper- 
iments which are strongly perturbed by anisotropic effects, or DNS isotropic 
flow at moderate Reynolds numbers. Therefore, one has to face the problem of 
how to disentangle isotropic from anisotropic fluctuations and how to extract 
information on the asymptotic scaling with a flnite -often short- inertial-range 
extension. Only recently, after many experimental and numerical confirma- 
tions of the results of [49], the situation became clearer [51]. We may affirm 
now with some degrees of certitude that the isotropic scaling exponents are 
anomalous, the K41 prediction ({n) = n/3 is wrong, except for n = 3 which 
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is fixed to be C(3) = 1 by the exact 4/5 law. Moreover, the possibihty to show 
analytically the existence of anomalous scaling in turbulent advection [24], 
definitely eliminated those arguments supporting the impossibility to have a 
Re-independent anomalous scaling in any hydrodynamic system. Prom a phe- 
nomenological point of view, it is easy to extend the K41 theory such as 
to include anomalous scaling. Already Kolmogorov noticed, after Landau's 
criticism in 1962, that it is unrealistic to expect the isotropic inertial range 
fluctuations to depend only on the mean energy dissipation, e. Kolmogorov 
proposed in 1962 [52] to employ the coarse-grained energy dissipation over a 
box of size r, 

= ^ / dye{x + y), (17) 

\y\<r 

to match the correct dimensions of structure functions in (14), the so-called 
Refined Kolmogorov Hypothesis: 

S^''\r)^C^''\€^{r))r^. (18) 

This hypothesis connects the deviation from the K41 prediction, ({n) —n/3 — 
r(n/3), to the anomalous scaling of the coarse-grained energy dissipation : 
(et(r)) ~ ^T(n./3)^ Anomalous scaling of isotropic structure functions is there- 
fore connected to the multifractal properties of the three dimensional energy 
dissipation field [5]. It should be noted however that the Refined Kolmogorov 
Hypothesis related inertial range scaling to scaling of dissipative quantities, 
and delicate issues connected to small distance expansions and fusion rules are 
being disregarded here [53,54]. At any rate, the relation presented by Eq.(18) 
did not advance the calculation of the scaling exponents beyond crude phe- 
nomenology. 

2.5 Persistence of Anisotropies 

A central issue of K41 phenomenology is the assumption of return-to-isotropy 
for smaller and smaller scales. Recently this assumption had been put to test 
in experiments and simulations [10,21,35,55]. A useful experimental set-up 
to test the return to isotropy is a homogeneous shear flow [6] where the large- 
scale mean- velocity has a linear profile: V = (V^)y,0,0). The shear is given 
by Sij — diVj = SiySjxVo. We thus have a homogeneous but anisotropic flow, 
close to the ideal case for studying the influence of large scale anisotropies 
on the small scale statistics. "Small scales" are defined here in comparison 
to the characteristic shear length, Ls = e^^^/S; for r -C we may expect 
that anisotropic fluctuations are sub-leading with respect to the isotropic ones. 
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The case r ^ L5 is of interest in situations where the shear is very intense, as 
very close to the walls in bounded flows. In such cases we expect a dramatic 
change from the K41 phenomenology [16,18,19]. Fortunately, it is not that 
difficult to design experiments or DNS possessing an almost perfect hnear 
profile with homogeneous shear [7, 16, 21, 55]. A popular way to measure small- 
scales anisotropics is to focus on the Re dependence of isotropic and anisotropic 
statistical observables built in terms of velocity gradients. For example, due 
to the symmetries of the mean flow, gradients of the stream-wise component 
in the shear direction, dyUx, may have a skewed distribution only due to the 
anisotropic fluctuations; they have a symmetric PDF in a perfectly isotropic 
flow. A natural measure of the residual anisotropy at small scales as a function 
of Re is the mixed generalized skewness based on gradients: 

M^'-^'\R>) f'''^^^Ui - (19) 



These objects vanish in isotropic ensembles. Of course, at flnite Reynolds 
numbers one expects that the large-scale anisotropy introduced by the shear 
still remains, even on the gradient scale. Therefore, the rate of decay of (19) as 
a function of Re is a quantitative indication of the rate of decay of anisotropy 
at small scales. In the next section we review Lumley's dimensional arguments 
[8] for anisotropic fluctuations, which predicts: 

^(2n+l)^^^^) _ y ^ (20) 



In fact, both numerical [20,21] (at low Reynolds numbers) and experimental 
tests (up to R\ ~ 1000) showed a clear disagreement with the dimensional 
prediction (20). For example in [7] the authors quote a decay in agreement 
with the prediction for M'^'^\Rx)^ an almost constant behavior as a function 
of Re for the fifth order, M^^^Rx) ~ 0(1) and an mcreasm^' behavior for the 
seventh order M^'^^Rx) ~ Rx^'^^ ! These results have cast a severe doubt on 
the fundamental assumption of the K41 theory. Similar results, with even more 
striking contradictions with the hypothesis of the return-to-isotropy, have been 
measured in the problem of passive scalar fluctuations, 9 — T — (T) , advected 
by an isotropic velocity fleld in the presence of a mean homogeneous scalar 
gradient, V(T) = {g,0,0). The equation of motions for the passive advected 
fleld in this case are: 

Both experimental and numerical data show a strong disagreement with the 
prediction that generalized skewness of temperature gradients becomes smaller 
upon increasing Reynolds and Peclet numbers [56-61]. 
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We will show below how the analysis based on S0(3) decomposition and its 
theoretical consequences settles this puzzle of strong persistence of anisotropies 
[35]. In fact, contrary to what appears, the K41 phenomenology with its as- 
sumption of return-to-isotropy and the above experimental results are not at 
all in contradiction (see section 6.2.1). 

2.6 Longitudinal and Transversal Isotropic Structure Functions 

Another debated issue concerning the K41 phenomenology and its multifrac- 
tal generalization (18) has to do with the observed discrepancies between 
the scaling properties of longitudinal, transversal, and mixed longitudinal- 
transversal structure functions in (supposedly) isotropic fully developed tur- 
bulence [30,62-65]. As previously stated, K41 theory, for isotropic flows, pre- 
dicts the same scaling behavior, in the limit of high Re, independent of the 
Cartesian components of the velocity increments in the structure functions. 
For a given order p — m-\- n only pref actors in (16) may depend on the par- 
ticular choice of n and m. Let us denote with (^('*'"*) the scahng exponent of 
the mixed structure function (16) made of n longitudinal increments and of 
m transversal increments in a isotropic ensemble: 

For p < 4 the issue does not exist; due to the incompressibility constraint 
all second and third order longitudinal or transversal structure functions have 
the same scaling in a isotropic ensemble. For p > 3 many experiments and 
numerical simulations found that (^("■''") < ) if n < n' and m> m' when 

n+m = n'+m'. It appears that with increasing m the scaling exponents reduce 
(the signal is more intermittent). The largest difference for a structure func- 
tion of order p is therefore achieved when we compare the purely longitudinal 
scaling C(P) 0) with the purely transversal scaling ({0,p). Other experimen- 
tal data suggests the possibility of a slow tendency of the longitudinal and 
transversal scaling exponents to coalesce for increasing Re [66-68]. 

We will argue below that the experimental measurements of different expo- 
nents stems from anisotropic corrections that affect differently the longitudinal 
and transverse components. In other words, by not removing the anisotropic 
contributions, one cannot expect pure power-law behavior. The situation is 
more complex for the analysis of data from DNS. There, one may implement 
highly isotropic forcing and boundary conditions, such that in most cases 
any residual anisotropic effects may safely be neglected even at moderate Re 
numbers. On the other hand, state-of-the-art numerical simulations are still 
strongly limited in the maximum Re achievable. Only very recently reliable 
data with high-statistics became available at resolution 1024^ [64], while most 
of the previous DNS where limited to lower resolutions. At resolution of 1024^ 
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one reaches a moderate R\ ~ 400, far below many experiments. Because of 
the consequent hmited extension of the inertial range, such DNS did not re- 
solve the puzzle of longitudinal vs transversal scaling. The numerical results 
oscillate between evidence for different scaling properties and for its oppo- 
site [62, 64, 69-71]. The issue is complicated by the fact that longitudinal and 
transversal structure functions possess different finite Re effects. For example, 
in [64] it was shown that structure functions of different order have differ- 
ent dependence on the viscous cut-off; this introduces ambiguity in defining 
a common inertial range where power law is expected. We thus propose that 
until high resolution isotropic measurements became available, all evidence for 
different scaling exponents for longitudinal and transverse structure functions 
should be considered with suspicion. 



2.7 Position Dependent Scaling Exponents 



In some inhomogeneous simulations and experiments it was claimed that the 
measured scaling exponents depended on the point of measurement within the 
flow domain [9,18,19,72-74]. If true, such flnding would deal a death blow 
to the idea of universality of the scaling exponents in turbulence. It should 
be stressed that in all the examples where such findings were reported the 
flow contained strong anisotropic and inhomogeneous components and/or the 
scaling range was not sufficient to actually present direct log-log plots for 
the structure functions vs. r. In some of these cases, scaling was extracted 
by using the method called "Extended Self Similarity" (ESS) [51]; the use of 
this method can be dangerous in presence of anisotropic and inhomogeneous 
effects. Whenever strong anisotropics are present one has to distinguish among 
two scaling ranges. At scales larger than the shear length, the energy cascade 
mechanism of the Kolmogorov theory is overwhelmed by shearing effects [16, 
19]. Only for scales smaller than the shear length the meaning of anisotropic 
corrections to the isotropic K41 scaling theory is well posed. We argue below 
that the reported position-dependent isotropic exponents in the latter case 
stem from anisotropic components which appear with different amplitudes 
at different points in the flow. The different "exponents" that were measured 
were not real exponents but the result of a crossover between the isotropic and 
anisotropic corrections. Once the data is projected onto the isotropic sector 
the leading exponents become position independent as expected. 
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3 Historical Review: Attempts at Anisotropy 

3.1 Bachelor's Approach 

The first systematic approacii to anisotropy in turbulence was suggested by 
G. K. Batchelor in 1946 [75]. Batchelor did not attempt to describe the most 
general form of anisotropy in turbulence, but instead confined himself to the 
easier case of axisymmetric turbulence. In axisymmetric turbulence the mean 
value of any function of the velocity field and its derivatives is invariant to all 
rotations of the axes in a given direction. Therefore, the anisotropy in axisym- 
metric turbulence is induced by a single direction in space. We denote this 
symmetry axis by the unit vector n. Being the easiest case of anisotropic tur- 
bulence, axisymmetric turbulence was the main model for studying anisotropy 
in subsequent years. 

Batchelor used the invariant theory in order to take the anisotropy vector n 
into account in the tensor representations. His method is simple: add the vector 
n to the list of vectors in Robertson's method. For example, suppose we wish 
to construct an axisymmetric representation of the second-order correlation 
function C"^(r) defined in (3). Then, just as in the isotropic case, we create 
a scalar function by contracting the two indices of C°'^{r) with two arbitrary 
vector a and b, with the difference that now we assume that the resultant 
scalar function depends on the unit vector n as well as on the other vectors 
r, a and h. We therefore look for an invariant representation of the scalar 
function C{r.in,a.ih).i which depends only on the different scalar products 
r • r, r • n, r ■ a and the various determinants [rna], [rnb] For some 
reason, Batchelor decided to ignore the skew-symmetric parts and considered 
only the scalar products. Using the fact that C(r, n, a, h) is a bilinear function 
of a and 6, Batchelor found that 

where A, B,C, D, E are functions of the amplitude r and of the scalar product 
r • n. Notice that in this expansion the number of unknowns has grown from 
two to five with respect to the isotropic expansion. It would have been nine, 
had we taken the skew-symmetric parts into account. Indeed, a prominent 
characteristic of anisotropic representations is that they are far more complex 
than their isotropic counterparts. 

Using this sort of representations, Batchelor was able to generalize Karman- 
Howarth results to the case of axisymmetric turbulence. That is, after repre- 
senting the second- and third-order correlation-functions in terms of few scalar 
functions, Batchelor used the solenoidal condition and the Navier-Stokes equa- 
tions to derive some linear differential relations among them. 
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3.2 Chandrasekhar and Lindborg's Approaches 

A somewhat more elegant approach to axisymmetric turbulence was offered a 
few years later by S. Chandrasekhar [76]. Chandrasekhar 's treatment is sim- 
ilar to Batchelor's in following Robertson's work [39]. Chandrasekhar took 
advantage of the skew-symmetric tensor e^"^ for creating a representation of 
solenoidal axisymmetric tensors. He noticed that the curl of an axisymmetric 
tensor automatically satisfies the solenoidal condition. Therefore, by repre- 
senting the second- and third- correlation-functions as a curl of auxiliary ten- 
sors, Chandrasekhar automatically solved the solenoidal equations, and was 
left with the dynamical equations (which are derived from the Navier-Stokes 
equation) only. Chandrasekhar's dynamical equations are considerably sim- 
pler than those of Batchelor. Nevertheless, they are still very comphcated and 
this, perhaps, explains why there was no serious attempt to continue Chan- 
drasekhar's work in subsequent years. 

In 1995 there was another attempt to formulate the kinematics of homogeneous 
axisymmetric turbulence in [77] . The representation in this paper was "exper- 
imentally oriented" , in the sense that the scalar functions that are used can be 
measured directly in experiment. To accomplish that, one defines two auxiliary 
unit vectors (that were also used in [78]): ei{r) = n x r /\n x r\ , e2{r) = 
ei(r) X n . The triplet (n, ei, 62) is an orthonormal basis of M.^. But since it 
is made out of r and n using invariant operations (i.e., vectorial products), 
it is invariant to simultaneous rotations of r and n, and thus it is invariant 
to rotations of r alone around n (because in these rotations n remains fixed 
anyhow). Therefore any tensor that is built from these unit vectors and the 
products r ■ r, r ■ n is necessarily an axisymmetric tensor. For example, in 
order to represent the second-order correlation- function C"^(r), one writes 

C"^(r) = i?in"n^ + i^aCsC^ + R^e^e^ + R^[n''e^ + n^e^] 
+ gi[n"ef + nf^el] + g2[ef ef + efe?] , 

where the six scalar functions . . ., i?4 and Qi, Q2 are functions of fi = r - n 
and p = \r X n\. Notice that this representation takes into account the skew- 
symmetric part of C"'^(r) using the scalar functions Qi and Q2. 

The major advantage in this representation is that the scalar functions have 
an immediate interpretation in terms of measurable quantities. With respect 
to the example above, if {u, v, w) are the velocity components in the direction 
of (n, 62, ei) respectively, then due to the orthonormality of triplet (n, 62, ei), 
we get 

Ri = {u{x)u{x + r)) , R2 = {v{x)v{x + r)) , 
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= {w{x)w{x + r)) , R4 = {u{x)v{x + r)) , 
Qi ^ {u{x)w{x + r)) , Q2 ^ {v{x)w{x + r)) . 

Next one uses the solenoidal condition to derive linear differential relations 
between the scalar functions. One can also consider the triple correlation func- 
tion and the velocity-pressure correlation function in the very same method. 
This way one derives a representation for the dynamical equation of C"^(r). 
Mathematically speaking, this representation is no better than Batchelor's rep- 
resentation, and may even be considered worse than Chandrasekhar's. This is 
because there is no reduction in the number of scalar functions and there is 
no simplification in the resulting equations. The motivation for this represen- 
tation is experimental without compelling physical or mathematical contents. 



3.3 Case Specific Approaches 



The above mentioned works can be viewed as systematic attempts to deal with 
the problem of anisotropic turbulence, where a general method to describe 
the anisotropic (or, more exactly, axisymmetric) quantities in turbulence is 
suggested. In that respect they differ from most research that followed on the 
subject, which was usually confined to a particular model or a specific problem 
related to anisotropy. 



3.3.1 Temporal Return to Isotropy 

One such problem is the temporal return to isotropy in which one tries to 
understand the mechanisms that drive a decaying turbulence which is ini- 
tially anisotropic into being statistically isotropic. As in the works of Batch- 
elor and Chandrasekhar, the statistics was usually assumed to be spatially 
homogeneous to simplify the problem. The theoretical attempts to explain 
this phenomenon can be roughly divided into two groups. The first group, 
which was initiated by Rotta [79] in 1951, consists of attempts to model the 
decay of anisotropy using one-point closures. In that framework, one usually 
considers the dynamical equation of the Reynolds stress, which is the same- 
point correlation-function of the velocity field C°'^{t) = (u°'(x,t)u^(x,t)^ . 
In a homogeneous, decaying turbulence, this correlation obeys the following 
equation 



with 

60"^ = (u^'d^p) + (u^dy) + 2 



(21) 
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Notice that C ■ is a second-rank r-independent tensor that contains an isotropic 
part (which is its trace) and an anisotropic (traceless) part. This explains the 
motivation behind the definition of 0"^, which is to capture the anisotropic 
part of the decay rate of C"'^. To solve Eq. (21), one must model the 0"^^ 
tensor. Usually this has been done on a phenomenological basis. A systematic 
treatment to this problem was ofi'crcd by in [80]. In that paper, the authors 
suggested that 0"^^ should depend on time implicitly through C"^, e and u. 
Additionally, since 0"^ is a dimensionless tensor, it must depend only on di- 
mensionless parameters. There are six such independent dimensionless quanti- 
ties. The authors chose to represent them in a way which isolates the property 
of anisotropy from other properties, and form the tensor 



The tensor h"'^ is proportional to the anisotropic, traceless part of C"^ and 
hence contains five independent components. It is often denoted as the "Reynolds 
stress anisotropy" or simply as the "anisotropy tensor". This tensor has be- 
come a central measure for anisotropy in turbulence and has been used ex- 
tensively in experimental and numerical analysis of anisotropy in turbulence. 
The sixth component was defined to be proportional to the isotropic part of 
C"^ (the energy) by i?; = With these dimensionless quantities, (f)"'^ can 
be written as = (j)"^{b,Ri) . They further simplified that expression by 
noticing that if 0"^^ depends solely on h and Ri then it must depend on them 
in an isotropic manner, since any anisotropic dependence necessarily means 
that (j)"'^ also depends on the boundary conditions. To represent this isotropic 
dependence explicitly, they used the invariant theory [81] and introduced the 
second- and third- principal-invariants of the traceless tensor b: 



According to the invariant theory, can be most generally written as 



and the problem is reduced to determining the functional form of f3{II, III, Ri) 
and 7(//, //, Ri). Based on this formalism, there have been many attempts to 
model the functions I3{II , III , Ri) and ■y{II,III,Ri) to match experimental 
results [82]. For example, Rotta's model is considered a linear model for the 
anisotropy decay because he used 




// = Tr[b^] , III = Tr[b^] . 




(22) 



P{II, III, Rj) = Ci « 3.0 , 



j{II,III,Rj) = . 



Consequently, the decay of the anisotropy tensor is given by 



at6"^ = -(e/g2)(c'i-2)6"^ 
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which is a Uncar equation in 6"^, provided that the isotropic quantities e, are 
independent of If-^ . This sort of equation predicts that \f^(t) is proportional 
to 6"'5(t = 0) and that every component of the tensor decays at the same rate. 
Experimentally however linearity is not supported. For example, it has been 
observed experimentally that the return to isotropy is relatively rapid, at least 
at the beginning of the process, when the invariant /// is negative, whereas 
the return to isotropy is fairly slow in the case where the invariant is positive 
[83]. 

The other line of research that was used to study the problem of the return- 
to-isotropy consists of attempts to model the decay with two-point closures. In 
these models one considers the different correlation functions of the velocity 
field across a separation vector r, instead of using the same-point correlations 
as in the one-point closures. The mathematical structure here is usually much 
more complicated than that of the one-point closures, but in return, the two- 
points models often provide a deeper understanding of the physics involved. 

An important example for such a model is given by [78] . In this work the author 
used the direct-interaction approximation (DIA) [84] to study the decay of an 
axisymmetric turbulence into isotropy. The DIA is a well-known truncation 
of the renormalized perturbation theory for turbulence. The perturbation is 
done in the interaction strength parameter (which is set to unity in the end), 
and is truncated at the second-order - i.e., at the direct-interaction terms. 
The calculations are done in Fourier space. They result in two coupled equa- 
tions for the time evolution of the two-time, second-order correlation-function 
C"^(fc,t,t') = (u'^{k,t)u^{-k,t')) and the response function G'^^{k,t,t'). The 
latter is defined as the average of the change in the velocity field at time t 
as a result of an infinitesimal change in the forcing at time t' . The equations 
for C"^{k, t, t') and G°'^{k, t, t') determine their time evolution. They are non- 
linear, nonlocal integro-differential equations and are therefore very hard to 
deal with. In the isotropic case the equations can be considerably simplified 
by noting that both C'^^{k,t,t') and G°'^{k,t,t') must satisfy the solenoidal 
condition, which in Fourier space means that both tensors must vanish once 
we contract any of their indices with the vector k. It is easy to see that under 
such a condition, C'^f^{k,t,t') and G''f^{k,t,t') can be represented in terms of 
one scalar function: 

G'"^{k, t, f) = c{k, t, f)D'"^{k) , G'"'^(fc, t, f) = g{k, t, t')D''^{k) , (23) 



where 

D"^(fe) = 5^*^ - k'^k^ . (24) 
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When turning to the axisymmctric case, the representation of the tensors 
become much more complex. Instead of one scalar function for each tensor, 
two functions must be used corresponding to the two scalar functions that 
where used by Batchelor [75] and Chandrasekhar [76] to describe the second- 
order correlation-function in real space. One uses the separating vector k and 
the anisotropy unit- vector n to create two unit vectors which are orthogonal 
to k: 

ei{k) =^ fc X n/\k x n\ , e2{k) =^ fe x (fc x n)/\k x (fc x n)| . (25) 

With these vectors, C°'^{k,t,t') was written as 

C°^(fe, t, t') = ci(fc, t, Oe?ef + C2(fc, t, t')e^e^ , (26) 

and G''^{k,t,t') was written in a similar way Tising the functions gi{k, t, t') and 
g2{k, t, t'). To parameterize the angular dependence of the scalar functions, one 
expands the four scalar functions in terms of Legendre polynomials, 

Ci{k,t,t') ^Y.Ci^j{k,t,t')Pj{k ■ n) , gi{k,t,t') ^J29i,j{k,t,t')Pj{k-n) , 
j j 

obtaining an infinite set of coupled equations for the infinite set of functions 

Cij{k, t, t') and gij{k, t, t'). These equations were solved numerically after trun- 
cating all the j > part of the expansion. Doing so, one finds the time evolution 
of Cifi{k,t,t') and gi^o{k,t,t') and connects them to the physical observables 
of the one-point closures, such as Rotta's constant. The conclusions (partially 
numerical and partially obtained after a long series of uncontrolled approx- 
imations) were that the return-to-isotropy is much stronger at small scales 
(large k) and that in some classes of initial conditions, the return-to-isotropy 
is indeed a linear phenomenon. 

This calculation was soon revised in [85]. In this paper, the authors compared 
the DIA calculation to the results of a numerical simulation of a homogeneous 
and axisymmetric turbulence. This time, however, the Legendre polynomials 
expansion in the DIA calculation was extended to include also the j — 2 
components. Their conclusions were that the DIA calculation that included the 
j = 2 parts were in a good agreement with the numerical simulation, especially 
the j = parts, provided that the initial anisotropy was small. Additionally, 
the authors found that the previous calculation [78], which considered the 
j — part only, was quite inadequate to describe the process of return-to- 
isotropy - even in the case of weak anisotropy. 

Other attempts to study the problem of return-to-isotropy in an axisymmetric 
turbulence used two-point closures. For example, eddy-damped quasinormal 
Markovian approximation (EDQNM) [86,87] have been used in [88]. In this 
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closure scheme, one approximates the fourth order cumulants of the velocity 
field by a linear damping term of the third-order correlation-function of the ve- 
locity field (the eddy dumping). Additionally, a "Markov" assumption is used 
that allows one to integrate the history integral in the equations and retain an 
equation for the same time second-order correlation-function C"^^(fc,t). As in 
the DIA model, this equation is formulated in Fourier space and is both non- 
linear and nonlocal. To parameterize the axisymmetric correlation function in 
Fourier space, the authors used the following representation 

C"*^(fc,t) = ^[D"'^{k)E{k,fi,t) + n''n''D''^{k)Df^''{k)F{k,fi,t)] , 

where here ji = k ■ n and D°'^{k) is defined in Eq. (24). This representation 
incorporates the solenoidality condition and is more elegant than (26) in the 

sense that the isotropic case is very easily recovered once we set F{k, fi, t) = 
and let E{k,fi,t) become /x- independent. Plugging this expansion to the 
dynamical equation of C^^ik^t), the authors obtained two coupled equa- 
tions ioT E{k, n,t) and F{k, fi,t) which they solved numerically for "medium" 
and "strong" anisotropic initial conditions. Their results indicate that in the 
medium anisotropy cases, Rotta's constant approaches a constant of the or- 
der of unity, qualitatively agreeing with Rotta's model and with the results 
[78]. On the other hand, in the strong anisotropy case this constant does not 
show any saturation, indicating the failure of Rotta's model. Additionally, 
their results support the idea that the decay isotropy strengthens when the 
III invariant is negative. 

3.4 Dimensional Analysis in the Presence of Strong Shear 

An important discussion of the effects of strong shear on the energy spectrum 
was presented in [8]. In this paper the author included anisotropic corrections 
in to the K41 framework extending the phenomenological dimensional rea- 
sonings leading to (14). He considered the dependence on anisotropic mean 
observables, like the large-scale Shear proportional to the large-scale mean 
gradient: S oc d{V): 

5W(r) = (er)t/(")(^,!^,5). 

By further assuming that anisotropic corrections are "small" and analytic in 
the intensity of the shear S, he proposed the following form for the anisotropic 
correction to the isotropic two-point longitudinal structure functions, in the 
inertial range [8]: 

^(2)(r.) ^ C(2)(er)i + D^''\r)Srl (27) 
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where the coefficient D^'^^{r) takes into account the dependence on the direc- 
tion r in the anisotropic term. The counterpart of (27) for the spectrum and 
co-spectrum in Fourier space is: 

{k\,{k)ui{-k)) - k-l{5u - ^)+Auk-l (28) 

where the first term on the RHS is the isotropic K41 scahng and the second 
term is the anisotropic contribution with An being a traceless matrix depend- 
ing on the details of the large scale shear. 

In the past, most of the measurements of the anisotropic contributions to 
S^'^\r) concentrated on the Fourier representation (28), [9,11,12,89]. In [9] 
the authors showed that the prediction (27 28) is well verified in a wind tun- 
nel flow. Later, many other experiments have confirmed this result in different 
experimental situations (see for example the recent results for an homogeneous 
shear in [7]). Only recently, a more extensive study of anisotropics has been 
carried out, considering also higher order statistical objects [7, 10, 29, 30]. The 
situation became immediately less clear: the prediction (28) is not the end 
of the story (see below the section on anomalous scaling for anisotropic fluc- 
tuations). We show later that in the jargon of the S0(3) decomposition the 
anisotropic part of a spherically averaged and solenoidal second-rank tensor is 
made from j = 2 contribution only, for this reason the dimensional analysis is 
often viewed as predicting a 4/3 exponent for the j = 2 sector of the second- 
order structure function. This result was later derived by several authors in 
terms of Clebsch variables, but again by dimensional reasoning [90-92]. An- 
other, more systematic attempt to derive the scaling behavior of the second 
order structure function in a weakly anisotropic turbulent flow was presented 
in [93] within a variable scale mean field theory. In that paper the authors 
reached the conclusion that all anisotropic contribution to the second order 
structure function must scale ~ r^/^. To reach this result the authors had to 
simplify the tensorial structure of the equations for the second order correla- 
tion functions; we argue below that this uncontrolled simplification biased the 
estimate of the anisotropic exponents. 



4 The Modern Approach to Anisotropy 

In the past 10 years, the subject of anomalous scaling in turbulence has gained 
a great deal of attention, as it became more and more accepted that in the 
infinite Re limit, the scaling exponents of the structure functions in the inertial 
range do not conform with the classical prediction of the Kolmogorov theory. 
The numerical values of these exponents, as well as the physical mechanism 
which is responsible for the anomalous scaling, have been the target of an 
extensive experimental, numerical, and theoretical research. 
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On the theoretical side, important progress was made by studying Kraichnan's 
model of passive scalar advcction [94]. This model describes the advection of 
a passive scalar field by a synthetic, solenoidal velocity field with a Gaus- 
sian, white-in-time statistics. The linearity of the equations for the passive 
scalar field and the white-in-time statistics of the velocity field make it pos- 
sible to write down a closed set of equations for the same-time correlation 
functions of the passive scalar [94]. In [95,96], it was shown that the solu- 
tion of these equations can lead to anomalous scaling. The key point is that 
the homogeneous solutions of these equations are scale invariant with non- 
trivial anomalous scaling exponents, which are different from the dimensional 
scaling exponents that characterize the inhomogeneous, "forced" solutions. 
Being usually smaller than the dimensional scaling exponent, the anomalous 
exponents dominate the small scales statistics of the passive scalar field. The 
homogeneous solutions are commonly referred to as "zero modes" , and have 
been calculated to first order perturbatively in refs. [95, 96] for the fourth or- 
der structure function and for all even structure functions in ref. [97]. Exact 
computer assisted calculations of the exponents of the third order structure 
functions were presented in [98] . Besides suggesting an elegant mechanism for 
anomalous scaling, Kraichnan's model also provided an example in which the 
scaling of the anisotropic parts of structure functions is different from the 
isotropic scaling. In the paper [99] it was shown how such a thing can happen, 
by expanding the second-order structure function of the passive field in terms 
of spherical harmonics Yj^mi'f')- It was found that this expansion leads to a 
set of decoupled j-dependent equations for the expansion prefactors. These 
equations can be easily solved by a power law whose exponent is an increas- 
ing function of j. These exponents are universal in the sense that they are 
independent of the forcing and boundary conditions. 

The authors in [99] also noticed that the fact that the anisotropic exponents 
are higher than the isotropic exponent neatly explains the isotropization of 
the statistics as smaller and smaller scales are probed. Based on this example, 
it was suggested in [100] that a similar mechanism may exist in a Navier- 
Stokes turbulence. The authors expanded the second-order structure function 
in terms of spherical harmonics 

S^'\r) = Y.Sl'J,{r)Yar), (29) 

and argued that in the case of weak anisotropy, one can linearize the equations 
for the anisotropic corrections of the second-order structure function around 
the isotropic solution. In such the kernel of the linearized equation 

is invariant under rotations (isotropic), and as a result the equations for the 
different (j, m) components decouple, and are m-independent - much as in 
the case of the second-order structure function in Kraichnan's model. In a 
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scale-invariant situation, this leads to anisotropic, j-dependent exponents 



If one accepts that homogeneous turbulence enjoys universal statistics in the 
inertial range, then the kernel of the above linearized equation is universal, and 
consequently so are the anisotropic scaling exponents Cj^^ ■ All of these state- 
ments could not have been proved rigorously (and still haven't been proved 
rigorously), yet they offered a new approach to understanding anisotropy in 
turbulence, an approach that is explored in the rest of this review. 

4.1 Mathematical Framework 

Experiments in fluid turbulence are usually limited to the measurement of the 
velocity field at one single spatial point as a function of time. This situation 
has begun to improve recently, but still much of the analysis of the statisti- 
cal properties of Navier-Stokes turbulence is influenced by this tradition: the 
Taylor hypothesis [101] is used to justify the identification of velocity differ- 
ences at different times with differences of longitudinal velocity components 
across a spatial length scale r. Most of the available statistical information 
is therefore about properties of longitudinal two-point differences of the Eu- 
lerian velocity field and their moments. Recent research [53] has pointed out 
the advantages of considering not only the longitudinal structure functions, 
but tensorial multi-point correlations of velocity field differences 



F"i-""(a;i, a;;, ii; ; . . . ; Xn, tn) = «(a;i, x[, ti) . . . w^{xn, tn))(30) 



where all the coordinates are distinct. When the coordinates fuse to yield 
time-independent structure functions depending on one separation only, these 
are the so-called tensorial structure functions, denoted as 



ga,...an ^ ^[^ai (j. + y.) _ (3.)] • • • [-U°" {x + v) - w"" (x)]) . (31) 



Needless to say, the tensorial information is partially lost in the usual mea- 
surements conducted at a single point. One of the main stresses of the present 
review is that keeping as much of tensorial information as possible can help 
significantly in disentangling different scaling contributions to the statisti- 
cal objects. Especially when anisotropy implies different tensorial components 




w{x, x', t) = u{x', t) — u{x, t), 



given by 
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with possible different scaling exponents characterizing them, careful control 
of the various contributions is called for. 



To understand why irreducible representations of the symmetry group may 
have an important role in determining the form of correlation functions, we 
need to discuss the equations of motion which they satisfy. We shall show that 
the isotropy of the Navier-Stokes equation and the incompressibility condition 
implies the isotropy of the hierarchical equations which the correlation func- 
tions satisfy. We will use this symmetry to show that every component of the 
general solution with a definite behavior under rotations (i.e., components of 
a definite irreducible representation of the SO (3) group) has to satisfy these 
equations by itself - independently of components with different behavior un- 
der rotations. This "foliation" of the hierarchical equations may possibly lead 
to different scaling exponents for each component of the general solution which 
belong to a different SO{3) irreducible representation. 

4.2 Tensorial Correlation Functions and S0{3) Irreducible Representations: 
General Theory 

The physical objects that we deal with are the moments of the velocity field 
at different space-time locations. In this section we follow Ref. [23] which 
suggests a way of decomposing these objects into components with a definite 
behavior under rotations [23]. It will follow that components with different 
behavior under rotation are subject to different dynamical equations, and 
therefore, possibly, scale differently. Essentially, we are about to describe the 
tensorial generalization of the well-known procedure of decomposing a scalar 
function ^(r) into components of different irreducible representations using 
the spherical harmonics: 



4.2.1 Formal Definition 

Consider the correlation function F*^"^ of Eq. (30). This n-rank tensor is a 
function of 2n spatial variables and n temporal variables. It transforms as a 
tensor field: if F^"^ is measured in two frames / and / which are connected by 
the spatial transformation (say, a rotation) x°' — A°'^x^ then, the measured 
quantities in each frame will be connected by the relation: 




(32) 



F"(xi,x;,ti;...;x„,x;,tO=A"i^i 
= A"i/?i . . .A°"/3"F^i-^"(A-^xi, A-^x;,ti; 



A"n/3n^/3l.../3n(^Xl,X;,ti;...; 

•A"'x A"'x' f ) 
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We see that as we move from one frame to another, the functional form of the 
tensor field changes. We want to classify the different tensor fields according 
to the change in their functional form as we make that move. We can omit the 
time variables from our discussion since under rotation they merely serve as 
parameters. We thus define T{{xi}) = F{{xi}, {U — 0}). Consider coordinate 
transformations which are pure rotations. For such transformations we may 
simplify the discussion further by separating the dependence on the amplitude 
of Xj from the dependence on the directionality of Xj: 

rpai...a„(-^ Y ^ Tion-.-anf^ ~ .-A- N 

J. V I' • • • ■> ■'^J — V'*'!) ■ ■ ■ ) -^pt ■'^l: • • • ■> ^p)-i 

where here we have p < n, i.e we consider also the possibility that n — p 
spatial locations in (33) coincide. For pure rotations we may treat the am- 
plitudes parameters: the transformations properties of 
under rotation are determined only by the dependence of T"^ ' "" on the unit 
vectors xi,...,Xp. Accordingly it seems worthwhile to discuss tensor fields 
which are functions of the unit vectors only. Notice that in the scalar case 
we follow the same procedure: by restricting our attention to scalar functions 
that depend only on the unit vector x, we construct the spherical harmonics. 
These functions are defined such that each one of them has unique transfor- 
mation properties under rotations. We then represent the most general scalar 
function as a linear combination of the spherical harmonics with x-dependent 
coefficients, see Eq. (32). 

The classification of the tensor fields T"i'"""(xi, . . . , Xp) according to their 
functional change under rotations follows immediately from group representa- 
tion theory [102, 103]. But in order to demonstrate that, we must first make 
some formal definitions. We define to be the space of all smooth tensor fields 
of rank n which depend on p unit vectors. This is obviously a linear space of 
infinite dimension. With each rotation A G 5*0(3), we may now associate a 
linear transformation 0\ on that space via the relation (33): 

[CaT]"^-"" (Xl, . . . , Xp) = A^l^l • • • A"n/3n7^/3l-/3n(^-l-j^^^ _ _ _ ^ J^l^^^ 

Using this definition, it is easy to see that the set of linear operators 0\ 
furnishes a representation of the rotation group SO (3) since they satisfy the 
relations: 

CaiC'az = C'AiAz, = Oa-1. 

General group theoretical considerations imply that it is possible to decompose 
Sp into subspaces which are invariant to the action of all the group operators 
0\. Moreover, we can choose these subspaces to be irreducible in the sense that 
they will not contain any invariant subspace themselves (excluding themselves 
and the trivial subspace of the zero tensor field). For the 5*0(3) group each 
of these subspaces is conventionally characterized by an integer j = 0, 1, 2, . . . 
and is of dimension 2j + 1 [102, 103]. It should be noted that unlike the scalar 
case, in the general space <S", there might be more than one subspace for each 
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given value of j. We therefore use the index q to distinguish subspaces with 
the same j. For each irreducible subspace {q,j) we can now choose a basis 
with 2j + 1 components labeled by the index m: 

BZm'^^i^i, . . . , Xp) ; m = -J, . . . , +j. 

In each subspace {q,j), the group operators Oa furnish a 2j + 1 dimensional 
irreducible representation of 5*0(3). Using the basis i?°j^''°"(xi, . . . ,Xp), we 

can represent each operator 0\ as a (2j + 1) x (2j + 1) matrix dI^,^(A) via 
the relation: 

[OaB];;-''" (xi, . . . , X,) = A-^^ . ■ • A-/^"S5^''(7t^xi, . . . , 

^ E ^HL(A)^,%'""(xi,...,ip). 

m'=-j 

It is conventional to choose the basis B^j^ such that the matrices Dl^,j^{(f)), 
that correspond to rotations of radians around the 3 axis, will be diagonal, 
and given by: D^^,^{(f)) — Smm'^^^'^- The space possesses a natural inner- 
product: 

(T, U)= |cixi . . . dxp • T"— (xi . . . %)ga,p, • • • 9.^.^^" "'" (^i • ■ • 

where is the 3-dimensional Euclidean metric tensor. By definition, the 
rotation matrices h."'^ preserve this metric, and therefore it is easy to see that 
for each A e SO{?>) we get: 

(C»aT,0aU) = (T,U) 

so that, Oa are unitary operators. If we now choose the basis 'Qq,jm fo be or- 
thonormal with respect to the inner-product defined above, then the matrices 
-^m'm(-A) ^i^^ be unitary. 

Finally, consider isotropic tensor fields. An isotropic tensor field is a tensor 
field which preserves its functional form under any arbitrary rotation of the 
coordinate system. In other words, it is a tensor field which is invariant to 
the action of all operators Oa- The one dimensional subspace spanned by this 
tensor-field is therefore invariant under all operators Oa and therefore it must 
be a j = subspace. 

Once the basis B^jm has been selected, we may expand any arbitrary tensor 
field F"^ ' ""(xi, . . . ,Xp) in this basis. As mentioned above, for each fixed set 
of amplitudes xi, . . . , Xp, we can regard the tensor field F"i---""(xi, . . . , Xp) as 
a tensor field which depends only on the unit vectors Xi, . . . ,Xp, and hence 
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belongs to Sp. We can therefore expand it in terms of the basis tensor fields 
^g,jm with coefficients that depend on xi, . . . ,Xp: 

F ^ "(Xi, . . . , Xp) = ^ }jm(^l) ■ ■ ■ ) •^p)BqJin (Xl, ■ ■ ■ , Xp) . (34) 

The goal of the following sections is to demonstrate the utility of such ex- 
pansions for the study of the scaling properties of the correlation functions. 
For the important case of tensorial structure functions (31) the basis function 
depend on one spatial vector only r, and we can expand 

S^-\r)=j:S^Ur)B^:Uf-)- (35) 



4.2.2 Construction of the Basis Tensors 

4.2.2.1 The Clebsch-Gordan machinery. A straightforward (although 
somewhat impractical) way to construct the basis tensors ^qjm is to use the 
well-known Clebsch-Gordan machinery. In this approach we consider the 
space as a direct product space of n 3-dimensional Euclidean vector spaces 
with p infinite dimensional spaces of single- variable continuous functions on 
the unit sphere. In other words, we notice that Sp is given by: 



1 n r n 
So] ® 5° 



and therefore every tensor T'°i- "^"(xi . . . Xp) can be represented as a linear 
combination of tensors of the form: 

(xi)-...-(/Pp(Xp). 

where v"' are constant Euclidean vectors and V9i(xj) are continuous functions 
over the unit sphere. The 3-dimcnsional Euclidean vector space, Sq, contains 
exactly one irreducible representation of 5*0(3) - the j = I representation - 
while S'l, the space of continuous functions <S> over the unit sphere, contains 
every irreducible representation exactly once. The statement that Sp is a direct 
product space may now be written in a group representation notation as: 



5^=i®i®...®i®(o®ie2...)(8)...(oeie2...) 

We can now choose an appropriate basis for each space in the product: 
• For the 3-dimensional Euclidean space we may choose: 

ei = (1, 0) , eo = (0, 0, 1) , e_i = (1, -z, 0) 
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• For the space of continuous functions over the unit sphere we may choose 
the well-known spherical harmonic functions. 

Once these bases have been chosen, we can construct a direct-product basis 
for S^: 

^h.'.'.Z{hlJ.l)-{lpHp) (•''■1' ■ ■ ■ ' = ■ • • • ■ ' ■ • ■ • ■ ^p,/ij,(^p)- 

The unitary matrix that connects the Ejj...j^(ij^^) ..(ip^p) basis to the '^q,jm ba- 
sis can be calculated using the appropriate Clcbsch-Gordan coefficients. The 
calculation is straightforward but very long and tedious. However, the above 
analysis enables us to count and classify the different irreducible representa- 
tions of a given j. By using the standard rules of angular-momentum addition: 

s®l^\s-l\® ...®{s + l) 

we can count the number of irreducible representations of a given j. For ex- 
ample, consider the space Si of second-rank tensors with one variable over the 
unit sphere. Using the angular-momentum addition rules we see: 

= 1 ® 1 (0 © 1 ® 2 © 3 © . . .) (36) 

= (0©1©2)®(0©1©2©3©...) 

= (3 X 0) © (7 X 1) © (9 X 2) © (9 X 3) © .. . 

We see that there are exactly three j = representations, seven j = 1 repre- 
sentations and 9 representations for each j > 1. It can be further argued that 
the symmetry properties of the basis tensors with respect to their indices come 
from the 1(8)1 = 0©1©2 part of the direct product (36) . Therefore, out of the 
9 irreducible representation of j > 1, 5 will be symmetric and traceless, 3 will 
be anti-symmetric and 1 will be trace-full and diagonal. Similarly, the parity 
of the resulting tensors (with respect to the single variable) can be calculated. 

Once we know how many irreducible representations of each j are found in 
iSp, we can construct them "by hand", in some other, more practical method 
which will be demonstrated in the following subsection. 



4.2.2.2 Alternative derivation of the basis functions. The method 
found most useful in application is based on the simple idea that contractions 
with r", 5"^, e"^''' and differentiation with respect to r" are all isotropic oper- 
ations [23]. Isotropic in the sense that the resulting expression will have the 
same transformation properties under rotation as the expression we started 
with. The proof of the last statement follows directly from the transformation 
properties of r", 5"^, e"'^'^ . 
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The construction of all ^qjm that belongs to Si now becomes a rather trivial 
task. We begin by defining a scalar tensor field with a definite j. m. An obvious 
choice will be the well-known spherical harmonics Yjjn[r), but a better one will 
be: 



The reason that we prefer %m(r) to Yj^{r), is that $jm(r) is polynomial in r 
(while Yjm{r) is polynomial in f) and therefore it is easier to differentiate it 
with respect to r. Once we have defined $j>„(r), we can construct the ^qjm by 
"adding indices" to $j^(r) using the isotropic operations mentioned above. 
For example, we may now construct: 

• r~^~^r"^jjn{r): etc... 

Notice that we should always multiply the resulting expression with an ap- 
propriate power of r, in order to make it r-independent, and thus a function 
of f only. 

The crucial role of the Clebsch-Gordan analysis is to tell us how many repre- 
sentations from each type we should come up with. First, it tells us the highest 
power of r in each representation, and then it can also give us the symmetry 
properties of ^qjm with respect to their indices. For example, consider the 
irreducible representations of Sf - second rank tensors which depend on one 
unit vector r. The Clebsch-Gordan analysis shows us that this space contains 
the following irreducible representations spelled out in (36). That is, for each 
j > 1 we are going to have 9 irreducible representations. The indices symme- 
try of the tensor comes from the Sq^Sq = 1(8)1 = 0®1®2 part of the direct 
product. This is a direct product of two Euclidean spaces, so its a second rank 
constant tensor. We can mark the representation number in this space with 
the letter s, and the representation number of the Si = 0©1©2©3©... 
space with the letter /. This way each representation in Sf of a given j will 
have two additional numbers (s,/), which actually serve as the index q that 
distinguishes irreducible representations of the same j. The s index will de- 
termine the indices symmetry of the tensor, while the / index will determine 
the highest power of r in the tensor. If we now recall that in the space of 
constant second-rank tensors, Sq^Sq = 0©1©2, the s = 0, 2 representations 
are symmetric while the s = 1 representation is anti-symmetric, we can easily 
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construct the -B^ f^: 
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*im(r), 



It should be stressed that these are not exactly the same one we would 



have gotten from the Clebsch-Gordan machinery. For example, they are not 
orthogonal among themselves for the same values of j, m (although, they are 
orthogonal for different values of j or m). Nevertheless, they are linearly in- 
dependent and thus span a given (j, m) sector in the Sf space. The set of 
eigenf unction, B^^^, can be further classified in terms of its properties under 
permutation of tensorial indices, af3 and in terms of their parity properties, 
i.e. how do they transform under the r —r operation. Taking in to account 
both properties we may distinguish: 

Subset I Symmetric in a, (3 and with parity (—1)-' : 

Bfjmir), B^Smir), <^(r), <„(f ) 
Subset II Symmetric to a,P exchange and with parity (— 1)-'+-'^: 

<,n(r),<„(r) 

Subset III Antisymmetric to a,P exchange and with parity (—1)-'''"^: 

B2U^),Bfj^{r) 
Subset IV Antisymmetric to a, (3 exchange and with parity (—1)-': 

BsUr) 

The reader may find more details on the algebra of SO (3) decomposition of 
second order tensor in the Appendix A. 
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4-3 The Isotropy of the Hierarchy of Equations and its Consequences 



In this section we follow Ref. [23] in deriving equations of motion for the sta- 
tistical averages of the velocity and pressure fields differences. We start from 
the Navier-Stokes equations, and show that its isotropy implies the isotropy 
of the equations for the statistical objects. Finally, we demonstrate the folia- 
tion of these equations to different sectors of j, m. Consider the Navier-Stokes 
equations (2) in a bounded domain Q. In principle these equations can be the 
basis for deriving infinite linear hierarchy of equations for the Eulerian corre- 
lation functions and to study its properties under rotation. Unfortunately the 
relevant dynamical time scales are revealed only when the effect of sweeping is 
removed. Therefore we choose to work here with the transformation proposed 
in [104] in which the fiow is observed from the point of view of one specific 
fiuid particle which is located at Xq at time to. Let p(xo, to\t) be the particle's 
translation at time t: 



We then redefine the velocity and pressure fields to be those seen from an 
inertial frame whose origin sits at the current particle's position: 



7r(xo, to\x, t) = p[x + p(xo, to\t),t] . 
Next, we define the differences of these fields: 

>V"(xo, tojx, x', t) = v'^i^o, to|x, t) - v'^i^o, to|x', t) , 

n(ro,to|x,x',t) = 7r(xo,to|x,t) - 7r(xo, to|x', t) . 
A straightforward calculation shows that the dynamical equations for W are: 



dtW^i^, x', t) = - {d^ + ffj n(xo, to|x, x', t) + u [d' + d") W"(xo, to|x, x' 

-9^W''(xo,to|x,Xo,t)>V"(xo,to|x,x',t)-9;>V^(xo,io|x',Xo,i)>V"(xo,io|^ 
a,>V"(xo,io|x,x',i) = a;>V"(xo,to|x,x',t) = . 



By inspection, to merely serves as a parameter, and therefore we will not 
denote it explicitly in the following discussion. Also, in order to make the 
equations easier to understand, let us introduce some shorthand notation for 
the variables (xfc,x^,ifc): 





to 



v(xo,to|x,t) 



u[x + p(xo,to|^),^] , 



Xfe = (xfc,x'^,tfc) , Xk = {xk,x'k,tk) , Xfc = (xfc,x'fc) . 
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Using (38), we can now derive the dynamical equations for the statistical 
moments of W,!!: Let (■) denote a suitable ensemble averaging. We define two 
types of statistical moments: 

^«i-""(xo|Xi,...,X.„) = (>V"Hxo|Xi)...W""(xo|X„)) , 
7i"^-""(xo|Xi,...,X„) = (n(xo|Xi)>V"^(xo|X2)...>V""(xo|X„)) . 

Equation (38) implies: 

at,.F"--""(xo|Xi, . . . , X„) = - (d^^^^ + )) H"--""(xo|Xi, . . . , X„) 
_Q(xO;rMai...a„ (^^^|x^ . . . , X„) - (xo|X', Xi, . . . , X„) 

+u {dl^^ + dl,)) (xolXi, . . . , X„) , (39) 

with X = (xo,x',i) ; X'= (x, Xo,i) , with the further constraint: 

C'^-^""-""(^o|Xi, . . . , x„) = , a(fi).F"i-""(xo|Xi, . . . , x„) = . 

Equations (39), are linear and homogeneous. Therefore their solutions form a 
linear space. The most general solution to these equations is given by a linear 
combination of a suitable basis of the solutions space. To construct a specific 
solution, we must use the boundary conditions in order to set the linear weights 
of the basis solutions. We shall now show that the isotropy of these equations 
implies that our basis of solutions can be constructed such that every solution 
will have a definite behavior under rotations (that is, definite j and m). But 
before we do that, note that in many aspects the situation described here is 
similar to the well-known problem of Laplace equation in a closed domain Q: 

y2^ ^ "^lan = a . 

The Laplace equation is linear, homogeneous and isotropic. Therefore its so- 
lutions form a linear space. One possible basis for this space is: 

*im(r) = r%m{r) , 

in which the solutions have a definite behavior under rotations (belong to an 
irreducible representation of 5*0(3) ). The general solution of the problem 
is given as a linear combination of the ^j>„(r), cf. Eq. (32). For a specific 
problem, we use the value of ^(r) on the boundary (i.e., we use cr{r)) in order 
to set the values of ai^m- 

To see that the same thing happens with the hierarchy equations (39), we 
consider an arbitrary solution {jF^"'^ Ti^"^ | n = 2,3, . . .} of these equations. 
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We may write the tensor fields ^^"^ , in terms of a basis Bg 

•^"^"■°'"(xo|Xi, . . . , X„) = ^ Fgj)„(a;o, Xi, . . . , Xn)-Bgj^(xo, Xi, . . . X„) , 



'q,3m- 



'^"^■■■°'"(xo|Xi, . . . , X„) = ^ i/^j|„(xo, Xi, . . . , Xn)B^^jj^\ro, Xi, . . . X„) 



where here and below we use the shorthand notation, to denote the 

S0(3) basis of nth order tensors, Now all we have to show is that the 

pieces of ^("^'H^"^ with definite j,m solve the hierarchy equations by them- 
selves - independently of pieces with different j, m. The proof of the last state- 
ment is straightforward though somewhat tedious. We therefore only sketch it 
in general lines. The isotropy of the hierarchy equations implies that pieces of 
^^'^\'H.^'^^ with definite j, m, maintain their transformations properties under 
rotation after the linear and isotropic operations of the equation have been per- 
formed. For example, if ^°i -""(xo|Xi, . . . , X„) belongs to the irreducible rep- 
resentation (j, m), then so will the tensor fields: d^.^^ T"^"'"'^ ^ Q^k) Q^k) jrai...a„ ^ 
although, they may belong to different Sp spaces (i.e., have one less or one 
more indices). Therefore, if we choose the bases |Bg"-^| to be orthonormal, 
plug the expansion (40) into the hierarchy equations equations (39), and take 
the inner product with B^*^-^, we will obtain new equations for the scalar 

functions 



dtiFqjmi^ 0, -^1, ■ ■ • , Xn) — 

- E (fe) + ^S)) ^&(-o, X,,..., X„)B(7J), Bjl) (42) 

g' 

~ E (^M^^'^-^i'jm^^l^O, X,Xi,..., Xn)B'g^^^\ B^"j^^ 

q' 

- E {dl^'^^Fi^JU X\ X„ . . . , X„)b(5;), bJ5^) 

q' 

+^ E ( + dU)) Fiti^o, X,,..., X^)B^l, Bjl) , 

q' 



E (^0, X,,..., X„)Bg„, bJXJ)) = , (43) 

q' 

E (C^^^t (-0, Xi, . . . , X.)Bg^, B^XJ)) = . 

q' 

Note that in the above equations, (•) denote the inner-product in the Sp 
spaces. Also, the sums over q\j',m' from (40) was reduced to a sum over q' 
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only - due to the isotropy. We thus see exphcitly from (42,43) the decouphng 

of the equations for different j. m,. 

At this point we remind the reader that in the case of the most used statistical 
objects in the analysis of experimental and numerical data are the longitudinal 
nth order structure functions: 

5(")(r) = {{5Mr)T)- 

For these objects the basis functions are simply the spherical harmonics and 
the SO (3) decomposition reads: 

A question of major interest for all that follows are the numerical values of 
the scaling exponents which are defined by the power laws 

4.. 4- Dimensional Analysis of Anisotropic Fluctuations 

The actual calculation of scaling exponents in the anisotropic sectors is dif- 
ficult, and will be considered in the rest of this review. It is worthwhile to 
have a phenomenological guess based on dimensional analysis. Unfortunately, 
once anisotropies are considered, dimensional considerations becomes tricky. 
Historically, the first successful attempt to introduce dimensional considera- 
tions in anisotropic turbulence was the approach discussed in Sect. 3.4. There, 
the key role was played by the large-scale mean-shear. However this work is 
limited to the analysis of second order correlations, without discriminating 
among different j sectors. In fight of the SO (3) decomposition it should be 
considered as a prediction for Q ' . Another dimensional argument was pre- 
sented in [15] extending the consideration of Sec. 3.4. This argument takes 
into account also the particular angular structure entering in the interaction 
between small-scale fiuctuations and large-scale shear. By decomposing the 
velocity field, u, in a small-scale component, v, and a large-scale anisotropic 
component, t/, one finds the following equation for the time evolution of v. 

dtVa + vpdpVa + UpdiiV^ + vpdpUa = -daP + i^Av^. 

The major effect of the large-scale field is given by the instantaneous shear 
d/sUa which acts as an anisotropic forcing term on small scales. We can write 
the balance equation for two point quantities {vs{x')va{x)) in the stationary 
regime: 

{vs{x')V/3{x)dpVa{x)) ~ {d^UaVs{x')v^{x)) . 
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The shear term is a large-scale "slow" quantity and therefore, as far as scaling 
properties are concerned, can be safely estimated as: {d^UaVs(x')v^{x)) ~ 
Dafj,{vs{x')vn{x)) . The tensor is associated to the joint probabihty to have 
a given shear and a given small scale velocity configuration. The being a 
constant tensor can possess at most angular momentum up to j = 2. Similarly 
for three point quantities we may write: (vvvdv) ~ (dUvvv), which can be 
easily generalized to velocity correlation of any order. One may therefore argue, 
by using simple composition of angular momenta, (j = 2© j — 2), the following 
dimensional matching for structure functions in different anisotropic sectors: 

5S(r) ~ r\D\ . S^l-,^{r), (45) 

where 5'j|J^(r) is a shorthand notation for the projection on the j-th sector of 

the n-th order correlation function introduced in the previous section, Fg"}^(r). 
In (45) with |D| we denote the typical intensity of the shear term Da/s in the 
j = 2 sector. From Eq. (45) one can obtain higher j exponents of the higher 
order structure functions from the lower order structure functions of lower 
anisotropic sectors which appear on the RHS. For example, the dimensional 

(3) 

prediction for the third order scaling exponent in the j — 2 sector, Q can be 
obtained by the matching: S2^li(r) ~ r |D|5'o^^(r) ~ r^2 . By using the same 
argument and the known scaling of the third order correlation for j = 0, 2, the 
scaling exponents of the fourth order correlation for j = 2, 4 can be estimated. 
The following expression is readily obtained for any order : 

^jn) _ ( — +i) (^dimensional prediction). (46) 
3 

This formula coincides with the prediction (27) for n = 2 and j = 2. We will 
see below that both measurements and closure calculations exhibit exponents 
which are anomalous, i.e. different from these dimensional predictions. 



5 Exactly Solvable Models 

In this section we review the work done on anomalous scaling in the anisotropic 
sectors of exactly solvable models. The first of these models is the Kraichnan 
model of passive scalar advection in which the velocity field is rapidly varying 
in time [24,94, 105]. This model offers detailed understanding of the anoma- 
lous scaling in all the anisotropic sectors both from the Lagrangian and the 
Eulerian points of view. The scaling exponents can be calculated however only 
in perturbation theory. The second model that we consider is of passive advec- 
tion of a magnetic field [106]. In this case one can compute non-perturbatively 
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the scaling exponents of the second order correlation function in all the sec- 
tors of the symmetry group. These two models show that the spectrum of 
scaling exponents is discrete and strictly increasing as a function of j. If this 
is true for systems with pressure, like the Navier-Stokes equation, it may lead 
to problems of convergence of the integrals induced by the existence of the 
pressure terms. To this aim we review below a third exactly solvable model 
in which pressure is used explicitly to keep an advcctcd vector solcnoidal.lt 
was shown that also here the spectrum is discrete and strictly increasing, and 
it was explained how the putative divergences are avoided. The mechanism 
discovered here is most likely also operating in the Navier-Stokes case. The 
last model reviewed in this section is the the second order structure function 
in the Navier-Stokes problem, linearized for small anisotropics. Also in this 
case we find a discrete spectrum of strictly increasing scaling exponents as 
a function of j. Most of the results here presented can also be reproduced 
within the Renormalization Group approach. We do not enter here in this 
subject which would deserve a whole review by itself. The interested reader 
can find the most important results for passive scalar advection in [107, 108], 
for magnetic fields in [109, 110] and for passive vectors in [111]. 



5.1 Anomalous Scaling in the Anisotropic Sectors of the Kraichnan Model of 
Passive Scalar Advection 

Kraichnan's model of passive scalar advection in which the driving (Gaussian) 
velocity field has fast temporal decorrelation turned out to be a very important 
case model for understanding the anomalous scaling behavior in turbulent ad- 
vection, including the anisotropic sectors of turbulent scalar fields. We review 
here the derivation that shows that the solutions of the Kraichnan equation for 
the n order correlation functions foliate into sectors that are classified by the 
irreducible representations of the SO((i) symmetry group. A discrete spectrum 
of universal anomalous exponents is found, with a different exponent charac- 
terizing the scaling behavior in every sector. Generically the correlation func- 
tions and structure functions appear as sums over all these contributions, with 
non-universal amplitudes which are determined by the anisotropic boundary 
conditions. The isotropic sector is always characterized by the smallest expo- 
nent, and therefore for sufficiently small scales local isotropy is always restored. 
We start by presenting the Eulerian calculation which results in actual values 
of the scaling exponents (in perturbation theory) [25]. The Eulerian calcula- 
tion of the anomalous exponents is done in two complementary ways. In the 
first they are obtained from the analysis of the correlation functions of gradient 
fields. The theory of these functions involves the control of logarithmic diver- 
gences which translate into anomalous scaling with the ratio of the inner and 
the outer scales appearing in the final result. In the second way one computes 
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the exponents from the zero modes of the Kraichnan equation for the correla- 
tion functions of the scalar field itself. In this case the renormalization scale is 
the outer scale. The two approaches lead to the same scaling exponents for the 
same statistical objects, illuminating the relative role of the outer and inner 
scales as renormalization scales. To clarify this further, Ref. [25] presented an 
exact derivation of fusion rules which govern the small scale asymptotic of the 
correlation functions in all the sectors of the symmetry group and in all di- 
mensions. The purpose of the Eulerian calculation is twofold. On the one hand 
we are interested in the effects of anisotropy on the universal aspects of scahng 
behavior in turbulent systems. On the other hand we are interested in clari- 
fying the relationship between ultraviolet and infrared anomalies in turbulent 
systems. The two issues discussed in this subsection have an importance that 
transcends the particular example that we treat here in detail. Having below 
a theory of anomalous scaling in all the various sectors of the symmetry group 
allows us to explain clearly the relationship between the two renormalization 
scales and the anomalous exponents that are implied by their existence. Since 
we expect that Kolmogorov type theories, which assume that no renormaliza- 
tion scale appears in the theory, are generally invalidated by the appearance of 
both the outer and the inner scales as renormalization scales, the clarification 
of the relation between the two is important also for other cases of turbulent 
statistics. 

The central quantitative result of the Eulerian calculation is the expression 
for the scaling exponent ^j"^ which is associated with the scaling behavior 
of the n-order correlation function (or structure function) of the scalar field 

in the j'th sector of the symmetry group. In other words, this is the scaling 
exponent of the projection of the correlation function on the j'th irreducible 
representation of the SO((i) symmetry group, with n and j taking on even 
values only, n — 0,2, . . . and j = 0,2, . . .: 



n{n + d) {d + l)j{j + d~2) 



The result is valid for any even j < n, and to 0(e) where e is the scaling 
exponent of the eddy diffusivity in the Kraichnan model (and see below for 
details). In the isotropic sector (j = 0) we recover the result of [97]. It is note- 
worthy that for higher values of j the discrete spectrum is a strictly increasing 
function of j. This is important, since it shows that for diminishing scales the 
higher order scaling exponents become irrelevant, and for sufficiently small 
scales only the isotropic contribution survives. As the scaling exponent ap- 
pear in power laws of the type (r/A)^, with A being some typical outer scale 
and r ^ A, the larger is the exponent, the faster is the decay of the contribu- 
tion as the scale r diminishes. This is precisely how the isotropization of the 
small scales takes place, and the higher order exponents describe the rate of 
isotropization. Nevertheless for intermediate scales or for finite values of the 
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Reynolds and Pcclct numbers the lower lying scaling exponents will appear 
in measured quantities, and understanding their role and disentangling the 
various contributions cannot be avoided. 



5.1.1 Kraichnan's Model of Turbulent Advection and the Statistical Objects 

The model of passive scalar advection with rapidly decorrelating velocity field 
was introduced in [94]. In recent years [95-97,99, 105, 112, 113] it was shown 
to be a fruitful case model for understanding multi-scaling in the statistical 
description of turbulent fields. The basic dynamical equation in this model is 
for a scalar field T{r,t) advected by a random velocity field u{x,t): 

[dt - «;o + u{x, t) ■ V]T{x, t) = f{x, t) . (48) 



In this equation f{x, t) is the forcing. In Kraichnan's model the advecting field 
u{x,t) as well as the forcing field f{x,t) are taken to be Gaussian, time and 
space homogeneous, and delta-correlated in time: 



f{x, t)f{x', V) = ^x-x')6{t-t') , {u"{x, t)u^{x', t')) = W'^{x-x')S{t-t'). 

Here the symbols and (■ ■ ■) stand for independent ensemble averages 

with respect to the statistics of / and u which are given a priori We will study 
this model in the limit of large Peclet (Pe) number, Pe= [/aA/kq, where U\ is 
the typical size of the velocity fluctuations on the outer scale A of the velocity 
field. We stress that the forcing is not assumed isotropic, and actually the main 
goal of this section is to study the statistic of the scalar field under anisotropic 
forcing. 

The correlation function of the advecting velocity needs further discussion. It 
is customary to introduce W"^(r) via its fe-representation: 

W'^^'ir) = ^ / ^ P'^^P) eM-iP ■ r) , P'^^ip) = - ^] (49) 



where P"^{p) is the transversal projector, Vtd = {d — l)Q{d)/d and Q{d) is 
the volume of the sphere in d, dimensions. Equation (49) introduces the four 
important parameters that determine the statistics of the driving velocity field: 
A and A are the outer and inner scales of the driving velocity field respectively. 
The scaling exponent e characterizes the correlation functions of the velocity 
field, lying in the interval [0,2]. The factor D is related to the correlation 
function as follows: 

W"^(0) = ^^^^(A^ - A^) . (50) 
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The most important property of the driving velocity field from the point of 
view of the scaling properties of the passive scalar is the "eddy diffusivity" 
tensor [94] 

X°^(r) = 2[>V"''(0) - W^{r)] . (51) 



The scaling properties of the scalar depend sensitively on the scaling exponent 
e that characterizes the r dependence of K°'^{r) oc [A^ — X^]Saf3, for r ^ A, 
namely: 



OC r 



d- 1 + e 



A < r < A . 



(52) 



5.1.2 The Statistical Objects 

In the statistical theory we are interested in the power laws characterizing the 
r dependence of the various correlation and response functions of T{x, t) and 
its gradients. We will focus on three types of quantities: 

1) "Unfused" structure functions of T{x, t) are defined as 

4")(a^i,a;;,...a^„,0 = {[T{x^,t) - T{x[,t)] (53) 
X [T{x,, t) - T(4, t)]... [T(a;„, t) - T«, t)]) , 

and in particular the standard "fused" structure functions are 

S^^\T)^{[T{x + r,t)-T{x,t)r) . 

In writing this equation we used the fact that the stationary and space- 
homogeneous statistics of the velocity and the forcing fields lead to a sta- 
tionary and space homogeneous ensemble of the scalar T. If the statistics is 
also isotropic, then Sj^^ (r) becomes a function of r only, independent of the 
direction of r. The "isotropic scahng exponents" ^q"^ of the structure functions 

characterize their r dependence in the limit of large Pe, when r is in the 

/ N l/e 

"inertial" interval of scales. This range is A, -C r <^ A, L where r) = Al^j 
is the dissipative scale of the scalar field,. When the ensemble is not isotropic 
we define the exponents (47) by expanding s!f\r) according to: 

jm 
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2) In addition to structure functions we are also interested in the simulta- 
neous nth order correlation functions of the temperature field which is time 
independent in stationary statistics: 

r(-)({a;i}) = {T{x,,t)T{x2,t)...T{xn,t)), (54) 



where we used the shorthand notation {xi} for the whole set of arguments of 
nth order correlation function T^"\ Xi,X2. ■ ■ Xn- 

3) Finally, we are interested in correlation functions of the gradient field VT. 
There can be a number of these, and we denote 

7^"({a;,}) = (n[v"^T(a;,,i)]), 

i=l 

The tensor 7Y"i- -"^" can be contracted in various ways. For example, binary 
contractions ai = a2,«3 = o^a, etc. with Xi = X2,Xs = x^ etc. produces the 
correlation functions of dissipation field |VT|^. Of particular interest is the 
coordinate independent tensor Jf obtained by taking all 

^ai...a„ ^ 7^"i-"n({a;. ^ a;}) . (55) 



When the ensemble is not isotropic we need to take into account the angular 
dependence of x, and the scaling behavior consists of multiple contributions 
arising from anisotropic effects. 



5.1.3 The Eulerian calculation 

The correlation functions of the gradient field JT^"^ of Eq. (55) are tensors 
independent of the coordinates. Nevertheless their calculation is somewhat 

heavy, and we do not reproduce it here; we refer the reader to [25] where the 
calculation is presented in full detail. The final result of the calculation is for 
the projection of iS"*-"^ onto the j sector of the SO (3) symmetry group reads 

Hr<-[f)^~' , (56) 



where the proportionality constant is a tensor in the limit 77 <^ A. The ex- 
ponents are the same as those found below for the correlation function 
in which all the scales are within the inertial range. The appearance of both 
renormalization lengths and the identity of the exponents in inertial and gradi- 
ent objects is a consequence of the fusion rules that were explored in [25] with 
some care. The correlation functions T^") satisfy the Kraichnan's equation [94] 
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n 1 

1 



2 

1=1 ^ Lk=l 



Y: Hx,-x,)T^^^-'\{xi}i^,^,), (57) 
^ {i^k}=i 

where {ccji^i,^ is the set off all Xi with I from 1 to n, except of I — i and 
I — j. Substituting K°-^{x) from Eqs. (50, 51) one gets: 



-'^EV?+ E W"'^(a.,-a.,)Vrvf]T(«)({xJ) 

i=l {i^k}=l 
1 " 

= - E *(a^^ - a^fe)T("-2H{a^/ W) , (58) 
^ {i^fe}=i 

where k = kq + D[h.'" — X^]. Here we used that in space homogeneous case 
Yh=i = and therefore 

n „ n n 

|Ev.| =Ev?+ E vfvf = o. 

i=l «=1 {i^it}=l 



In this section we consider the zero-modes of Eq. (57). In other words we 
seek solutions Z^'^\{xi]) which in the inertial interval solve the homogeneous 
equation 

E K^^ixi - Xfe)Vrvfz(")({a;,}) = . (59) 

We allow anisotropy on the large scales. Since all the operators here are 
isotropic and the equation is linear, the solution space foliate into sectors 
j, m corresponding the the irreducible representations of the SO((i) symmetry 
group. Accordingly we write the wanted solution in the form 



^^"H{n})=E^S(W), 

3,m 



where zj^^ are functions composed of irreducible representations of SO((i) with 
definite j', m. Each of these components is now expanded in e. In other words, 
we write, in the notation of Ref . [97] , 

For e = Eq. (59) simphfies to 

n 

EV.^<nW) = 0, (60) 

i=l 
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for any value of j, m. Next we expand the operator in Eq. (59) in e and collect 
the terms of 0(e): 



E v?Gg({n}) = K£;S({n}) , (61) 

i=l 



where eV^ is the first order term in the expansion of the operator in (59): 



i^k=l 



[rik 



{d - ly 



ik 



v?v 



k ; 



(62) 



where r^k = Xi - Xk- 

In solving Eq. (60) we are led by the following considerations: we want scale 
invariant solutions, which are powers of r^fc. We want analytic solutions, and 
thus we are limited to polynomials. Finally we want solutions that involve 
all the n coordinates for the function Ej'^; solutions with fewer coordinates 
do not contribute to the structure functions (53). To see this note that the 
structure function is a linear combination of correlation functions. This linear 
combination can be represented in terms of the difference operator Si{x,x') 
defined by: 

6iix,x')r^-\{x,}) ^ T(-\{x,})\x,=x - T^''\{x,})\x,=x' • (63) 



Then, 

S^^\x,,x[...Xr.x'J^ll5i{xi,x'i)T^''\{x,}) . (64) 

I 



Accordingly, if T^^\{xk}) does not depend on Xi, then 5i{ri, i''j)T^'^\{rk}) — 
identically. Since the difference operators commute, we can have no contri- 
bution to the structure functions from parts of T*^") that depend on less than 
n coordinates. Finally we want the minimal polynomial because higher order 
ones are neghgible in the hmit rik ^ A. Accordingly, E^'^ with j < n is a 
polynomial of order n. Following the procedure outlined in Appendix (B) we 
can write the most general form of E^^, up to an arbitrary factor, as 

Eg =a;r...<"i?:,Vr" + [•••] , (65) 



where [. . .] stands for all the terms that contain less than n coordinates; these 
do not appear in the structure functions but maintain the translational invari- 
ance of our quantities. Note that in this case we carry the index n in the tensor 
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basis functions since the theory mixes basis functions of different orders. The 
appearance of the tensor -B^jm"" is justified by the fact that E'pl must be 
symmetric to permutations of any pair of coordinates on the one hand, and it 
has to belong to the jm sector on the other hand. This requires the appearance 
of the fully symmetric tensor (B.5). In light of Eqs. (61-62) we seek solution 
for ^^^({rfe}) of the form 

Gfliirk}) = E HUi^u}) \og{Ta) + H,^{{r,}) , (66) 

where i/j'^drfc}) and Hjm{{rk}) are polynomials of degree n. The latter is 
fully symmetric in the coordinates. The former is symmetric in r^, r; and 
separately in all the other {rk}k^i,i- Substituting Eq. (66) into Eq. (61) and 
collecting terms of the same type yields three equations: 

EV?i^iJ;. = V,-V,£;S, (67) 

i 

[d- 2 + rik - [Vi - Vk)\njm H 2d -2 — ~ 2 — ' ^ ^ 

E ^i^jm = E ^jm ■ 

i l^k 

Here Kj'^ are polynomials of degree n — 2 which are separately symmetric in 
the I, k coordinates and in all the other coordinates except l,k. In Ref. [97] it 
was proven that for j — these equations possess a unique solution. The proof 
follows through unchanged for any j ^ 0, and we thus proceed to finding the 
solution. 

By symmetry we can specialize the discussion to Z = 1, A; = 2. In light of 
Eq. (68) we see that Hj^ must have at least a quadratic contribution in ri2. 
This guarantees that (66) is nonsingular in the limit — > 0. The only part 
of Hj^ that will contribute to structure functions must contain ra . . . r„ at 
least once. Since Hj^ has to be a polynomial of degree n in the coordinates, 
it must be of the form 

Hg = rt^r'^irr ■ ■ ■ C C"^"^-"" + [■ ■ .]i,2 , (69) 

where [. . .]i^2 contains terms with higher powers of ri2 and therefore do not 
contain some of the other coordinates . . . r„. Obviously such terms are unim- 
portant for the structure functions. Since Hj^ has to be symmetric in ri,r2 
and . . .Vn separately, and it has to belong to an jm sector, we conclude 
that the constant tensor C must have the same symmetry and must belong 
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to the same sector. Consulting Appendix (B) the most general form of C is 



Substituting in Eq. (68) one find 

(rl I 9^177^^ I 12 12 ' 3 pai...a„ , „ai„Q;2raia2 &'1>2 _ r 1 

V" + ^J-"jm ^ 2(^-2 -°'»,J>" + 2 12^2^ -"-jm — L" ■ ■Jl,2 ■ 

Substituting Eq. (69) and demanding that coefficients of the term ^ . . . r^" 
will sum up to zero, we obtain 

2 

-2{d + 2)a-— — ^ = 0, -2(rf + 2)c=0 ;^c = . 

The coefficient b is not determined from this equation due to possible con- 
tributions from the unknown last term. We determine the coefficient b from 
Eq. (67). After substituting the forms we find 

Recalling the identity (B.6) we obtain b = ^^[1 — 4a] . Finally wc find that a 
is n, j-independent, a = ~ 2{d+2){d-i) ' 'whereas b does depend on n and j, and 
we therefore denote it as 6„j 

"■' 4(<i + 2)(<i-l) ■ 

In the next Subsect. we compute from these results the scaling exponents in 
all the sectors of the SO((i) symmetry group. 



5.1. 4 The Scaling Exponents of the Structure Functions 

We now wish to show that the solution for the zero modes of the correlation 
functions Fj^^ (i.e Z^")) result in homogeneous structure functions 5''^"^. In 
every sector j, m we compute the scaling exponents, and show that they are 
independent of m. Accordingly the scaling exponents are denoted and we 
compute them to first order in e. Using (63) and (64), the structure function 
is given by: 



S^tUtuTu ■ ■ ■ ; r-„, rn) = ■ ■ ■ K^BZ 



no i,i 




e ^ Ar . . . A^" (r„ n, r,, n) \aBZ^- + M^'^'Ki--^, , 
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where ' = rf ' — rf' , and the function / is defined as: 



r"^'{ri,ri, ri,ri) = {n - n)"^(ri - r,)"' In |n - r^l 

+ {Ti - Tir^if, - rip In |r, - n| (71) 

- in -riY'{ri - TiY' In In - n| 

- (r^ - riY\n - In In - r,| . (72) 

The scahng exponent of -S'^^j^ can be found by multiplying all its coordinates 
by A direct calculation yields: 



no i,i 



S^^]J^^r„ ^^r^■, . . .) = /."^^^(ri, n; . . .) - 2e/."ln/.^ ■ ■ ■ A^ Af'A^ 

no i,! 

no i,i 

= /."5j:2(n, n; . . .) - 2e/.Mn/.Ar ■ ■ ■ A- x + + 0{e') . 

no i,i 

Using (B.8), we find that E[aS^;^-"" + Vi'^'^^'^SSl = [^(^ - 1)« + 
and therefore, we finally obtain: 



S!f\iir,,im; . . .) = {1 - 2e[n{n - l)a + 6,,,] In/x} 4"Vi, n; • • •) + 0{e') 

= /.^.^"^4"^(n,n;...) + 0(6^), 



The result of the scaling exponent is: 

2(rf + 2)(d-l)+4(d + 2)(d-l)^"'^J+^^'^ 



from which follows (47). This is the final result of this calculation. 
It is noteworthy that this result is in full agreement with (56), even though 
the scaling exponents that appear in these result refer to different quantities. 
The way to understand this is the fusion rules that are discussed next. 



5.1.5 Fusion Rules 

The fusion rules address the asymptotic properties of the fully unfused struc- 
ture functions when two or more of the coordinates are approaching each other, 
whereas the rest of the coordinates remain separated by much larger scales. 
A full discussion of the fusion rules for the Navier-Stokes and the Kraichnan 
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model can be found in [53, 54, 99]. In this section we quote the fusion rules that 
were derived in Rcf. [25] directly from the zero modes that were computed to 
0(e), in all the sectors of the symmetry group. In other words, we arc after the 
dependence of the structure function ^''^^''(ri, ri; . . .) on its first p pairs of co- 
ordinates ri, ri; . . . ; Tp, Tp in the case where these points are very close to each 
other compared to their distance from the other n — p pairs of coordinates. 
Explicitly, we consider the case where ri, ri; . . . ; Tp, Tp <^ r^+i, Tp+i; . . . ; r„, r„. 
(We have used here the property of translational invariance to put the cen- 
ter of mass of the first 2p coordinates at the origin. The full calculation is 
presented in [25], with the final result (to 0(e)) 

p 

i=/max m' 

In this expression the quantity V'«,m' is a tensor function of all the coordinates 
that remain separated by large distances, and 

^max = max{0, p + j-n}, j < n. 

We have shown that the LHS has a homogeneity exponent . The RHS is a 
product of functions with homogeneity exponents ^j^^ and the functions i/ji^m' ■ 
Using the linear independence of the functions S^J^, we conclude that ■0i,m' 

must have homogeneity exponent ^j"^ — This is precisely the prediction 
of the fusion rules, but in each sector separately. One should stress the intu- 
itive meaning of the fusion rules. The result shows that when p coordinates 
approach each other, the homogeneity exponent corresponding to these co- 
ordinates becomes simply ^j^^ as if we were considering a p-order correlation 
function. The meaning of this result is that p field amplitudes measured at p 
close-by coordinates in the presence of n — p field amplitudes determined far 
away behave scaling-wise like p field amplitudes in the presence of anisotropic 
boundary conditions. 



5.1.6 The Lagmngian Approach to Anomalous Scaling 

An elegant approach to the correlation functions is furnished by Lagrangian 
dynamics [114-118]. In this formalism one recognizes that the actual value of 
the scalar at position x at time t is determined by the action of the forcing 
along the Lagrangian trajectory from t — — oo to t: 

to 

T{xo,to)= J dt{f{x{t),t))^ , (73) 

—oo 
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with the trajectory x{t) obeying 



x{to) = Xo , 



dtx{t) = u{x{t),t) + V2iir]{t) , 



(74) 



and is a vector of zero-mean independent Gaussian white random variables, 
{r]''{t)r]^{t')) = 5'^l^5{t-t'). With this in mind, we can rewrite -S^?"^ of Eq. (53) 
by substituting each factor of T{x.j) by its representation (73). Performing the 
averages over the random forces, we end up with 



To understand the averaging procedure recall that each of the trajectories Xi 
obeys an equation of the form (74), where u as well as {^71)^=1 are independent 
stochastic variables whose correlations are given above. Alternatively, we refer 
the reader to section II of [118], where the above analysis is carried out in 
detail. Here we follow the derivation of Ref. [114]. In considering Lagrangian 
trajectories of groups of particles, we should note that every initial configura- 
tion is characterized by a center of mass, say R, a scale s (say the radius of 
gyration of the cluster of particles) and a shape Z. In "shape" we mean here 
all the degrees of freedom other than the scale and R: as many angles as are 
needed to fully determine a shape, in addition to the Euler angles that fix the 
shape orientation with respect to a chosen frame of coordinates. Thus a group 
of 2n positions {Xi} will be sometime denoted below as {R, s, Z}. 

One component in the evolution of an initial configuration is a rescaling of all 
the distances which increase on the average like t^^^^ ; this rescaling is analogous 
to Richardson diffusion. The exponent ,^2 which determines the scale increase 
is also the characteristic exponent of the second order structure function [94]. 
This has been related to the exponent e of (52) according to ^2 = 2 — e. After 
factoring out this overall expansion we are left with a normalized 'shape'. It 
is the evolution of this shape that determines the anomalous exponents. 

Consider a final shape Zq with an overall scale Sq which is realized at t = 0. 
This shape has evolved during negative times. We fix a scale s > Sq and 
examine the shape when the configuration reaches the scale s for the last time 
before reaching the scale Sq. Since the trajectories are random the shape Z 
which is realized at this time is taken from a distribution 'y{Z; Zq, s Sq). 
As long as the advecting velocity field is scale invariant, this distribution can 
depend only on the ratio s/sq- 




(75) 



—00 




(76) 
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Next, wc use the shape-to-shapc transition probability to define an operator 
7(s/so) on the space of functions ^(^) according to 

[7(s/so)^](Zo) = / dZ'y{Z-Zo,s^so)^{Z). 

We will be interested in the eigenfunction and eigenvalues of this operator. 
This operator has two important properties. First, for an isotropic statistics 
of the velocity field the operator is isotropic. This means that this operator 
commutes with all rotation operators on the space of functions "^{Z). In other 
words, if Oa is the rotation operator that takes the function 'if{Z) to the 
new function ^{A~^Z), then Ca7 = 'JOa. This property follows from the 
obvious symmetry of the Kernel j{Z;Zq.s Sq) to rotating Z and Zq 
simultaneously. Accordingly the eigenfunctions of 7 can be classified according 
to the irreducible representations of SO (3) symmetry group. Because in this 
section we are not computing explicitly the exponents we do not need to 
present the precise form of the eigenfunctions and we will denote them for 
simplicity as Bgj„i{Z). The second important property of 7 follows from the 
5-correlation in time of the velocity field. Physically this means that the future 
trajectories of n particles are statistically independent of their trajectories in 
the past. Mathematically, it implies for the kernel that 

7(Z; Zo, So) ^ j dZi-f{Z; Zi, s si)-f{Zi; Zq, si ^ sq) , s > si > sq 
and in turn, for the operator, that 

7{s/so) =7(s/si)7(si/so) . 
Accordingly, by a successive application of 7(s/so) to an arbitrary eigenfunc- 

d-(2n) 

tion, we get that the eigenvalues of 7 have to be of the form agj — {s/sgp : 
(-)«. B,jrn{Zo) - / dZ-f{Z; Zo, s ^ So)B,j,^{Z) (77) 



Prom Schur's lemmas one can prove that the eigenvalues do not depend on 
m. On the other hand they can still be a function of q but for simplicity of 

notation we do not explicitly carry the q index in 

To proceed we want to introduce into the averaging process in (76) by av- 
eraging over Lagrangian trajectories of the 2n particles. This will allow us 
to connect the shape dynamics to the statistical objects. To this aim con- 
sider any set of Lagrangian trajectories that started at t — —00 and end up 
at time t = in a configuration characterized by a scale Sq and center of 
mass Rq = 0. a full measure of these have evolved through the scale L or 
larger. Accordingly they must have passed, during their evolution from time 
t — —00 through a configuration of scale s > sq at least once. Denote now 



53 



//2n(i, R, Z;s ^ So, Zo)dtdRdZ as the probability that this set of 2n trajec- 
tories crossed the scale s for the last time before reaching Sq, Zq, between t 
and t + dt, with a center of mass between R and R + dR and with a shape 
between Z and Z + dZ. 

In terms of this probability we can rewrite Eq.(76) (displaying, for clarity, 
Ro = and t = 0) as 



S^t"\Ro = 0, So, Zo,t = Q) = jdZ j dtj dRn2n{t, R,Z;s^ Sq, Z^) 

— oo 

xl^j dti---dtn (f){xi{ti) - X2{ti)) ■ ■ ■ (f){x2n~iitn) " X2nitn)) + perms (s; R, Z, t)) (78) 



The meaning of the conditional averaging is an averaging over all the realiza- 
tions of the velocity field and the random 77 j for which Lagrangian trajectories 
that ended up at time t — in R — 0, sq, Zq passed through R, s, Z at time 
t. 

Next, the time integrations in the above equation are split to the interval 
[— oo,f] and [t,0\ giving rise to 2" different contributions: 

t tot t 

j dti--- j dtn + J dti J dt2--- j dtn + ... 

—00 —00 t —00 —00 

Consider first the contribution with n integrals in the domain [— oo,t]. It 
follows from the delta-correlation in time of the velocity field, that we can 
write 



dtf- dtn (p{xi{ti) - X2{ti)) ■ ■ ■ (j){x2n-i{tn) " X2n{tn)) + perms (s; R, Z, t) 



I * 



dti--- dtn (f>{Xl{ti) - X2{ti)) - - - (f){x2n-l{tn) " X2n{tn)) + permS 



= S^t''\R, s, z, t) = ^^'(s, z) . 

The last equality follows from translational invariance in space-time. Accord- 
ingly the contribution with n integrals in the domain [— cxD,t] can be written 
as 



J dZS^T"'\s,Z) J dtJdR ii2n{t,R,Z;s^ sq,Zq) . 



(2n) 



(79) 



54 



We identify the shape-to-shape transition probabihty: 



7(Z; Zo, s ^ So) = I dtjdR //2n(i, R,Z-s^ so, -^o) ■ (80) 

— oo 

Finally, putting all this added wisdom back in Eq.(78) we end up with 

4'") (so, Zo) = / + y dZ7(Z; Zo, s ^ So)4'"^(s, Z) . (81) 

Here / represents all the contributions with one or more time integrals in the 
domain [t, 0]. The key point now is that only the term with n integrals in the 
domain [— oo,i] contains information about the evolution of 2n Lagrangian 
trajectories that probed the forcing scale L. Accordingly, the term denoted by 
I cannot contain information about the leading anomalous scaling exponent 
belonging to ^2^, but only of lower order exponents. The anomalous scaling 
dependence of the LHS of Eq.(81) has to cancel against the integral containing 
without the intervention of /. 

Representing now 



4 \sOi Zq) — ag .jjniSo)Bgjjn{Zo) , 
qjm 

4 Z) = X/ ^q,jmis)Bqjjn{Z) , 



(2n), 

qjm 
qjm 

= X] IqjmBqjmiZo) (82) 
qjm 

and substituting on both sides of Eq.(81) and using Eq.(77) we find, due to 
the linear independence of the eigenfunctions Bqj^n 

J2n) 

>-7 



O'qjmiso) — Iqjm + ( ) (^q,jm{s)- 



So. 

To leading order the contribution of Igj^ is neglected, leading to the conclu- 
sion that the spectrum of anomalous exponents of the correlation functions 
is determined by the eigenvalues of the shape-to-shape transition probability 
operator. Calculations show that the leading exponent in the isotropic sector 
is always smaller than the leading exponents in all other sectors. This gap 
between the leading exponent in the isotropic sector to the rest of the expo- 
nents determines the rate of decay of anisotropy upon decreasing the scale of 
observation. 

The derivation presented above has used explicitly the properties of the ad- 
vecting field, in particular the 5-correlation in time. Accordingly it cannot be 
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immediately generalized to more generic situations in which there exist time 
correlations. Nevertheless we find it pleasing that at least in the present case 
we can trace the physical origin of the exponents anomaly, and connect it to 
the underlying dynamics. In more generic cases the mechanisms may be more 
complicated, but one should still keep the lesson in mind - higher order corre- 
lation functions depend on many coordinates, and these define a configuration 
in space. The scaling properties of such functions may very well depend on how 
such configurations are reached by the dynamics. Focusing on static objects 
like structure functions of one variable may be insufficient for the understand- 
ing of the physics of anomalous scaling. Important confirmation of this picture 
have been found recently also for the case of passive scalars advected by a 2d 
turbulent flow in the inverse cascade regime [119] and for the case of shell 
models for passive scalars advection [120]. 

5.1.7 Summary and Discussion 

The main lesson from this subsection is that the scaling exponents form a 
discrete and strictly increasing spectrum as a function of j. This is the first 
example where this can be shown rigorously. The meaning of this result is 
that for higher j the anisotropic contributions to the statistical objects de- 
cay faster upon decreasing scales. The rate of isotropization is determined by 
the difference between the j dependent scaling exponents, and is of course 
a power law. The result shows that to first order in e the j-dependent part 
is independent of the order of the correlation function. This means that the 
rate of isotropization of all the moments of the distribution function of field 
differences across a given scale is the same. This is a demonstration of the fact 
that, to 0(e) the distributions function itself tends toward a locally isotropic 
distribution function. We note in passing that to first order in e the j de- 
pendent part is also the same for ^^^^ , a quantity whose isotropic value is not 
anomalous. For all j > 1 also ^j^^ is anomalous, and in agreement with the 

n = 2 value of Eq. (47). Significantly, for we have a nonperturbative result 
that was derived in [99], namely 

valid for all values of e in the interval (0,2) and for all j > 2. This exact result 
agrees after expanding to 0(e) with (47) for n = 2 and j — 2. 

The second lesson from this first exactly solvable example was the correspon- 
dence between the scaling exponents of the zero modes in the inertial interval 
and the corresponding scaling exponents of the gradient fields. The latter do 
not depend on any inertial scales, and the exponent appears in the combina- 

tion {A/rjp where r) is the appropriate ultraviolet inner cutoff. We found 
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exact agreement with the exponents of the zero modes in all the sectors of the 
symmetry group and for all values of n. The deep reason behind this agree- 
ment is the linearity of the fundamental equation of the passive scalar (48). 
This translates to the fact that the viscous cutoff rj is n and j independent, 
and also does not depend on the inertial separations in the unfused correlation 
functions. This point has been discussed in detail in [99,121]. In the case of 
Navier-Stokes statistics we expect this "trivial" correspondence to fail. Never- 
theless, many attempts have been done to describe the matching between the 
inertial and dissipative scaling properties [53, 122, 123] for the isotropic sector. 
Finally we note that in the present case we have displayed the fusion rules in 
all the j sectors, using the 0{e) explicit form of the zero modes. We expect 
the fusion rules to have a nonperturbative validity for any value of e. 
An interesting modification of Kraichnan models has been recently proposed in 
[124] where the scaling properties of a passive scalar advected by a Kraichnan- 
hke shear flow are investigated. The anisotropy introduced by the shear breaks 
the foliations of the correlation functions equations. Nevertheless, the authors 
have been able to explain the existence of a scaling range in the passive spec- 
trum with anomalous slope (i.e. different from the result obtained in absence 
of shear), for scales larger than the typical shear- length in the system. This 
anomalous slope is due to the fast advection of passive particles in the mean 
shear direction. 



5.2 Passively Advected Magnetic Field 



Another exactly solvable system of some interest is the case of passively ad- 
vected magnetic held. This model was flrst proposed in [106] . It describes the 
advection of a magnetic field B{x, t) by the same Kraichnan stochastic veloc- 
ity field described in Eq. (49). The equation of motion for the magnetic field 
is 



dtB{x, t) + [u{x, t) ■ Vj B{x, t) - [B{x, t) ■ Vj u{x, t) 
= KV'^B{x,t) + f{x,t) , 



(83) 



which has to be supplemented by the solenoidality condition V ■ B{x, t) — 0. 
The source ("forcing") term f{x,t) is a solenoidal vector field that is re- 
sponsible for injecting the magnetic field into the system at large scales. The 
second-order moment of the source field here is a second-order solenoidal ten- 



sor 



{nx + r,t')f{x,t)) = 5it-t')A'^^ 



(84) 
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instead of a scalar. The tensor A°''^{y) is nsed to mimic large-scale anisotropic 
boundary conditions and is therefore taken to be anisotropic, analytic in y 
and vanishing rapidly for y ^ 1. Finally, the dissipative term KV'^B{x,t) 
dissipates the magnetic field out of the system at small scales. 

Notice that in order to keep the magnetic field solenoidal, Eq. (83) contains a 
"stretching" term B{x,t)-'V u{x,t). This term may cause a "dynamo effect" , 
which is what happens when the magnetic field amplifies itself by extracting 
kinetic energy from the velocity field [125]. Such effect can destabilize the 
system, and prevents it from reaching a stationary state. 

Just as in the Kraichnan passive scalar case, we can use the fact that both 
the velocity and source fields are white-noise Gaussian processes, and derive a 
closed set of equations for the simultaneous nth-order correlation-functions of 
the magnetic field. For example, the equation of motion for the second order 
magnetic correlation function 



C°^(r, t) = (B''{x + r, t)B^{x, t 
can be easily derived [106]: 



+ (9^9,X"^)C"^'^ + 2kV'C"^ + A"f^ = f'^lpC'"' + , (85) 

where one has to add also the solenoidal condition for the magnetic field, 
daC"'^ — and the tensor K^'^ is the two-point velocity correlation (52). The 
solution of (85) was found in [106]. It was shown there that for < e < 1 
no dynamo occurs, while for e > 1 a dynamo is developed. Consequently for 
< e < 1 the system may reach a stationary state where the correlation 
function of the magnetic field behaves like a power law in the inertial range. 
In [106] the zero modes of the second-order correlation-function was calculated 
and its anomalous scaling in the isotropic sector was found for any < e < 1. 
Notice that for this passive vector model, the absence of any conservation law 
for the magnetic energy allows for anomalous scaling already for the second 
order correlation in the isotropic sector, at difference from what happens in the 
passive scalar case discussed in section (5.1.1). This was the first case where 
a fully non-perturbative analytical solutions was presented demonstrating the 
possibility to have anomalous scaling in hydrodynamic problems. 

In ref. [26, 27] this analysis was generahzed to all the sectors of the SO (3) group 
using the S0(3) decomposition. Here we review the results presented in ref. [26] 
where a systematic non-perturbative study of the solutions of (85) was given 
in all (j, m) sectors of the SO (3) group. As usual, it is advantageous to decom- 
pose the covariance C"'^ in terms of basis functions that block-diagonalize the 
angular part of the operator T,which is invariant to all rotations. In addition. 
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T is invariant to the parity transformation r — r, and to the index permu- 
tation {PjU). Accordingly, T can be further block-diagonahzed into 
blocks with definite parity and symmetry under permutations. 

In light of these consideration we seek solutions in terms of the decomposition 
given in (34): 

C-^{r,t)= J: C^iUr,t) Bf^^ir). (86) 

As discussed in sec. (4.2.2) the nine basis functions can be grouped in four 
sub-groups depending on their symmetries under parity and index permuta- 
tion (37). It should be noted that not all subsets contribute for every value 
of j. Space homogeneity implies the obvious symmetry of the covariance: 
C°-^{r, t) — C^"'{—r, t). Therefore representations symmetric to a, (3 exchange 
must also have even parity, while antisymmetric representations must have odd 
parity. Accordingly, even j's are associated with subsets I and III, and odd j's 
are associated with subset II. Subset IV cannot contribute to this theory due 
to the solenoidal constraint. 

5.2.1 The Matrix Representation of the Operator T 

Having the angular basis functions we seek the representation of the operator 
T in this basis. In such a representation T is a differential operator with respect 
to r only. In Appendix A of [26] it is shown how T mixes basis functions 
within a given subset, but not between the subsets - as is expected in the 
last section. In finding the matrix representation of T we are aided by the 
incompressibility constraint. Consider first subset I made of the four symmetric 
and with (— )•' parity basis functions: Bgj^{r) with g = 1,5,7, 9 in a given 
j, m sector. To simplify the notation we denote the a's coefficients according 
to a(r) = Ci%{r), b{r) = cif^{r), c(r) = c[%,{r) and d{r) = cif^{r). Primes 
will denote differentiation with respect to r. 
In this basis the operator T takes on the form 
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On the RHS we have matrix products. In addition, the solenoidal condition 
implies the following two constrains on a, b, c and d (cf. the Appendix of 
[23]): 
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= a' + 2- + jh' - f- + c' - j- 

= 6' + 3- + - + (j-lK-(j-l)(j-2)- 
r r r 

Using these conditions one can bring Ti and T2 to diagonal forms, 

Ti = 2(Dr' + 1 ; T2 = - [(L>r^ + + eDrl 1 

r 

where 1 is the unit matrix. T3 can be written in the form 

T,^Dr'-''Q(j,e) + Kr-^Q(j, 0). 

The exphcit expression for the four columns of Q{j, e) can be found in [26] In 
Appendix B of ([26]) the two remaining blocks (subsets II, III after the list 
(37) ), in the matrix representation of T as a function of j have been also 

investigated. The single basis -Bsjm (subset IV) cannot appear in the theory 

(2) 

since C^j„^ = by the solenoidal condition (cf. Appendix of [23]). Lastly, there 
are no solutions belonging to the j = 1 sector. This is due to the fact that 
such solutions correspond to subset II. In this subset the j = 1 solenoidal 

condition implies the equation: JrCg^^ H ^ = 0, or Cgim ^ which is 

not an admissible solution. 



5.2.2 Calculation of the Scaling Exponents 

Before turning to the computation of the exponents one should consider the 
existence of a stationary solution for t 00. In [106] it was showed that there 
is not dynamo in the isotropic sector as long as e < 1. In [26] it has been 
demonstrated that for the same values of e, the dynamo effect is absent also 
in the anisotropic sectors. The reader is referred to [26] for details on this 
subject. In the absence of a dynamo effect, we can consider a stationary state 
of the system, maintained by the forcing term /(r, t). The covariance in such 
a case will obey the following equation: 

Deep in the inertial range we look for scale invariant solutions of the above 
equation neglecting the dissipative terms. The most general scale invariant 
solution can be expressed as a linear superposition of homogeneous (zero- 
modes) and non-homogeneous solutions of the above equation: 

C'''{r)^CZm{r) + C:^on-K{r). 

In particular, only zero-modes can carry anomalous scaling, being the scaling 
properties of the non-homogeneous solutions fixed by the dimensional match- 
ing T°'^pC"'f ~ A"'^. Therefore, the existence of a leading anomalous scahng 
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contribution to small scales magnetic fluctuations is connected to the existence 
of one, some, zero-modes with scaling exponents smaller than the dimensional 
estimate. 



The calculation of the scale-invariant solutions becomes rather immediate once 
we know the functional form of the operator T in the basis of the angular ten- 
sors Bqjm- Using the expansion (86), and the fact that T is block diagonalized 
by such an expansion, we get a set of 2nd order coupled ODE's for each block. 
To demonstrate this point, consider the four dimensional block of T, created 
by the four basis tensors Bgjrn of subset I. According to the notation of the 
last section, we denote the coefficients of these angular tensors in (86), by the 
four functions a(r), 6(r), c(r), d{r): 



al3 



+ d{r)B^^j^ + ... 



where (. . .) stand for terms with other j,m and other symmetries with the 
same j, m. Let us first consider the case where ^ > 0. According to (87), well 
within the inertial range, these functions obey: 



Ti(« = 0) 



fa"\ 



b" 



+ T2{k = 0) 



(a'\ 
b' 



+ T3(« = 0) 



c 



(88) 



Due to the scale-invariance of these equations, we look for scale-invariant 
solutions in the form: 

a(r) = ar^, b{r) = br^, d{r) = cr^, d{r) = dr^ . (89) 



Where a, b, c, d are complex constants. Substituting (89) into (88) results in a 
set of four linear homogeneous equations for the unknowns a, b,c,d : 



[e(e - l)Ti(« = 0) + ^T,{k = 0) + T3(«: = 0)] 



b 

c 



. 



The last equation admits non-trivial solutions only when 

det [e(e - 1)Ti{k = 0) + ^T^iK = 0) + TsiK = 0)] = . 
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This solvability condition allows us to express ^ as a function of j and e. Using 
MATHEMATICA one finds eight possible values of ^, out-of-which, only four 
are in agreement with the solenoidal condition: 



= 1,2,3,4, 



(90) 



K{e, j) = - 2e^ + 2e^j + 2e^f - Ae^j - 3e^ - Ae^f - Sef - Sej + 4e + 16j + 16 f + 4 
H{e, j) = -e' - 8e + 2ef + 2ej + Af + 4 j + 5 . 



Not all of these solutions are physically acceptable because not all of them 
can be matched to the zero mode solutions in the dissipative regime. To see 
why this is so, consider the zero-mode equation for e = 0: 



The main difference between the e = case and the e > case is that in the 
former the same scale- invariant equation holds both for the inertial range and 
the dissipative range. As a result, for e = 0, the zero modes scale with the 
same exponents in the two regimes. These exponents are given simply by (90) 
with e = 0, because for e = the zero modes equation with k = is the same 
as (91) up to the overall factor which does not change the exponent. For 
e = th solutions should be valid for the dissipative regime as well as for the 
inertial regime, ruling out the two solutions with negative exponents in (90), 
for they will give a non-physical divergence as r — > 0. Assuming now that the 
solutions (including the exponents) are continuous in e, (and not necessarily 
analytic!), one finds that also for finite e only the positive exponents appear 
in the inertial range (an exception to that is the j = 0, to be discussed below). 
Finally there are two branches of solutions corresponding to the (— ) and (+) 
in the square root. 



Note that for j = 0, only the branch with the + sign under the square root 
exists since the other exponent is not admissible, being negative for e — > 0, and 
therefore excluded by continuity, however becomes negative as e increases. 
For j > 2 both solutions are admissible, and the leading is that one with the 
minus sign in the square root. 

Let us also discuss the behavior of the zero modes in the dissipative regime 
for e > 0. Here the dissipation terms become dominant and we can neglect all 
other terms in T. The zero mode equation in this regime becomes 2kW^C°'^ — 
0, which is again, up to an overall factor, identical to the zero mode equation 
with K = 0, e = 0. The solutions in this region are once again scale invariant 
with scaling exponents ^]^'*|e=o = j,j — 2. As expected, the correlation function 
C°'^{r) becomes smooth in the dissipative regime. 



{2n + 2D)V^C = . 



(91) 




subset I. 
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In [26] the computation of the exponents corresponding to subsets II and III 
is also presented. The result is: 

= "2 ~ 2^ ^ 2 \/l ~ + + + 2je + 4j + 4f , subset II 

= 6+ -Je^ + 2e+l + 4j2 + 2fe + 4j + 2ej , subset III. 

2 2 2 ^ 

For j = there is no contribution from this subset, as the exponent is negative. 
After matching the zero modes to the dissipative range, one has to guarantee 
matching at the outer scale L. The condition to be fulfilled is that the sum of 
the zero-modes with the inhomogeneous solutions (whose exponents are 2-e) 
must give C(r) —> as |r| —> L. Obviously this means that the forcing must 
have a projection on any sector -Bqjm for which C^j^ is nonzero. 

5.2.3 Summary and Conclusions 

The results of this section should be examined in the light of the previous 
section on passive scalars. That passive scalar case afforded only perturbative 
calculation of anomalous exponents in all anisotropic sectors. The present ex- 
ample offers exact, non-perturbative calculations, of the whole spectrum of 
scaling exponents that determines the covariance of a vector field in the pres- 
ence of anisotropy. The main conclusions are: (i) scaling exponents of the sec- 
ond order magnetic correlation functions are anomalous; (ii) they are strictly 
increasing with the index of j of the sector, meaning that there is a tendency 
toward isotropization upon decreasing the scales of observation. The equa- 
tions for the magnetic covariance foliate into independent closed equations for 
each set of irreducible representations of the SO (3) group. Moreover, scaling 
properties of the zero-modes do not show any dependence on the q index la- 
beling projections on different irreducible representations of the SO (3) groups 
for each fixed {j, m) . The consequence of the latter property is that transver- 
sal and longitudinal correlation have the same scaling exponents within each 
anisotropic sector. 

5.3 The Linear Pressure Model 

In this subsection we discuss the scaling exponents characterizing the power- 
law behavior of the anisotropic components of correlation functions in turbu- 
lent systems with pressure, exploring the fundamental question whether also 
for such systems the scaling exponents increase as j increases, or they arc 
bounded from above. The equations of motion in systems with pressure con- 
tain nonlocal integrals over all space. One could argue that the requirement of 
convergence of these integrals bounds the exponents from above. It is shown 
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here on the basis of a solvable model (the "Linear Pressure Model"), that this 
is not necessarily the case. The model described here is of a passive vector 
advection by a rapidly varying velocity field [28]. The advected vector field 
is divergent free and the equation contains a pressure term that maintains 
this condition. The zero modes of the second order correlation function are 
found in all the sectors of the symmetry group. We show that the spectrum 
of scaling exponents can increase with j without bounds, while preserving 
finite integrals. The conclusion is that contributions from higher and higher 
anisotropic sectors can disappear faster and faster upon decreasing the scales 
also in systems with pressure. To demonstrate that, consider a typical integral 
term of the form. 



Here G{r) — — l/(47rr) is the infinite domain Green function of the Laplacian 
operator, and C(r) is some statistical object which is expected to be scale 
invariant in the inertial range. If C{r) has an infrared cross over at scale L 
(or equivalently, the integral has an infrared cutoff at scale L), then the above 
expression will not be a pure power law of r, not even inside the inertial range. 
Then how is it possible that such an expression will cancel out a local term of 
C(r), as is required by the typical equations of motion? This puzzle has led in 
the past to the introduction of the concept of "window of locality" [126, 127]. 
The window of locality is the range for the scaling exponents in which no 
divergence occurs, even if the cross over length L is taken to infinity. For these 
exponents integrals of type (92) are dominated by the range of integration 
y ^ r and are therefore termed "local" . In a "local" theory no infrared cutoff 
is called for. 

In this subsection we present solutions for the scaling exponents in the anisotropic 
sectors of a linear model of turbulence with pressure. This model reveals two 
mechanisms that allow an unbounded spectrum of scaling exponents. First, a 
careful analysis of the window of locality in the anisotropic sectors shows that 
it widens as j increases. We always have a leading scaling exponent within the 
window of locality. Secondly, there is a more subtle mechanism that comes 
to play when sub-leading exponents exist outside the window of locality. In 
these cases we show that there exist counter-terms in the exact solution (not 
the zero modes!) which maintain the locality of the integrals. The bottom line 
is that in these models the anisotropic exponents are unbounded from above 
leading to a fast decay of the anisotropic contributions in the inertial range. 
The Linear Pressure model captures some of the aspects of the pressure term 
in Navier-Stokes turbulence, while being a linear and therefore much simpler 
problem. The non linearity of the Navier-Stokes equation is replaced by an 
advecting Kraichnan field u{x,t) and an advected field v{x,t). The advect- 
ing field u{x,t) is taken, as before, the Kraichnan field (49). Both fields are 




(92) 
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assumed incompressible. The equation of motion for the vector field v°'{x,t) 
is: 

dtv" + u^d/.v'^ + (9> - Kd\'' = P , (93) 

daV'' = daU"^0. (94) 

In this equation, f{x,t) is the same as the one in Eq. (83). Analyticity of 
f{x,t) is an important requirement. It means that A°'f^[x) can be expanded 
for small \x\ as a power series in x"; as a result its leading contribution in the 
j-sector is proportional to x^^^, given by d^df^x^Yj^a^x^). To see that this is 
the leading contribution the reader can consult the general discussion of the 
construction of the irreducible representations in Ref.[23]. All other analytic 
contributions contain less derivatives and are therefore of higher order in x. 

In order to derive the statistical equations of the correlation function of v°-{x, t) , 
we need a version of (93) without the pressure term. Following the standard 
treatment of the pressure term in Navier-Stokes equation, we take the diver- 
gence of (93) and arrive at. 

The Laplace equation is now inverted using the Green function of infinite 
domain with zero-at-infinity boundary conditions: 

p{x) = - JdyG{x - y)d,dX{yK{y) , 

with G{x) = —l/iiix. With this expression for p(a;), Eq. (93) can be rewritten 
as: 

dtv''{x,t)+u^ix,t)d,v''{x,t) ^ d^x^ JdyG{x - y)d,d,u^\y)v'' {y) 
- Kd%''{x,t) = r{x,t) . 

In [28] the equation of motion for the 2-point correlation function, C°^(r) = 
v'^(x + r)v^{x)^ was found: 

dtC^ir) - T''^{r) - T^'^i-r) + jdyG{r - y)df^d,T''%y) 
+ JdyG{-r - y)d''d,T^''{y) - 2Kd^C'"^{r) 
^(v''{x + r)f{x)) + (v^{x)r{x + r)) . 

where to simplify the equations we have defined an auxiliary function T"^{r): 
T'^^ir) = (v'^ix + r)u''{xy{x)) . 
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This equation is identical to the equation for the second order correlation 
function in the usual Navier-Stokes turbulence, provided that -u^ is replaced 
with v'^ in the expression above. Indeed, the vexing problem that we face 
is being made very clear: if the triple correlation function has a power law 
dependence on r with an arbitrarily large exponent, how can the integral 
converge in the infrared? One possibility is that the scaling exponent of T"^(r) 
is sufficiently low, making the integral convergent. The other possibility is that 
the correlation function is scale invariant only in the inertial range and vanishes 
quickly after that, which is equivalent to the introduction of an infrared cutoff. 
However the integral terms in the equation probe the correlation function 
throughout the entire space. Therefore, a cross over behavior of the correlation 
function at the outer scale L, seems to contradict a pure scaling behavior of 
the correlation function in the inertial range itself. This in turn implies the 
saturation of the anisotropic scaling exponents. 

To proceed, we use the fact that the field u{x,t), as well as the forcing, are 
Gaussian white noises with correlation given by Eq. (52). This enables us to 
express T"'^(r) and the correlation of the force in terms of C"^(r) and A°'^[r). 
One can use the well known method of Gaussian integration by parts [5] which 
leads to the final equations (see also appendix of ([28]): 

dtC'''^{r) = T''f^{r) + T'^"(-r) - JdyG{r - y)d^d,T''%y) 

- Jdy G{-r - y)d''d,T^''{y) + 2Kd''C''^{r) + A^'^ir) , (95) 
T'^^i.r) = ~K^^d^d,C^P{r) + U'(r) jdyG{r - y)d,[K^^-{y)d^d,C-\y) 

- ^ jdyG{y)d^dr[K^^''{y)d,d,C^^{r - y)] . (96) 

These equations have to be supplemented with two more equations that follow 
directly from the definition of C"'^(r): 

d^C^r) = , C"'^(r) = C^"(-r) . 

Finally we note that Eqs. (95,96) can be interpreted in a transparent way, 
utilizing two projection operators which maintain the RHS of Eq. (95) diver- 
gence free in both indices. To define them, let us consider a tensor field X'^^{r) 
which vanishes sufficiently fast at infinity. Then the two projection operators 
Vl and Vr are defined by: 

PL^"^(r) = X"^(r) - ^[J.) I dy G{r - y)d^X^P{y) , 
VKX-^{r) = X''^{r) - I dyG{r - y)d^X''''{y) . 
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We observe that Vj^X^f^ and V-rX'^I^ are divergence free in the left and right 
indices respectively. Using these operators we can rewrite Eqs. (95-96) in the 
form 

dtC^^ir) = VnT'^^ir) + PRr^"(-r) + 2Kd^C'"^{r) + A°^(r) , (97) 

T'^^ir) = -^Vi^K^'d^d^C'-^ir) - ^ jdyG{y)d^dr[K^''{y)d,d,C^^{r - y)] .(98) 

The projection in Eq. (98) guarantees that T^^^r) is divergence free in its left 
index, while the projection in Eq. (97) guarantees divergence freedom in the 
right index. 

Not all the terms in these equations are of the same nature. The integrals 

due to the projection operator are easy to deal with by applying a Laplacian 
on them. For example, d'^V^T'^'^ {r) = d'^T'^^{r) - d^d^T'"'{r) . On the other 
hand, there seems to be no way to eliminate the last integral in Eq. (98), 
and therefore we shall refer to it as the "non-trivial integral" . Only when the 
velocity scaling exponents in (52) are e = and e = 2 it triviahzes: the integral 
vanishes when e = and is proportional to C"^(r) when e = 2. Unfortunately, 
in these extreme cases also the projection operator trivializes, and the effect 
of the pressure cannot be adequately assessed. We prefer to study the problem 
for a generic value e for which the incompressibility constraint and the pressure 
terms are non-trivial. 

We deal with the this problem head-on in Sect. 5. 3. 4. Due to the non-trivial 
integral, we will not be able to provide a full solution of C"^(r), but only of 
the zero modes. However before doing so we would like to study a model that 
affords an exact solution in order to understand in detail the issues at hand. 
In the next section we therefore consider a simplified model of the Linear 
Pressure model, yet posing much of the same riddle. 



5.3.1 An exactly Solvable Toy model 

We construct a toy model which is inspired by equations (95, 96) for the cor- 
relation function in the Linear Pressure model. Within this model we demon- 
strate the strategy of dealing with the non-local pressure term. Since it is a 
simplification of the statistical equation of the Linear Pressure model, the toy 
model has no obvious underlying dynamical equation. 

In the toy model, we are looking for a "correlation function" C"(r), whose 
equations of motion are: 

9tC"(r) = -Ki'''{r)d^d^C'^{r) - d'(^^ jdxG{r - x)drK^"'{x)df,d^C^{x) 
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a„C"(r) = . 



(99) 



Here A°'{x) is a one-index analog of the correlation function of the original 
forces A°'f^{x). Accordingly, we take it anisotropic, analytic in x° and rapidly 
vanishing for a; ^ 1. As in the previous model, also here analyticity requires 
that the leading contribution for small x is proportional to d"x^Yj^[x) in the 
j-sector. Accordingly it is of order x^~^. 

The toy model is simpler than the Linear Pressure model in two aspects: 
First, the "correlation function", C"(r) has one index instead of two and 
therefore can be represented by a smaller number of scalar functions. Second, 
the unpleasant non-trivial term of the Linear Pressure model is absent. This 
will allow us to solve the model exactly for every value of e. Nevertheless, the 
toy model confronts us with the same conceptual problems that exist in the 
Linear Pressure model and in NS: can a scale invariant solution in the inertial 
range with a cross over to a decaying solution at scale L, be consistent with 
the integral term? If not, is there a saturation of the anisotropic exponents? 

Eq. (99) can be rewritten in terms of a new projection operator V, which 
projects a vector X"(r) on its divergence free part: 



where 



VX^{r) = X-{r) - d"^ jdyG{r - y)d^X^{y) 

We shall solve this integro-diffcrcntial equation by first turning it into a PDE 
using the Laplacian operator, and then turning it into a set of decoupled ODE's 
using the SO (3) decomposition. As in the Linear Pressure model, the non 
locality of the projection operator can be removed by considering a differential 
version of the operator: 

d^VT'^ir) = d^T'^ir) - a"a^T^(r) . 

In stationary condition 9^(7" = 0, and therefore the differential form of the 
toy model is given by: 



= K9'9'C"(r) + 9'A"(r) , 
9«C"(r) = . 



(100) 



We have reached a linear PDE of order 4. This PDE will be solved by exploiting 
its symmetries, i.e., isotropy and parity conservation, as demonstrated in the 
next subsection. 
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Eq. (100) and the incompressibility condition of C"(r) arc both isotropic and 
parity conserving. Therefore, if we expand C"(r) in terms of spherical vectors 
with a definite behavior under rotations and under reflections, we would get 
a set of decoupled ODE's for their coefficients. 

For each sector (j, m), j > of SO (3) we have three spherical vectors: 



Here ^jm{r) = r^Yjjn{r), and see [23] for further details. The flrst two spherical 
vectors have a different parity than the third vector, hence the equations for 
their coefficients are decoupled from the equation for the third coefficient. In 
the following, we shall consider the equations for the first two coefficients only, 
as they have a richer structure and larger resemblance to the Linear Pressure 
model. Finally note that the isotropic sector, i.e., j = 0, is identically zero. To 
see why, notice that in this special sector there is only one spherical vector, 
-^K)o(^) = r~^r". Hence the isotropic part of C°(r) is given by c{r)r~^r", c{r) 
being some scalar function of r. But then the incompressibility condition (99) 
implies that c(r) ~ r~^, which has a UV divergence. We therefore conclude 
that c(r) = 0, and restrict the calculation to j > 0. 

By expanding C"(r) in terms of the spherical vectors Bijmi B2jm, we obtain 
a set of ODEs (decoupled in the (j, m) labels) for the scalar functions that 
are the coefficients of these vectors in the expansion. The equations for these 
coefficients can thus be written in terms of matrices and column vectors. To 
simplify the calculations, we find the matrix forms of the Kraichnan operator 
and of the Laplacian of the projection operator separately, and only then 
combine the two results to one. 



5.3.2 The Matrix Form of the Operators and the Solution of the Toy Model 

In this subsection we derive the matrix form of the Kraichnan operator and of 
the Laplacian of the Projection operator in each j sector. To obtain the matrix 
of the Kraichnan operator in the basis of Bij^, B2jm, we expand C"(r): 



in appendix (E) we show how to find the operator on C"(r) in a matrix form 
which results in the final equation for Ci(r) and C2(r): 




r-^-V"$,„,(r) , 



C^{r) = ci(r)S" (f) + C2(r)B" (f ) . 
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,(2) 
V^2 



1= 









In addition also the incompressibility constraint daC°'{r) 
pressed as a relation between ci(r) and C2(r): 



(101) 
0, can be ex- 

(102) 



This constraint has to be taken into account when solving Eq. (101). The 
solution of Eq. (101) is somewhat tricky due to the additional constraint (102). 
Seemingly the two unknowns ci(r),C2(r) are over determined by the three 
equations (101, 102), yet this is not the case for the two equations (101) 
are not independent. To sec that this is the case and find the solution, it is 
advantageous to work in the new basis 

di^Ci+ jC2 , d2 = -2ci + j{j - 1)C2 . 

In this basis the incompressibility constraint becomes very simple: d2 = rd[, 
allowing us to express ^2 and its derivatives in terms of di. To do that in the 
framework of the matrix notation, we define the transformation matrix U: 



U 



1 j 
\-2j{j-l)J 



U 



-1 



J(J + 1) 



so that. 



f'h] 










\c2j 



. The equations of di{r) are the same as the equa- 



Pi 



tions for Q(r), with the matrices Mj replaced by Nj = UMjU ^ , and the 
sources pi replaced by 

= U 

Notice that a divergence free forcing A"{r) will cause Pi(r), P2(^) to be related 
to each other in the same way that rfi(r), d2{r) are related to each other, i.e., 
P2 — "^{Pi)' ■ Next, we perform the following replacements: 



do 



rdf^ + d\'', 



(1) 



rdf^ + 2d\'', 



(2) 



df 



rd?+3df\ 



df 



rdf^ +Adf\ 
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Wc get an equation written entirely in terms of the function di{r) and its 
derivatives: 



Pi 
[P2) 



,(103) 



where Vi are two dimensional vectors given by: 



14 = N4 



V2 = N2 



Va = N4 




+ N2 



Vo = No 







Their explicit values can be found in [28]. The Eq. (103) are for a column 
vector, and can be regarded as two scalar differential equations that we refer 
to as the "upper" and the "lower" . The upper ODE is of the fourth order, 
while the lower ODE is of fifth order. Non surprisingly, the lower equation is 
the first derivative of the upper equation, provided that A'^lr) is divergence 
free. Hence the two equations are dependent, and we restrict the attention to 
the upper equation. To simplify it, we divide both sides by Dr'^, replace di{r) 
by ip{r) and define the RHS to be the function S{r): 



S(r) = D-\-'p\{r) . 



(104) 



After doing so, we reach the following equation: 



03- 



^(3) ^(2) ^(1) 



+ 02 



+ ai^r- + ao^ = S{r) . 



(105) 



Its homogeneous solution is easily found once we substitute, il){r) — ipor^. The 
scaling exponents are the roots of the polynomial. 



^(0 = m -m- 2)(e - 3) + asm - l)(e - 2) + - 1) + + ao 

The polynomial roots are found to be real and non-degenerate. Two of them 
are positive while the other two are negative. They are given by: 



e.-(i) = - ^6 ± y^(j,e)±ysf7) i = 1, 2, 3, 4, 



(106) 
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0,5 1 1.5 2 0.5 1 1.5 2 

Fig. 1. Scaling exponents of the first few js as a function of e. Top panels show: set 
1 (left); set 2 (right). Bottom panel: set 3 (left); set 4 (right) 

where, 

A{j, e) = + ef + ej-2e + 5 + Aj + Af , 

and 

1 

B{j, e) = -Se^ j - 7e'^f + 16e^ + 2e^ + e^f - Sef - Sej - 32e + 16 + 64j + 64/. 

In the limit e ^ the roots become, in decreasing order: 

e,(l)=j + l, 0(2) = J -1, e.-(3) = -i, e.-(4) = -i-2. 

Fig. (1) displays the first few exponents as a function of e. We note that the 
spectrum has no sign of saturation as j increases. Before we discuss the mean- 
ing of this observation we will make sure that these solutions are physically 
relevant and participate in the full (exact) solution including boundary con- 
ditions. The general solution of Eq. (105) is traditionally given as the sum of 
a special solution of the non-homogeneous equation plus a linear combination 
of the zero modes. However when attempting to match the solution to the 
boundary conditions it is convenient to represent it as: 

i^{r)^y rr-. ttt dxx^-^i^'^ Six) , (107) 

' h few - 0(1)) • • • few - 0(4))/ ^ ^ ' 

all different roots 
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where the free parameters of the solution are the four constants rrii. Indeed 
a change in is equivalent to adding to the solution a term proportional to 
In the next subsection we find the values of rrii to match the boundary 
conditions, and discuss the properties of the solution. 



5.3.3 Boundary Conditions and Ineriial-range Behavior 

Prom Eq.(107) it is clear that the only values of rrii that guarantee that the 
solution remains finite as r — > and that it decays as r — > oo are mi — m2 — 
+00, 777.3 = ?774 = 0: 



^(r) = -j^rjT^ ^ ^ jdxx^-^^^^^ S{x) 

fdxx^-^^'-^^ S{x) 





(6(1) 


-0(2))(0(1)-0(3))(0(1) 


-0(4)) 


(e.-(2) 


-0(l))(0(2)-e,(3))(0(2) 


-0(4)) 








' (e,(3) 


-0(l))(0(3)-0(2))(0(3) 


- 0(4)) 








(0(4) - 


0(l))(0(4)-0(2))(0(4)- 


0(3)) 1 



r 

jdxx^-^i^^^ S{x) 


r 

Jdxx^-^^^^^ >S(x)(108) 


To understand the asymptotic of this solution we find from Eq. (104) that 
for X <^ L, S{x) has a leading term which goes like x-'"^"'^, whereas for for 
X ^ L, S{x) decays rapidly. It is now straightforward to prove that for r <^ L, 
the ^j(3),^j(4) terms scale like r^+^~''^ the 0(2) term scales like r^^^^^ and the 
^j(l) term scales like r^^^^) for values of e for which 0(1) < J + 3 — e and like 
r^+^~'^ otherwise. In addition it is easy to see that for r ^ L, ilj{r) exhibits 
an algebraic decay: the ,^^(1), ,^j(2) terms decay rapidly due to the decay of 
S{x) whereas the ,^j(3), ,^^(4) terms decay algebraically like r^" respectively. 
The asymptotic of the full solution are thus given by 



ip{r) ~ < 



(109) 



The obvious conclusion is that there is no saturation in the anisotropic scaling 
exponents as j increases. The lack of contradiction with the existence of an 
integral over all space has two aspects. The main one is simple and obvious. 
The Integro- differential equation (99) for has a differential version (100). 
Solving the differential version we are unaffected by any considerations of 
convergence of integrals and therefore the solution may contain exponents that 
increase with j without limit. Nevertheless the full solution (108) exhibits a 
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cross over at L: it increases in the inertial range r <^ L and decays for r ^ L. 
Thus plugging it back to the Integro-differential equation we are guaranteed 
that no divergence occurs. 

The question why the cross-over length L does not spoil the scale invariance 
in the inertial range still remains. The answer is found in differential form of 
the equation of motion, given by Eq. (100). From this equation we find that 
the integrand is a Green's function times a Laplacian of a tensor. By definition 
such an integral localizes, i.e. it is fully determined by the value of the tensor 
at the external vector r. In the language of Eq.(92) A(y) — V^-B(t/)! 

The second and less obvious aspect is that the window of locality widens 
up with j. This is due to the cancellations in the angular integration of the 
anisotropic solutions that are due to the orthogonality of the Yj^{r) and 
their generahzations B^j^{r). To demonstrate this consider again the simple 
integral (92), and assume that C{y) belongs to (j, m) sector, i.e. C{y) — 
a{y)Yjm{y) ■ For y ^ r, we may expand the Green function in r/y: 



1 1 oo 

G(r -y)^ ^ = —Y 

Here a„ are Taylor coefficients. Obviously the dangerous terms for the infrared 
convergence are those with low values of n. However all these terms will vanish 
for n < j due to the angular integration against Yj^iy)- The reason is that all 
these terms are of the form r^^y"-'^{r ■ y)"-'^ with < j. The angular part here 
has projections only Y^v^' with j' < < j. The first term to contribute comes 
when n — j, and is proportional to the amplitude integral J^dy y^ajmiy)y~^~^ ■ 
For a power law ajm{y) ~ y^ this implies locality for A < j — 2,, instead of 
A < — 2 in the isotropic sector. The lower bound of the window of locality is 
also extended, and a similar analysis for y ^ r leads to A > —j — 3. For the 
toy model this translates to the window of locality 

-j-e< <j + l-e. 

From the previous analysis we find that the leading power law of the full 
solution in the inertial range is r^^^'^\ which is inside this "extended" window 
of locality. Nevertheless, the subleading power r^J*^^-* originating from the first 
term in Eq. (108) is above this window, and its presence in the solution can 
be explained only using the first mechanism. 

We will see when we turn back to the Linear Pressure model that both these 
mechanisms operate there as well, leading again to a lack of saturation in the 
exponents. 
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5.3-4 Solving the Linear Pressure Model 



We now return to the Linear Pressure model. The methods used to solve it 
follow very close those developed for the toy model and therefore will not be 
described in the full way. Contrary to the toy model where we can have the full 
solution in the present case we can solve only for the zero modes. These are 
scale invariant solutions which solve an equation containing an integral. Their 
exponent must therefore lie within the "extended" (j dependent) window of 
locality. Finally one can argue that these zero modes are a part of the full 
solution that decays for r ^ L, and therefore solve the original equation as 
well. We start from Eqs (97) and (98). In the Appendix of ([28]) Eq. (98) was 
brought to the form: 

T-P^r) = -^Vi.K^-'d.d^C-^ir) - 1^-^^^-^^ Jdy G{y)y^-'d'C^P{r - t/),(110) 



which is true for every e ^ 1. The e = 1 case will not be treated here explicitly. 
Nevertheless, in [28] it was argued that that the results for e = 1 can be 
deduced from the e ^ 1 results by continuity. 

Looking at Eq. (110), we note that when e = 2, the integral on the RHS of 
the above equation trivializes to a local term C"^(r). In this limiting case the 
model can be fully solved utilizing the same machinery used in the previous 
section. The solution can then be used to check the zero modes computed 
below for arbitrary values of e. 

To proceed, we substitute Eq. (110) into Eq. (97), noting that the projector 
leaves the non-trivial integral in (110) invariant since it is divergence-free 
in both indices. Setting dtC'^^{r, t) — Q \n the stationary case, we arrive to 
following equation 



= 



(111) 



12eD 



(e-3)(e-5) 



jdyG{y)y'-^d^C''^{r - y) + 2Kd^C''f^{r) + A''^{r) 



As in the toy model, we apply two Laplacians to the above equation in order 
to get rid of the integrals of the projection operators, and obtain 



= -d" [VnVLK^''d^d,C^^] (r) 
12€D 



(112) 



(e-3)(e-5) 



dy G{y)y'-^d^C''^{r - y) + 2Kd''C'"^{r) + d^A'^ir) 



Here and in the sequel, the operator 9^" should be interpreted as (5^)". We now 
seek the homogeneous stationary solutions of C"^(r) in the inertial range (zero 
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modes). These satisfy the equations obtained by neglecting the dissipation, 
and setting the forcing and time derivative to zero: 



= d^K^'''d^d,C'^^(r) + d"d^drd^K^"'d^d,C^''{r) (113) 



+ 



Let us now define the RHS of the above equation as the "zero modes operator" 
0(e), and write the zero modes equation compactly as 







0{e)C'^^ (r) 



The solutions of this problem is obtained as before by expanding C"^(r) in a 
basis that diagonalizes 0{e). Full detail of this procedure are available in [28]. 
We turn now to discuss the results. In Fig. (2) we show the leading scaling 
exponents of the Linear Pressure model for j = 0,2,4,6,8,10. The results 
are shown for From Fig. (2), we see that in the isotropic sector and in the 
j — 2 sector, the leading exponent is ^f'^ — 0, corresponding to the trivial 
C"f^{r) = const solution. These zero modes will not contribute to the second 
order structure function, which is given by 

^"/^(r) = 2[C"/^(r) -C"^(0) 

and so we have to consider the zero mode with the consecutive exponent. In 

(2) 

the isotropic sector this exponent is exactly Q = 2 — e, as can be proven by 
passing to Fourier space. This special solution is a finger-print of the existence 
of a constant energy flux in this model. Returning to the main question of this 
subsection, we see that no saturation of the anisotropic exponents occurs since 
the leading exponent in every j > 2 sector is ~ j — 2. These exponents 
are within the window of locality of Eq. (Ill) which is given by — j — 3 < 

(2) 

Q < j — e. However the next-to- leading exponents (which are the leading 
ones in the structure function for j — 0,2), are already out of this window, 
and their relevance has to be discussed, in [28] it was proposed that the same 
mechanism that works in the toy model also operates here, and that all these 
higher exponents can be found in the full solution. To understand this, let us 
write a model equation for the correlation function in the spirit of Eq.(92): 

VC{r) + jdy K{r - y)C{y) = F{r) , (114) 



with K being some kernel, and V being some local differential operator. In view 
of Eq. (Ill), the differential operator T) should be regarded as the Kraichnan 
operator, and the integral term should be taken for all integral terms in the 
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Fig. 2. Leading scaling exponents for the first few js. The dashed hne indicates the 
upper bound of the window of locah^y 
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equation, including integrals due to the projection operators. These integrals 
create a window of locality that we denote by A^q^ < A < Aj^j. Any pure 
scaling solution C (r) ~ with A outside the window of locahty will diverge 
and hence will not solve the homogeneous part of Eq. (114). Nevertheless, 
we will now demonstrate how this zero mode can be a part of a full solution 
without breaking scale invar lance . For this we act with a Laplacian on both 
sides of Eq. (114), in order to get rid of the projection operators integrals. Of 
course, like in the Linear Pressure model, this will not eliminate all integral 
terms, and thus we can write the resultant equation as 

d^VC{r) + J dyK{r - y)cFC{y) = cFF{r) . (115) 

The main assumption, which was proven analytically in the simple case of the 
toy model, is that the above equation has a solution which is finite for all 
r, and decays for r ^ L. Let us now consider the zero modes of Eq. (115); 
their exponents have to be within the "shifted" window of locality Ajq^ + 2 < 
A < Ajjj + 2. Suppose now that with Aj^j < A < Aj^j + 2 is such a solution, 
which is therefore part of the full solution of Eq. (115). We now claim that 
this solution also solves the original equation Eq. (114), hence allowing the 
existence of scaling exponents outside of its window of locality. To see that, 
we first notice that since the full solution decays for r ^ L, then all integrals 
in Eq. (114) converge, and are therefore well defined. All that is left to show 
is that the equation is indeed solved by C(r). But this is a trivial consequence 
of the uniqueness of the solution for Laplace equation with zero at infinity 
boundary conditions. Indeed, if we denote the integral term in Eq. (114) by 

/(r)= j dyK{r-y)C{y) , 

then from Eq. (115) we have 

9'7(r) =a'[F(r)-PC(r)] , 

and since both I{r) and F{r) — VC{r) decay as r — >• oo, then they must be 
equal. Of course no breaking of scale invariance occurs because the equation is 
satisfied and F{r) — T>C{r) is a sum of an inhomogeneous solution and power 
laws. 

Returning to the Linear Pressure model, we have shown that not only the 
first, leading exponents in every sector are legitimate, but also the next few 
exponents. These exponents are inside the shifted window of locality of the 
"Laplaced" equation (113), which is given by — j + 1 < A < j + 4 — e. At this 
point, we may ask whether this is also the case for the other exponents, which 
are outside this shifted window of locality. In light of the above discussion, it 
is clear that all of them may be part of the full solution, for we can always 
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differentiate Eq. (Ill) sufficient number of times, tlius sfiifting tfie window of 
locality to include any of these exponents. However this procedure is unneces- 
sary once we have written the prefactor ^(A; j, e) as an infinite sum of poles in 
A. In that case the equation is defined for all values of A except for a discrete 
set of poles, enabling us to look for exponents as high as we wish. 



5.3.5 Summary and Conclusions 

The main question raised and answered in this subsection is whether the ex- 
istence of the pressure terms necessarily leads to a saturation of the scaling 
exponents associated with the anisotropic sectors. Such terms involve integrals 
over all space, and seem to rule out the existence of an unbounded spectrum. 
We have discussed a mechanism that allows an unbounded spectrum without 
spoiling the convergence of the pressure integrals. The mechanism is demon- 
strated fully in the context of the simple toy model, and it is proposed that 
it also operates in the case of the Linear Pressure model. The mechanism is 
based on two fundamental observations. The first one is that the window of 
locality widens up linearly in j due to the angular integration. The second, 
and more important, is that a scaling solution with an unbounded spectrum 
can exist as a part of a full solution, which decays at infinity. Indeed pure scal- 
ing solutions cannot solve themselves the zero modes equation if their scaling 
exponent is out of the window of locality. However the zero modes are always 
part of the full solution which decays to zero once r ^ L, and we have shown 
that if such a solution solves a differential version of the full equation, it must 
also solve the original equation. Therefore by differentiating the full equation 
sufficiently many times, we can always reach a differential equation with a 
window of locality as high as we wish. In that equation we can find zero mode 
solutions with arbitrarily high exponents (notice that in the toy model, it was 
sufficient to differentiate once to get rid of all integrals, thus obtaining an 
"infinitely wide" window of locality). But since these zero modes are part of 
a full solution that decays at infinity, then this solution is also valid for the 
original equation, hence showing that in the full solution there can be power 
laws with arbitrarily high exponents. Finally we want to comment about the 
relevance of these calculations to Navier-Stokes turbulence. If we substitute 
blindly e = 4/3 in our results, we predict the exponents 2/3, 1.25226, 2.01922, 
4.04843, 6.06860 and 8.08337 for j = 0, 2, 4, 6, 8 and 10 respectively. It would 
be tempting to propose that similar numbers may be expected for Navier- 
Stokes flows with weak anisotropy, and indeed for j — and 2 this is not 
too far from the truth. We return to this issue after analyzing the Navier- 
Stokes case in next section. The closeness of the Linear Pressure Model with 
Navier-Stokes equations has also been used in [128] to propose a closure for 
the non-linear turbulent problem. 
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5-4 A Closure Calculation of Anisotropic Exponents for Navier-Stokes Tur- 
bulence 



In this subsection we start from the Navier-Stokes equations, and write down 
an approximate equation satisfied by the second order correlation function, in 
a closure approximation (renormalized perturbation theory in 1-loop order) 
[129,130]. This equation is nonlinear. For a weakly anisotropic system we 
follow [130] in linearizing the equation, to define a linear operator over the 
space of the anisotropic components of the second order correlation function. 
The solution is then a combination of forced solutions and "zero modes" which 
are eigenf unctions of eigenvalue zero of the linear operator. 



5.4-1 Model Equations for Weak Anisotropy in the Closure Approximation 

It is customary to discuss the closure equations in k, t representation. The 
Fourier transform of the velocity field u{r,t) is defined by 

u{k,t) = j dr exp[—i{r ■ k)]u{x,t). 

The Navier-Stokes equations for an incompressible fluid then read 



^^ + uk']u"{k,t) = ^r"^^(fc) / ^^S{k + q + p)u*f'{q,t)u*\p,t). 



{2nY 

The interaction amplitude r"^T(fc) is defined by V'^^ik) = - [P°'^{k)kf^ + P'^l^{k)k^ 

with the transverse projection operator P"^ defined as = 5"^ — The 
statistical object that is the concern of this subsection is the second order 
(tensor) correlation function F{k,t), 

{2TifF'"^{k,t)S{k-q) = {u''{k,t)u*^{q,t)) . 

In stationary conditions this object is time independent. Our aim is to find its 
/c-dependence, especially in the anisotropic sectors. 

It is well known that there is no close-form theory for the second order simul- 
taneous correlation function. We therefore need to resort to standard closure 
approximations that lead to model equations. Such a closure leads to approx- 
imate equations of motion of the form 

^^^^M ^ laP^k, t) - ueP-^ik, t) , (116) 
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where 



/"^(fc) = / ^j0S{k + p + q)^'^^{k,q,p) . (117) 



In this equation ^°'l^{k,q,p) = ^[^"'^{k, q,p) + and 

= e{k,q,p)T''''\k)[T"''^' {q)F^^' (pf^ik) + 
r'rl^'^'{p)F^^'{q)A'^''^{k) + r^^'^'(fc)F^'5'(q)A^^'(p)] . (118) 

In stationary conditions and for k in the inertial range we need to solve the 
integral equation I°'^{k) — 0. 

The process leading to these equations is long; one starts with the Dyson- 
Wyld perturbation theory, and truncates (without justification) at the first 
loop order. In addition one asserts that the time dependence of the response 
function and the correlation functions are the same. Finally one assumes that 
the time correlation functions decay in time in a prescribed manner. This is 
the origin of the "triad interaction time" Q{k,q,p). If one assumes that all 
the correlation functions involved decay exponentially (i.e. like exp(— 7^|t|), 
then 

e{k,q,p)= \ . (119) 

For Gaussian decay, i.e. like exp[— (7^t)^/2], 

Qik,q,p)= \ (120) 



All these approximations are uncontrolled. Nevertheless this type of closure 
is known to give roughly correct estimates of scaling exponents and even of 
coefficients in the isotropic sector. 

Eq. (117) poses a nonlinear integral equation which is closed once 7^ is mod- 
eled. One may use the estimate 7^^ ~ kUk where Uk is the typical velocity 
amphtude on the inverse scale of k, which is evaluated as C/| ~ k^F°'"(k). 

7^ = C^k^'^^F<^<^{k) . (121) 



In isotropic turbulence Eqs. (117) and (121) have an exact solution with K41 
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scaling exponents, 

Ff{k)^P''^{k)F{k), F{k) = Ce^l^k-^^l\ ^k = C^e^'^k'"\ (122) 

Note that the scahng exponents in fe-representation, denoted as C, have a 
(i-dependent difference from their numerical value in r- representation. In 3- 
dimensions C^^^ = C^^^ + 3, and the exponent 2/3 turns to 11/3 in Eq.(122). 
For weak anisotropic turbulence Eq.(117) will pose a linear problem for the 
anisotropic components which depends on this isotropic solution. 



5.4-2 Closure with Weak Anisotropy 

In weakly anisotropic turbulence one has to consider a small anisotropic cor- 
rection f"^{k) to the fundamental isotropic background 

The first term vanishes with the solution (122). Linearizing the integral equa- 
tion with respect to the anisotropic correction leads to: 



I'^^{k)=J^^5ik+p + q)[S^f'^\k,q,p)r\k) + 2T-f'^\k,q,p)r\q)] = 0, 

S'^^^Hk av)^ S^"'ik,q,p) S^-P{k,q,p) 

S {k,q,p}- ^^^^^^^ , T {k,q,p)^ ^^^^^^^ . (123) 

We reiterate that the functional derivatives in Eq.(123) are calculated in the 
isotropic ensemble. In computing these derivatives we should account also for 
the implicit dependence of Q{k, q,p) on the correlation function through Eq. 
(121). We can rewrite Eq. (123) in a way that brings out explicitly the linear 
integral operator L, 

L\f) ^ <l)r'(<l) = ' (124) 

where the kernel of the operator is 

£-^^'{k,q)^Sik-q) J ^^S'^^~'\k,p,-k - p) + 2T-^^\k,q,-k - q) .(125) 



5.4-3 Symmetry Properties of the Linear Operator 

The first observation to make is that the linear operator is invariant un- 
der all rotations. Accordingly we can block diagonalize it by expanding the 
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anisotropic perturbation in the irreducible representation of the S0(3) sym- 
metry group. These have principal indices j with an integer j going from to 
oo. The zeroth component is the isotropic sector. Correspondingly the integral 
equation takes the form 

oo 

I-P(k) = lf{k) + £ lf{k) = . (126) 

The block diagonalization implies that each j'-block provides an independent 
set of equations (for every value of k): Ij^{k) = 0. The first term of (126) 
vanishes with the solution (122). For all higher values of j we need to solve 
the corresponding equation 

L\ f^) = . (127) 

We can block diagonalize further by exploiting additional symmetries of the 
linear operator. In all discussion we assume that the turbulent flow has zero 
helicity. Correspondingly all the correlation functions are invariant under the 
inversion of k. Consequently there are no odd j components, and we can write 

oo 
j=2,4,... 

We also note that in general u{—k) = u*{k). Accordingly, the correlation 
functions arc real. From this fact and the definition it follows that the corre- 
lation functions are symmetric to index permutation, F^^^k) = F,f°(fe) and 
ff^{k) = fj"{k). As a result the hnear operator is invariant to permuting 
the first {a, (5) and separately the second (7, 5) pairs of indices. In addition, 
the operator is symmetric to fc ^ —k together with q —q. This follows 
from the inversion symmetry and from the appearance of products of two 
interaction amplitudes (which are antisymmetric under the inversion of all 
wave- vectors by themselves) . Finally, the kernel is a homogeneous function of 
the wavevectors, meaning that in every block we can expand in terms of basis 
functions that have a definite scaling behavior, being proportional to k~'^. 

5.4-4 SO (3) Decomposition 

As a result of the symmetry properties the operator L is block diagonalized 
by tensors that have the following properties: 

• They belong to a definite sector {j,m) of the S0(3) group. 

• They have a definite scaling behavior. 

• They are either symmetric or antisymmetric under permutations of indices. 

• They are either even or odd in k. 
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We have already explicitly presented the tensors involved for the case of pas- 
sive vector advcction. Here we only quote the final results translated into k 
space. In every sector (j, m) of the rotation group with j > 1, one can find 9 in- 
dependent tensors that scale hke k~^. They are given by k~^Bjjj^{k), 
where the index j runs from 1 to 9, enumerating the different spherical ten- 
sors. The unit vector k = k/k. These nine tensors can be further subdivided 
into four subsets exactly like the real-space decomposition of Sect. 4.2.2.2: 

• Subset I of 4 symmetric tensors with (— parity. 

• Subset II of 2 symmetric tensors with {—Y^^ parity. 

• Subset III of 2 antisymmetric tensors with {—y^^ parity. 

• Subset IV of 1 antisymmetric tensor with (— parity. 

Due to the diagonalization of L by these subsets, the equation for the zero 
modes foliates, and we can compute the zero modes in each subset separately. 
In this subsection, we choose to focus on subset 1, which has the richest struc- 
ture. The four tensors in this subset are given here by 



Bf^M = k-^-'k^k('<P,Uk), 

B^jmik) = k-^[k-d^ + k^d^](l>im{k) , 

mU^)^k-^'^^ct>,m{k) , 

Bfjjk) = k-^+'d-d^(j>j^{k) , (128) 
where (pjmik) are the standard spherical harmonics. 

The last property to employ is the incomprcssibility of the target function 
f°'^{k). Examining the basis (128) we note that we can find two linear com- 
binations that are transverse to k and two linear combinations that are lon- 
gitudinal in k. We need only the former, which have the form 

Bf^^{k) = k-^P-^{k)<i>,m{k), (129) 
B^^^^{k)^k-=[k^<Pdf' - {j - l){k^d- + k'^d^)+j{j - l)5"^]0,„.(fc) . 

Using this basis we can now expand the target function as 

ff{k) = k-^^' [ci<^(fe) + C2<^(fe)] . (130) 

5.4-5 Calculation of the Scaling Exponents 
Substituting Eq.(130) into Eq.(127) we find 

Lq S \Bij^)ci + Lq ''^ |B2jm)c2 = . (131) 
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Projecting this equation on the two function of the basis (129) we obtain for 

~(2) " 

the matrix Li^i{j,Q ) = {Bijm\Lq \Bijrn) the form: 

Cf ) = /(03 Brjjk)C'^'''\k, q)q-'f B]^^^{q). (132) 

Here we have full integration with respect to qr, but only angular integration 
with respect to k. Thus the matrix depends on /c as a power, but we are not 
interested in this dependence since we demand the solvability condition 

detL,/j,cf ) = 0. (133) 



It is important to stress that in spite of the explicit m dependence of the 
basis functions, the matrix obtained in this way has no m dependence. In the 
calculation below we can therefore put, without loss of generality, m = 0. This 
is like having cylindrical symmetry with a symmetry axis in the direction of 
the unit vector n. In this case we can write the matrix Bij[k) (in the vector 
space a, j3 = X, y, z) as 

Br,f{k) = k-^Bf AyP,{k ■ n)) , (134) 
where are matrix operators, acting on wave vector k: 



Ual.13 

K,k^'^'--^^ ^135) 

'^2,j,k- dk'-dk^ 'yok^^dk'^ ^ ^ 

and Pj{x) denote j-th order Lcgcndre polynomials. The technical details of the 
calculations were presented in [130]. Here we present and discuss the results. 



5.4-6 Results and Concluding Remarks 

The determinants det[Lj ;(j, Q ] were computed as functions of the scaling 

~(2) 

exponents Q in every j-sector separately, and the scaling exponent was de- 
termined from the zero crossing. The procedure is exemplified in Fig. 3 for the 
isotropic sector j — 0. We expect for this sector (^'^ — 11/3, in accordance 
with Q = 2/3. Indeed, for both decay models, i.e the exponential decay 
(119), shown in dark line, and the Gaussian decay (120) shown in light line, 
the zero crossing occurs at the same point, which in the inset can be read as 
3.6667. For the higher j-sectors the agreement between the exponential and 
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Fig. 3. determinant and zero crossing for the sector j = 0. The scahng exponent 

(2) 

computed from the zero crossing is Q ' « 0.667. 




Fig. 4. determinant and zero crossing for the sector j = 2. The scahng exponent 
computed from the zero crossing is Cf ^ ^ 1.36 - 1.37. 

Gaussian models is not as perfect, indicating that the procedure is not exact. 
In Fig. 4 we present the determinant and zero crossings for j = 2. From the 
inset we can read the exponents ^2^' = 4.351 and 4.366 for the exponential 
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Fig. 5. determinant and zero crossing for the sector j = 4. The scahng exponent 

(2) 

computed from the zero crossing is 



and Gaussian models respectively. This is in correspondence with C,f^ — 1.351 
and 1.366 respectively. These numbers are in excellent correspondence with 
the experimental measurements reported in [13,29], cf. the next Section. The 
results for ?' = 4 are presented in Fig. 5. Here the zero crossing, as seen in 
the inset, yields very close results for Q between the exponential and Gaus- 
sian decay models, i.e. Q ~ 4.99. Note that this result is very close to the 
boundary of locahty as discussed in [130]. Nevertheless the zero crossing is 
still easily resolved by the numerics, with the prediction that C,f^ ^ 1.99. The 
simulation estimate of this number in [36] was 1.7 ± 0.1. We note that while 
the result Q ^ 1.99 is not within the error bars of the simulational estimate, 
it is very possible that the closeness of the exponent to the boundary of the 
window of locality gives rise to very slow convergence to asymptotic scaling. 
We therefore have to reserve judgment about the agreement with simulations 
until larger scaling ranges were available. 



Similar results arc obtained for j = 6, see Fig. 6. Also this case exhibits zero 
crossing close to the boundary of locality, with Q pa 6.98. Again we find close 
correspondence between the exponential and Gaussian models. In terms of 
C^^^ this means Ce^^ ~ 3.98. This number appears higher than the simulational 
result from [36], which estimated Q ~ 3.3 ± 0.3. We note however that for 
J = 6 the log-log plots measured in DNS [36] possess a short scaling range. 
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Fig. 6. determinant and zero crossing for the sector j = 6. The scaUng exponent 

(2) 

computed from the zero crossing is Q f» 3.98. 

Interestingly enough, the set of exponents 1.36, 1.99 and 3.98 for 

j =0, 2, 4 and 6 respectively are in close agreement with the numbers obtained 
for the linear pressure model, ^f^=2/3, 1.25226, 2.01922, 4.04843, for j = 
0, 2, 4 and 6 respectively. We reiterate at this point that the latter set is exact 
for the hnear pressure model, whereas the former set is obtained within the 
closure approximation. In fact, the close correspondence is not so surprising 
since he linearization of Navier-Stokes equations for small anisotropy results 
in a linear operator which is very close to the one that exists naturally in the 
Linear Pressure Model. Numerical results [15, 36, 131] obtained at moderate 
Re and with strong anisotropics show a small disagreement with the numbers 
calculated in the closure approximation. We do not expect a much more precise 
theoretical evaluation of these exponents before numerical and experimental 
data at higher Re are obtained and the intermittency problem in the isotropic 
sector is fully settled. 



6 Analysis of Experimental Data 

The major difficulty in applying the S0(3) decomposition to experimental data 
lies in the fact that one never has the whole field u{x). We thus cannot project 
the statistical objects onto chosen basis functions Bgj^ and simply integrate 
out all other contributions. Rather, we need to extract the wanted information 



88 



laboriously by fitting partially resolved data, or to measure quantities that do 
not have projections on the isotropic sector, to see right away the anisotropic 
contributions. We begin with the first approach. 



6.1 Anisotropic Contribution to the Statistics of the Atmospheric Boundary 
Layer 

The atmospheric boundary layer offers a natural laboratory of turbulence that 
is unique in offering extremely high Re number. Students of turbulence inter- 
ested in the scaling properties that arc expected to be universal in the limit 
Re oo are thus attracted to atmospheric measurements. On the other hand 
the boundary layer suffers from strong inhomogeneity (explicit dependence of 
the mean velocity on the height) which leads to strong anisotropies such that 
the vertical and the horizontal directions are quite distinguishable. In addition, 
one may expect the boundary layer to exhibit large-scale quasi 2-dimensional 
eddys whose typical decay times and statistics may differ significantly from 
the generic 3-dimensional case. The aim of this section is to review systematic 
methods of data analysis that attempt to resolve such difficulties, leading to 
a useful extraction of the universal, 3-dimensional aspects of turbulence. 

Obviously, to isolate tcnsorial components belonging to other than isotropic 
sectors one needs to collect data from more than one vector component of 
the velocity field. Having two probes is actually sufficient to read surprisingly 
rich information about anisotropic turbulence. In the experiments discussed 
in this subsection two types of geometry were employed, one consisting of two 
probes at the same height above the ground and the other with the two probes 
separated vertically. In both cases the inter-probe separation is orthogonal to 
the mean wind. 



6.1.1 Experiments, Data Sets and the Extraction of Structure Functions 

The results presented in this subsection are based on two experimental setups 
[10,13,29], which are denoted throughout as I and II respectively. In both 
setups the data were acquired over the salt fiats in Utah with a long fetch. In 
set I the data were acquired simultaneously from two single hot wire probes 
at a height of 6 m above the ground, with a horizontal separation of 55 cm, 
nominally orthogonal to the mean wind. The Taylor microscale Reynolds num- 
ber was about 10,000. Set II was acquired from an array of three cross-wires, 
arranged above each other at heights 11 cm, 27 cm and 54 cm respectively. 
The Taylor microscale Reynolds numbers in this set were 900, 1400 and 2100 
respectively. Table 1 lists a few relevant facts about the data records analyzed 
here. The various symbols have the following meanings: U — local mean ve- 
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Table 1 



Data sets I (first line) and II (second-fourth lines). 

locity, u' — root-mean-square velocity, e = energy dissipation rate obtained 
by the assumption of local isotropy and Taylor's hypothesis, rj and A are the 
Kolmogorov and Taylor length scales, respectively, the microscale Reynolds 
number Rx = u'X/u, and fs is the sampling frequency. 

For set 1 it is important to test whether the separation between the two probes 
is indeed orthogonal to the mean wind. (We do not need to worry about 
this point in set II, since the probes are above each other). To do so one 
computes the cross-correlation function {ui{t + T)u2{t)) . Here, ui and U2 refer 
to velocity fluctuations in the direction of the mean wind, for probes 1 and 
2 respectively. If the separation were precisely orthogonal to the mean wind, 
this quantity should be maximum for r = 0. Instead, for set I, it was found 
that the maximum shifted slightly to r = 0.022 s, implying that the separation 
was not precisely orthogonal to the mean wind. To correct for this effect, the 
data from the second probe were time-shifted by 0.022 s. This amounts to a 
change in the actual value of the orthogonal distance. The effective distance 
is A ^ 54 cm (instead of the 55 cm that was set physically). The coordinates 
were chosen such that the mean wind direction is along the 3-axis, the vertical 
is along the 1-axis and the third direction orthogonal to these is the 2-axis. We 
denote these directions by the three unit vectors n, m, and p respectively. The 
raw data available from set I is u^^\t) measured at the positions of the two 
probes. In set 11 each probe reads a linear combination of u^'^\t) and u^-^\t) 
from which each component is extractable. From this raw data we would like to 
compute the scale-dependent structure functions, using the Taylor hypothesis 
to surrogate space for time. This needs a careful discussion. 

6.1.2 Theoretical constructs: the Taylor Hypothesis, Inner and Outer Scales 

Decades of research on the statistical aspects of thermodynamic turbulence 
are based on the Taylor Hypothesis [101], which asserts that the fluctuating 
velocity field measured by a given probe as a function of time, u{t) is the same 
as the velocity u{r/U) where U is the mean velocity and r is the distance to 
a position "upstream" where the velocity is measured at t — 0. The natural 
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limitation on the Taylor hypothesis is provided by the typical decay time of 
fluctuations of scale r. Within a K41 scaling theory this time scale is the 



turn-over time r/JS{r) where S{r) = 5'""(r). With this estimate the Taylor 



Hypothesis is expected to be valid when ^JS{r)/U — > 0. Since S{r) when 

r — > 0, the Taylor hypothesis becomes exact in this limit. We will use this to 
calibrate the units when we employ two different probes and read a distance 
from a combination of space and time intervals. 

The Taylor hypothesis has also been employed when the mean velocity van- 
ishes, and instead of U one uses the root-mean-square u'. Ref.[132] has pre- 
sented a detailed analysis of the consequences of the Taylor hypothesis on the 
basis of an exactly soluble model. In particular ways were proposed there to 
minimize the systematic errors introduced by the use of the Taylor hypothesis. 
In light of that analysis we will use here an "effective" wind Ueff which for 
surrogating the time data of a single probe is made of a combination of the 
mean wind U and the root-mean-square u', 



where 6 is a dimensionless parameter. Evidently, when we employ the Taylor 
hypothesis in log-log plots of structure functions using time series measured 
in a single probe, the value of the parameter b is irrelevant, changing just the 
(arbitrary) units of length (i.e the arbitrary intercept). When we used data 
collected from two probes, we mix read distance and surrogated distance, and 
the parameter b becomes a unit fixer. The numerical value of this parameter is 
found in [132] by the requirement that the surrogated and directly measured 
structure functions coincide in the limit r ^ 0. When we do not have the 
necessary data we will use values of b suggested by the exactly soluble model 
treated in [132] . The choice of these values can be justified a posteriori by the 
quahty of the fit of to the predicted scahng functions. 

When we have two probes placed at different heights the mean velocity and u' 
as measured by each probe do not coincide. In applying the Taylor hypothesis 
one needs to decide which value of C/g// is most appropriate. This question has 
been addressed in detail in ref.[132] with the final conclusion that the choice 
depends on the velocity profile between the probe. In the case of linear shear 
the answer is the precise average between the two probes. 



where the subscripts 1,2 refer to the two probes respectively. In all the sub- 
sequent expressions we will therefore denote separations by r, and invariably 






(136) 




(137) 



91 




Fig. 8. Raw log-log plot of the transverse component of the 2nd order structure 
function. 

this will mean Taylor-surrogated time differences. The effective velocity will be 
(136) or (137) depending on having probes at the same height or at different 
heights. The value of b will be 6 = 3 following rcf.[132]. It can be shown that 
the computed scaling exponents are not sensitive to the changing b. (They 
change by a couple of percents upon changing b by 30%.) In seeking scaling 
behavior one needs to find the inner and outer scales. Below the inner scale all 
structure functions have an analytic dependence on the separation, S{r) ~ r^, 
and above the outer scale the structure functions should tend to a constant 
value. We look at the longitudinal structure functions 

computed from a single probe in set I and from the probe at 0.54m in set II , 
see Fig.7. We simultaneously consider the transverse structure function 

S^\r) = {{u\x + r) - u\x)f) 

computed from the probe at 0.54m in set II, see Fig. 8. The spatial scales 
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are computed using the local mean wind in both cases since the scaling ex- 
ponent for the single-probe structure function are not expected to be affected 
by the choice of convection velocity. This choice does determine the value of 
r corresponding to a particular time scale however. One may expect that any 
correction to the numerical value of r is small for a different choice of con- 
vection velocity, and not crucial for the qualitative statements that follow. In 
Fig. 7 we clearly see the behavior characterizing the transition from the 
dissipative to the inertial range. As is usual, this behavior persists about a 
half-decade above the "nominal" Kolmogorov length scale rj. There is a region 
of cross-over and then the isotropic scahng ~ r°-^^ expected for small scales in 
the inertial range begins. We thus have no difficulty at all in identifying the 
inner scale, it is simply revealed as a natural crossover length in this highly 
resolved data. We understand by now that we cannot expect to be able to fit 
with this single exponent for larger scales and must include scaling contribu- 
tions due to anisotropy. We expect that the contributions due to anisotropy 
will account for scaling behavior up to the outer scale of a 3-dimensional flow 
patterns. The question therefore is how to identify what this large scale is. 
One approach would be to simply use the scale where the structure function 
tends to a constant, which corresponds to the scale across which the velocity 
signal has decorrelated. It becomes immediately apparent that this is not a 
reasonable estimate of the relevant large scale. Fig. 7 shows that the structure 
function stays correlated up to scales that are at least an order of magnitude 
larger than the height at which the measurement is made. On the other hand, 
if we look at the transverse structure function computed from the probe at 
0.54m, Fig. 8 we see that it ceases to exhibit scaling behavior at a scale that 
is of the order of the height of the probe. It appears that we are observing 
extremely flat eddys that are correlated over very long distances in the hor- 
izontal direction but have a comparatively small correlation lengths in the 
direction perpendicular to the boundary. Since we know that the presence of 
the boundary must limit the size of the largest 3-dimensional structures, the 
height of the probe should be something of an upper bound on the largest 3- 
dimensional flow patterns that we can detect in experiments. Thus we arrive 
at a qualitative understanding of the kind of flow that is observed in these 
atmospheric measurements. The size of the largest 3-dimensional structures 
is determined by the decorrelation length of the transverse structure func- 
tion. This is because the transverse components of the velocity are unaffected 
by the extended, persistent, 2-dimensional eddys that govern the behavior of 
the longitudinal components. The theory of scaling behavior in 3-dimensional 
turbulence can usefully be applied to only those flow patterns that are truly 3- 
dimensional. The extended flat eddys must be described in terms of a separate 
theory, including maybe notions of 2-dimensional turbulence which has very 
different scaling properties [133]. Such considerations are outside the scope of 
this review. Rather, in the following analysis, the outer-scale was chosen to 
be of the order of the decorrelation length of the transverse structure function 
(where available) or of the height of the probe. We will see below that up to 
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a factor of 2 these are the same; taking L to be as twice the height of the 
probe is consistent with all data. We use this estimate in the study of both 
transverse and longitudinal objects. 

6.1.3 Eoctracting the universal exponents of higher j sectors 
We consider the second order structure function 

5"^(r) = ((k"(x + r) - k"(x))(k^(x + r) -u^{^))). (138) 

The lowest order anisotropic contribution to the symmetric (in indices), even 
parity (in r due to homogeneity), second-order structure function is the j = 2 
component of the SO (3) symmetry group. Ref. [29] presents a derivation of the 
m — axisymmetric (invariant under rotation about the 3-axis) part of the 
j — 2 contribution to this structure function in homogeneous turbulence. The 
derivation of the full j = 2 contribution to the symmetric, even parity structure 
function appears in Appendix C. Fig. 9 b shows the fit to the structure function 
computed from a single probe in set I 

S^'(r, e = 0) = {(u't\x + r) - uf\x))^), (139) 

where the subscript 1 denotes one of the two probes, with just the j — 
contribution. The best-fit exponent for the range < r/A < 4.5 is Q — 
0.68 ± 0.01 (Fig. 9 a). Above this range, was impossible to obtain a good fit 
to the data with just the isotropic exponent and Fig. 9 b shows the peel-off 
from isotropic behavior above r/A = 4. 



To find the j — 2 anisotropic exponent one needs to use data taken from the 
two probes. To clarify the procedure we show in Fig. 10 the geometry of set I. 
What was computed is actually 

S'%r, 0) = {[u^\u,fft + U^^ft,) - u^^\u,fft)f) 

Here 6 — arctan(A/C/e//if), tf — f/Ueff^ and r = /S? -\- {pefftfY. U^ff was 
as in Eq.(136) with 6 = 3. We will refer from now on to such quantities as 

S^'^r.e) = {{uf{x + r)-vS\x)f) . (140) 

Next, one may fix the scaling exponent of the isotropic sector as 0.68 and 
find the j = 2 anisotropic exponent that results from fitting to the full j = 2 
tensor contribution. Finally, one needs to fit the objects in Eqs. (139) and 
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Fig. 9. The single-probe structure function computed from data set I. (a) shows 
the minimization by the best-fit value of the exponent in the isotropic sector 

/OA 

Q ~ 0.68 for the single-probe structure function in the range < r/A < 4.5. (b) 

(2) 

shows the fit using the best value of Q obtained in (a) , indicating the peel-off from 
isotropic behavior at the end of the fitted range. 




Rv/hO = A 



Fig. 10. Diagrammatic illustration of the experimental set-up. Shown is the posi- 
tioning of the probes with respect to the mean wind and how the Taylor hypothesis 
is employed 

(140) to the sum of the j = (with scaling exponent Q = 0.68) and the 
j = 2 contributions (see Appendix C) 



S^\t^ 9) = 5|io(r, e) + Sfi,(r, 9) = c, ( ^ ) ° [2 + cf - ci'^ cos^ 9 



+ al- 



A2) 

r \ 



(^(2) + 2f - Cr'(3Cr^ + 2) cos^ 9 + 2Cr^(Cr^ - 2) cos^ 9ig] 



.(2).oa{2) 



{2)/^{2) 



+ 6 



a(2) 



(d'^ + + 3) - cr^(3cr^ + 4) cos'^ 9 + (2cr^ + - 2) cos^ 9 



.(2) 



'•(2)/qa(2) 



(2) 



.(2) 
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Fig. 11. The ^ minimization by the best-fit value of the exponent in the j = 2 
anisotropic sector from the fit to both the 9 = and the ^-dependent structure 
function in the range < r/A < 25. 
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A. 
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A. 





- 2Cf ^(d^^ + 2) sin ^ cos ^ + 2CP{ci^^ - 2) cos^ 6 sin 6 
-2Cr(Cr-2)cos2 0sin2^; 
-2Cf(Cr-2)sin^4 

(2) 



The above fit was performed using values of Q ranging from 0.5 to 3. The 
best value of this exponent is the one that minimizes the for the fits. From 
Fig. 11 one may read the best value to to be 1.38 ± 0.15. The fits with this 
choice of exponent are displayed in Fig. 12. The corresponding values of the 
5 fitted coefficients can be found in the paper [29]. The range of scales that 
are fitted to this expression is 1 < r/A < 25. We thus conclude that that 
structure functions which is symmetric in r exhibits scaling behavior over 
the whole scaling range, but this important fact is missed if one does not 
consider a superposition of the j = and j = 2 contributions. Finally, let 
us note that the value of the exponent is perfectly in agreement with the 
analysis of numerical simulations [34] in which one can comfortably integrate 
the structure function against the basis functions, eliminating all contributions 
except j = 2 (see next section). 



(141) 



6.1.4 Extracting the j=l component 

In homogeneous fiows, it follows from symmetry and parity that the often com- 
puted and widely analyzed structure function as defined in Eq. (138) does not 
possess any contribution from the J = 1 sector. The lowest order anisotropic 
contribution belongs to the j = 2 sector. In order to isolate the scaling be- 
havior of the J = 1 contribution in atmospheric shear fiows we must either 
explicitly construct a new tensor object which will allow for such a contribu- 
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5 10 15 20 25 30 5 10 15 20 25 30 



r/A r/A 

Fig. 12. The structure functions computed from data set I and fit with the j = and 

(2) 

full j = 2 tensor contributions using the best fit values of exponents Q = 0.68 and 
Q ' = 1.38 in the range < r/A < 25. Panel (a) shows the fit to the single-probe 
(6 = 0) structure function and panel (b) shows the fit to the 6'-dependent structure 
function. 



tion, or see if it can be extracted from the structure function itself computed 
in the case of inhomogeneity. We have pursued both avenues. In the former, 
we construct the tensor 



T"/5(r) = {u'^{x + r) - u'^{x)){u^{x + r) + u^ix))). (142) 



It is easily seen that the function vanishes both in the case of a = (3 and 
when r is in the direction of homogeneity. From data set II we can calculate 
this function for non- homogeneous (in the shear direction) scale-separations. 
In general, this will exhibit mixed parity and symmetry and therefore, to 
minimize as far as possible the final number of fitting parameters we look at 
only the antisymmetric contribution. We derive the tensor contributions in 
the J = 1 sector for the antisymmetric case in Appendix D and use this to 
fit for the unknown j = 1 exponent. Below we describe the results of this 
analysis. Next, we computed the 6'-dependent structure function from set II. 
We expect that this could exhibit the j = 1 component as inhomogeneity 
does not allow us to apply incompressibility in the different symmetry and 
parity sectors to eliminate this contribution as in the case of the homogeneous 
structure function. This structure function is symmetric but of mixed parity. 
We derive the tensor contributions in the j = 1 sector for the symmetric case 
in Appendix D and use this to fit for the j = 1 exponent. 
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Fig. 13. Right: the minimization by the best-fit value of the exponent C)^'' of the 
J = 1 anisotropic sector from the fit to ^-dependent T^^(r, 0) function in the range 
< r/A < 2.2. Left: The fitted f'^^{r,Q) function. The dots indicate the data and 
the hne is the fit. 

6.1.5 Antisymmetric Contribution 

We consider the tensor object in Eq. (142). In order to have as few parameters 
as possible in the fitting procedure, we take the antisymmetric part 



which will only have contributions from the antisymmetric j = 1 basis ten- 
sors. An additional useful property of this object is that it does not have any 
contribution from the isotropic helicity-free j = sector due to its antisymme- 
try. This allows us to isolate the j = 1 contribution and determine its scaling 



(2) 

exponent Ci starting from the smallest scales available. Using data from the 
probes at 0.27m (probe 1) and at 0.11m (probe 2) we calculate 



where again super-scripts denote the velocity component and sub-scripts de- 
note the probe at which this component is measured. We want to fit this object 
to the tensor form derived in Appendix (D), namely: 



T^\r, ^,0 = 0) = -a3,i,or^i sin 6 + a2,i,ir^i'' + a3,i,_ir^i^' cos 9. 

Fig. 13 gives the minimization of the fit as a function of ([^^ and we use 
the best value of 1 ± 0.15 for the final fit. This is shown in the left panel. The 
fit in Fig. 13 peels off at the end of the fitted range. The maximum range over 
which one can fit is of the order of the height of the probes and again, this is 
consistent with the considerations presented above. 



(2) 



.(2) 
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Fig. 14. Right: the minimization by the best-fit value of the exponent (^[^^ of 
the j = 1 anisotropic sector from the fit to ^-dependent inhomogeneous structure 
function in the range < r/A < 2.2. Left: The fit to the inhomogeneous structure 
function computed as in Eq.(143). The dots indicate the data and the hne is the fit. 



6.1.6 Symmetric Contribution 



Finally, we compute the structure function Eq. (140) where the subscripts 
denote probe 1 at 0.27m and probe 2 at 0.11cm. As discussed in Appendix 
(D), since the scale separation has an inhomogeneous component, we expect 
a contribution from the j = 1 anisotropic sector and we would like to extract 
what the scaling exponent in this sector is. Note that the j = sector con- 
tributes two independent tensor forms with coefficients we will denote by Ci 
and C2, since incompressibility does not provide a constraint to relate them. 
This fact combined with Eq. (D.8) gives us the the tensor form to which we 
must fit our function 



S^^{r,9) = cir'^o + C2/0 cos^ 6* + ai,i_or'^i cos6' + a7,i,or^i 2cos6' (143) 

a(2) o a{2) ,(2) ^(2) „ 

+'j^9,i,o'" cos^ + ag^i ir^i (— 2 cos^^sin^) + ai i__ir^i sin 6' + ag^i^^ir''! cos 6'sin^ 
(2) 

We fix the exponent Q to be 0.68 and perform fits with varying values of 
(i^'' for 8 unknown coefficients. The best value of ([^^ is obtained for the range 
< r/A < 4.2 and is 1.05 ±0.15 as is shown in Fig. 14. In the left panel we 
show the fit to the data using this value of the exponent. -8.2 The fit peels off 
at the end of the fitted range at the scale on the order of twice the height of 
the probe, consistent with the earlier discussion. There does not exist a well- 
defined ([^^ as given by the standard minimization procedure for ranges 
smaller or larger than that fitted for in Fig. 14. The quality of the fit is good 
although, as was expected from the large number of parameters in the fitting 
function Eq. (143), ^is a function of the ([^^ is not as smooth as for all 
previous fits and its minimum is a relatively weak one. Therefore, we present 
this result mainly as it provides support to that of the antisymmetric case in 
the previous section, dots 
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6.1.7 Summary and conclusions 

In summary, we considered the 2nd order tensor structure functions of veloc- 
ity differences in the atmospheric boundary layers. The following conclusions 
appear important: 

(1) The atmospheric boundary layer exhibits 3-dimensional statistical tur- 
bulence intermingles with activities whose statistics are quite different. 
The latter are eddys with quasi-two dimensional nature, correlated for 
hundreds of meters, having little to do with the three-dimensional fluc- 
tuations discussed above. We found that the 

(2) We found that the "outer scale of turbulence" as measured by the three- 
dimensional statistics is of the order of twice the height of the probe. 

(3) The inner scale is the the usual dissipative crossover, which is clearly seen 
as the scale connecting two different slopes in log-log plots. 

(4) Between the inner and the outer scales the sum of the components up 
to j = 2 appears to offer an excellent representation of the structure 
function. 

(5) The scaling exponents Cj^^ are measured as 0.68±0. 01, 1±0.15, 1.38±0.10 
for j' = 0, 1, 2 respectively. 

We note that as far as the low order j sectors are concerned, the picture that 
emerges for Navier-Stokes turbulence is not different from the linear advec- 
tion problems that were treated in the previous section. If the trends seen here 
continue for higher j values, we can rationalize the apparent tendency toward 
isotropy with decreasing scales. If indeed every anisotropic contribution intro- 

duced by the large scale forcing (or boundary conditions) decays as {r/Lyi 
with increasing Q as a function of j, then obviously when r/L ^ only 
the isotropic contribution survives. This is a pleasing notion that justifies the 
modeling of turbulence as isotropic at the small scales. 

6.2 Homogeneous Shear 

In this subsection we discuss recent experiments in which anisotropy is created 
without inhomogeneity [7,30,31]; such experiments are particularly appealing 
for our purposes. Homogeneous-shear flow can be realized in a wind tunnel by 
using a variable solidity screen followed by flow straighteners. Such a set-up 
results in a shear flow that remains approximately constant for the length 
of the tunnel. To produce high Reynolds numbers one places an active grid 
before the shear generating screen [7]. In this way Re\ can be as high as 1000. 
To assess directly the effects of anisotropy it is useful to measure statistical 
objects that vanish identically in the isotropic sector. A possible choice is the 
set of skewness and hyper-skewness [14,38,134] as explained in Sect. (2.5). 
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Other purely anisotropic inertial range observables can be defined by mixing 
longitudinal and transversal increments with an odd number of transversal 
components: 



Systematic measurements of these anisotropic mixed correlation functions was 
reported in [10,30,31]. From the experimental data it is not possible to ex- 
actly disentangle different anisotropic projections in different sectors. This is 
because SO (3) projection requires the knowledge of the whole velocity field in a 
3d sub-volume, something clearly out of reach in any experimental apparatus. 
The simplest working hypotheses one can make is that, due to the hierar- 
chical organization of anisotropic scaling exponents, the statistical behavior 
of quantities as (144) is dominated, at scales small enough, by the leading 
j — 2 sector. In other words, the experimental measurements of the scaling 
properties of (144) is the best estimate of the exponent ^i"-P+^i+'^\ [30, 31] 
the plots of purely anisotropic quantities like (144) up to order n — S with 
n — p+2q+l were shown. The data clearly shows that these purely anisotropic 
structure functions have quite good power-laws behavior with exponents that 
are sub-leading with respect to the exponents of the isotropic structure func- 
tions of the same order, n. For example S^^'^\r) ~ r^'^^ while the fourth-order 
longitudinal structure function in isotropic ensembles is known to scale as 
S^^'^\r) ~ 7-1-27 Similar qualitative and quantitative results were obtained by 
analyzing data from an atmospheric boundary layer in [10] and in the bound- 
ary layer close to a wall [33]. In the latter two works, a phcnomcnological fitting 
procedure to the large scale behavior allowed the authors to find a power law 
for the anisotropic structure functions which pertain to a much larger range 
of scales. We draw the reader's attention to the discrepancy in the best fit for 
the scaling exponents founds for S^^'^\r) and S'^^'^\r) in [30,31]. Similar dis- 
crepancies are also reported for higher order structure functions. In our view, 
this cannot be taken as evidence that there is a g-dependence of the scaling 
exponents of the S0(3) projections. First, the anisotropic exponents are rel- 
atively inaccurate due to statistical errors; the amplitudes of the anisotropic 
fiuctuations are relatively small. Second, as already said, the experimental 
data cannot disentangle exactly the contribution of the j = 2 sector. There- 
fore, it may well be that contributions from the j = 4 (and higher) sectors 
affect differently the correlation functions with different tensorial structure. 
Similarly, other experimental investigation focused on the SO (3) decomposi- 
tion [32, 135, 136] have found results depending on the geometric set-up of the 
analyzing probes. The experimental analysis of anisotropic turbulence via the 
SO (3) decomposition is at its infancy; more refined experimental techniques 
are needed before a firm conclusion can be reached on these issues. 




(144) 
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6.2.1 Explanation of Persistence of Anisotropics 

As discussed in Sect. (2.5) tliere are numerical and experimental evidences 
of the persistence of small-scale anisotropic fluctuations in various instances 
[7,10,21,35,55]. The issue has many important consequences. We would like 
to refer to the violation of the return-to-isotropy in different meanings [35]. A 
strong violation would be implied if the following set of inequalities between 
different anisotropic exponents of the same correlation function were broken: 

Ct^ < Cf ) < • • • < Cf \ (145) 



i.e. if one, or more, anisotropic sector becomes leading with respect to the 
isotropic one. This would destroy the phenomenology of turbulence as de- 
veloped since Kolmogorov's theory in 1941. Turbulence would become more 
and more anisotropic at smaller and smaller scales. As a result, strong non- 
universalities in small-scales statistics would show up depending on which 
anisotropic sector is switched on/off by the large-scale forcing. A strong vio- 
lation of the return-to-isotropy postulate has never been observed in Navier- 
Stokes turbulence. On the other hand, when the hierarchy (145) holds, any 
dimensionless anisotropic observables made of ratios between anisotropic and 
isotropic projections of the same correlation function vanishes in the small- 
scales limit. For example, focusing on the decomposition of longitudinal struc- 
ture functions (44) we may write: 

lim4^ ^ /i"^-^o"' ^ 0. (146) 
-o(5S)(r)) 



A new phenomenon occurs when anisotropic fluctuations are assessed by using 
dimensionless observables made of c?ij9^eren^ correlation functions. For instance, 
by using again the S0(3) decomposition of longitudinal structure function (44) 
one may build up anisotropic observables defined as: 

RirriVf) = — 74{ ~ ^ With Xi — Ci TrCn ■ (147) 



This is the nth order moment of the velocity probability density function, 
normalized by its isotropic second order moments. The quantities defined in 
(147) must be exactly zero in isotopic ensembles, and should go to zero as 
power laws, i?^^(r) ~ r-'/^, in an anisotropic ensemble in which the dimen- 
sional scaling (46) is satisfied. On the other hand, results from experiments 
and numerics show a much slower decay, and, in some cases, no decay at 
all [7,35]. We refer to this phenomenon as weak violation of the return-to- 
isotropy. Such a weak violation is not in contradiction with the inequalities 
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(145); there the relative importance of anisotropic fluctuations with respect 
to isotropic fluctuation of the same correlation function are implied. The vi- 
olation of the dimensional recovery-of-isotropy is simply due to the existence 
of anomalous scahng in the anisotropic sectors. Indeed, in this case, the expo- 
nents, xf\ governing the LHS of (147) can assume values much smaller than 
the dimensional estimate (including negative values!). This is exactly what is 
observed in the experiments and numerics. From Table 3 one realizes that due 
to the presence of anomalous scaling in the anisotropic sectors we have a slow 
recovery-of-isotropy, in agreement with what was explained before. 

The anisotropic observables built in terms of the generalized flatness or skew- 
ness discussed in section (2.5) are nothing but Eq. (147) evaluated at the 
dissipative length scale, r = rj. Therefore, the "persistence-of-anisotropies" 
discussed in [7, 21] can be explained invoking the very same reasoning. 



6.2.2 Summary of experimental results: universality of the anisotropic sectors 

Comparing the results obtained in [10,30,31] the following picture emerges. 
First, all the correlation functions up to n = 10, show anomalous scaling be- 
havior, where anomalous is meant with respect to the dimensional Lumley-like 
prediction discussed in Sect. (4.4). Second, the values of scaling exponents ex- 
tracted from the two different experiments [31] and [10] are in good qualitative 
agreement (see Table 2). This is an important first confirmation of the uni- 



(p,2g+l) 


(1,1) 


(1,3) 


(3,1) 


(5,1) 


(3,3) 


(1,5) 


(7,1) 


(5,3) 


(3,5) 


WS 


1.05 


1.56 


1.42 


2.02 


1.89 


1.71 


2.33 


2.22 


1.99 


KS 


1.22,1.12 


1.58,1.54 






2.14,2.00 











Table 2 



Measured scaling exponents for S^P''^'^^^\r) of various orders in two experiments. 
WS corresponds to [31] and KS to [10] 

versality of scaling exponents in the j = 2 anisotropic sector. Finally, there 
exists a clear hierarchy between isotropic and anisotropic scaling exponents, 
the latter being always larger for any given order, n of the correlation function. 
This hierarchical organization is the necessary and sufficient requirement for 
the return-to-isotropy to hold, i.e. the small scales statistics of any correlation 
function is dominated by the isotropic fluctuations. Nevertheless the gap be- 
tween isotropic and anisotropic exponents, Co""* ~ C'2'\ tends to shrink when 
n increases, implying that anisotropic contributions may exhibit important 
sub-leading effects also at very high Re. 
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7 Analysis of DNS Data 



Direct numerical simulations of turbulence are natural grounds where the util- 
ity of the S0(3) decomposition can be exploited to its maximum benefit. The 
reason for this is that numerical simulations, in contrast to current exper- 
iments, provide access to the full velocity field at all points of the turbu- 
lent domain. Therefore, the full SO (3) decomposition can be realized, without 
the constraints of best-fits to partial data. Given a tensor structure function 
S^^\r), cf. Eq.(31), we can integrate it against the spherical tensors, S^^^(f) 
[e.g., (35)], on a sphere of radius r. These integrations yield the projection of 
the structure function on the different sectors of the SO (3) group, by virtue 
of the orthogonality of the basis tensors. On the other hand DNS suffer from 
limited Reynolds numbers; consequently they have relatively short inertial 
ranges. 

Prior to the introduction of the SO (3) decomposition the numerical investiga- 
tions of anisotropic flows were focused on either single point or two-points cor- 
relations, limited, often, to the analysis of the Fourier transforms in wavevector 
space. The most recent, highly-resolved, numerical investigation of this kind 
was reported in [37]; there the full tensorial properties of the Fourier transform 

dcf ■ 

of the two-point velocity correlation, Q"^{k) = J dre^ {u'^lx + r)u^{x)), 
were calculated in a homogeneous shear [6, 137]. The main result is a confirma- 
tion of Lumley's prediction for the scaling exponent of the purely anisotropic 
co-spectrum: 



where a ^ 13 to eliminate the isotropic contribution. Only recently DNS were 
performed to probe the anisotropic component in a systematic way by exploit- 
ing the S0(3) decomposition [15,34,36,131,138]. Here we review the main 
findings, showing that 

(1) The scahng laws (log-log plots) at moderate Reynolds numbers are signif- 
icantly improved by projecting the raw correlation functions onto each j- 
scctors. The improvement is particular noticeable whenever strong anisotropics 
are present in the system, as in the case of channel flows [34, 138] ; 

(2) Anisotropic sectors with j >2 (inaccessible in present experimental data) 
possess good scahng laws [15,36]; 

(3) The scaling exponents are discrete and increasing as a function of j. 

(4) The exponents are anomalous; i.e. they differ from the dimensional pre- 
diction (46). 

(5) There exists preliminary evidence that also for j > 2, the anomalous 





(148) 



|<|p|<2A; 



104 



exponents are universal, i.e. the scaling properties are independent of the 
external forcing mechanism [131]. 

DNS were performed both in wall-bounded flows and in homogeneous (but 
anisotropic) turbulence. In wall-bounded flows the anisotropies are accompa- 
nied by inhomogeneous effects. The presence of such effects may spoil the very 
meaning of scaling, and the SO (3) decomposition should be supplemented by 
some tool to project on the homogeneous components. Otherwise, the SO (3) 
decomposition must be used carefully, and locally, only in those regions of the 
tested flow where inhomogeneous effects are confined mostly to large-scales 
[34, 138]. In the second part of this section, we discuss numerical experiments 
built such as to have a perfectly homogeneous and anisotropic statistics at all 
scales. One such example is homogeneous shear flows [17,38]. More recently, 
other homogeneous anisotropic flows have been invented and simulated, in par- 
ticular the random-Kolmogorov-flow [15, 36, 139] and a convective cell with an 
imposed linear mean profile of temperature [131]. 

7.1 Anisotropic and inhomogeneous statistics: Channel flows 

In this section, we discuss the analysis of a DNS of a channel flow using 
the SO (3) decomposition. The coordinates are chosen such that x, y and z 
are the stream-wise, span-wise, and wall-normal direction respectively. The 
simulation was done on a grid with 256 points in the stream-wise direction 
and 128 x 128 points in the two other directions. The boundary conditions were 
periodic in the span-wise and stream-wise directions and no-slip on the walls. 
The Reynolds-number based on the Taylor micro-scale was quite moderate, 
Rx 70 at the center of the channel {z — 64). The simulation was fully 
symmetric with respect to the central plane. For more details on the averaged 
quantities and on the numerical code, see refs. [34, 140, 141]. 

The analysis focused on longitudinal second-, forth- and sixth-order structure 
functions: 

S(^\r^, r) = {[6ui{r^, r)f) , Sue{r^, r) = r ■ [u{r^ + r,t)- u{r' - r, t)] . 

The r'^ coordinate specifies the location of the measurement (i.e., the center 
of mass of the two measurement points), and 2r is the separation vector. 
Previous analysis of the same data-base [140] as well as of other DNS [142] and 
experimental data [9, 72] in anisotropic flows found that the scaling exponents 
of energy spectra, energy co-spectra and of longitudinal structure functions 
exhibit strong dependence on the position r,,. For example, in [18] the authors 
studied the longitudinal structure functions at fixed distances from the walls: 

S^''\r,z) = {{u^{x + r,y,z) -u^{x,y,z))''). 
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where (■ ■ ■)z denotes a spatial average on a plane at a fixed height z.l<z< 
64. For this set of observables they found that: (i) These structure functions 
did not exhibit clear scaling behavior as a function of the distance r. Conse- 
quently, one needed to resort to Extended-Self-Similarity (ESS) [51] in order 
to extract a set of relative scahng exponents Q'^\z) = (^"\z)/(^^\z); (ii) the 
relative exponents, Q^\z) depended strongly on the height z. Moreover, only 
at the center of the channel and very close to the walls the error bars on the 
relative scaling exponents extracted by using ESS were small enough to claim 
the very existence of scaling behavior in any sense. Similarly, an experimental 
analysis of a turbulent flow behind a cylinder [72] showed a strong dependence 
of the relative scaling exponents on the position behind the cylinder for not 
too big distances from the obstacle, i.e. where anisotropic efi'ects may still 
be relevant in a wide range of scales. In the following we present an inter- 
pretation of the variations in the scaling exponents observed in non-isotropic 
and non-homogeneous flows upon changing the position in which the analy- 
sis is performed. In particular, we will show that decomposing the statistical 
objects into their different (j, m) sectors rationalizes the findings, i.e. scaling 
exponents in given (j, m) sector appear quite independent of the spatial lo- 
cation; only the amplitudes of the SO (3) decomposition depend strongly on 
the spatial location. The analysis showed three major results. The first was 
the vast improvement in scaling behavior of the structure functions as a result 
of the decomposition. A typical example is found in Fig. (15 ) where the raw 
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Fig. 15. Log-log plot of the isotropic sector of the 4th order structure function 
S'Q^Q(r), vs. r at the center of the channel = 64 (+). The data represented by (x) 
correspond to the raw longitudinal structure function, S^'^\r^ = 64, rx) averaged 
over the central plane only. The dashed line corresponds to the intermittent isotropic 
high-Reynolds numbers exponents (^^^ = 1.28. 
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fourth-order structure function, evaluated on the central plane, is compared 
to its j = component. Without the S0(3) decomposition there is no scaling 
behavior at all and one needs ESS to estimate the scaling exponents. On the 
other hand, the j — component of the structure function shows a clear scal- 
ing behavior with the expected exponent, Co^^ = 1.28. This strengthens the 
foliation hypothesis, according to which, the raw structure function is a super- 
position of power laws from different sectors of the SO (3) group. Such a sum 
looses its scale invariance once the weights of the different exponents are of the 
same order and the inertial range is small. In such cases, one needs the SO (3) 
decomposition to isolate the different sectors and retain scale invariance. 
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Fig. 16. Logarithmic local slopes of the ESS plot of raw structure 
function, ^^^^^Ll^^^^, of order 4 versus raw structure function of order 2 at rf^ = 64 

dlog{b'-^> (r,rc) ^ 



(x), at = 32 (*) and of the j = projection, 



^^^^^4§r^, centered at r,- 



64 (+), 



and at = 32(n). Also two horizontal lines corresponding to the high- Reynolds 
number limit, 1.82, and to the K41 non-intermittent value, 2, are shown. 



A second prominent result is the apparent universality of the isotropic expo- 
nents. To show this in [34] the local slopes, '^'"^^'^no ('') ^ ^j-^^ jrjgg curves of 

dlog(Sl;J (r) 

the isotropic forth-order structure-function versus the isotropic second-order 
structure function were calculated at varying the distance from the wall. De- 
spite their different locations, all curves show the same ESS slope 1.82, which 
is the expected (anomalous) value. In Fig. (16) one picture is presented for the 
logarithmic local slopes at two different distances from the channel boundary. 
To appreciate the improvements in scaling and universality, also the slopes 
of the ESS on the raw structure functions are presented Finally, the analysis 
provided another evidence that the j — 2 scaling exponent of the second order 
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structure function is about 4/3, which is the dimensional theoretical prediction 
given in (27) (see also [8, 15, 90-92]). Considering the relatively low Reynolds- 
number and the fact that the prefactors aj^m in the SO (3) decomposition (145) 
are non- universal, together with the experimental result reported in [10,29- 
31], these findings give strong support to the view that the scaling exponent in 
the j = 2 sector is universal. Before concluding this section we cite that SO (3) 
and SO (2) decomposition have also been exploited in the analysis of channel 
flow data to highlights the importance of structures as streaks and hairpin 
filaments typical of many wall bounded flow [138]. Preliminary investigation 
of the importance of SO (3) decomposition to evaluate the performance of sub 
grid models used in Large Eddy Simulations [143] have also been reported in 
[139]. 



7.2 Anisotropic-homogeneous flows 

Direct numerical simulations offer the unique opportunity to study the physics 
of anisotropy in ideal situations, that is in perfectly homogeneous fiows. Re- 
cently, considerable effort has been spent on simulating a Random- Kolmogorov- 
Flow (RKF) [15,36, 139]. The RKF is fully periodic, incompressible and with 
anisotropic large-scale energy injection. A convenient choice for the forcing is 
/ = (0, 0, f,{x)) with f,{x) = Fi cos[27rx/L^ + 0i(t)] + cos[4ttx/L^ + 02(t)], 
with constant amplitudes Fi^2 and independent, uniformly distributed, 5- 
correlated in time and with random phases 0i,2(^)- The random phases lead to 
a homogeneous statistics. To give a first validation of the statistical properties 
of the RKF fiow we plot in fig. 17 the instantaneous energy spectrum, 

E{k)^ J {u{q)-u*{q))dq. 

\q\=k 

It exhibits a scaling law in close agreement with the K41 isotropic behav- 
ior k~^/^. Also purely anisotropic quantities as the co-spectra (148), show a 
good agreement with the Lumley DNS of the RKF were reported in 

[15, 36, 139]. The resolution was 256'^ reaching Re\ ~ 100, collecting up to 70 
eddy turn over times. A long time average is necessary because of the for- 
mation of persistent large-scale structures inducing strong oscillations of the 
mean energy evolution. This is typical to many strongly anisotropic fiows. The 
viscous term was replaced by a second-order hyper-viscosity, — z/A^n. Thanks 
to both the high degree of homogeneity and to the high number of indepen- 
dent samples, a quantitative analysis of scaling laws of longitudinal structure 
functions up to the anisotropic sector j = 6 and up to order n = 6 was possi- 
ble. In other words, the longitudinal structure functions could be decomposed 
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Fig. 17. Log-log plot of instantaneous energy spectrum in the isotropic sector E{k) 
(top). The straight line is the reference isotropic k~^^^ power law. Instantaneous 

co-spectrum Eyz{kz) (bottom). Here the straight line gives the reference kz 
anisotropic Lumley prediction. The two spectra have been shifted along the ver- 
tical direction for the sake of presentation. 

according to 



6 j 

E E 

j=0 m=-j 



for n < 6. 



In fig. 18 we present the results for the isotropic sector. Here we compare the 
raw structure functions in the three directions with the projection S^{r), and 
their logarithmic local slopes (inset). Only for the projected correlation it is 
possible to measure (with 5% of accuracy) the scaling exponent by a direct 

(2) 

log- log fit versus the scale separation. The best fit gives Q = 0.70 ± 0.03. 
On the contrary, the undecomposed structure functions are overwhelmed by 
the anisotropic effects present at all scales, and the scahng law is completely 
spoiled. We stress the accuracy of these results; already at these modest 
Reynolds numbers it is possible to ascertain the isotropic scaling laws if the 
anisotropic fluctuations are disentangled properly. 

In figure 19 there is an overview for the second order structure functions in 
all the sectors, isotropic and anisotropic, for which the signal-to-noise ratio 
is high enough to ensure statistically stable results. Sectors with odd j are 
absent due to the parity symmetry of the longitudinal structure function. We 
conclude from figure 19 a clear foliation in terms of the j index : sectors with 
the same j but different m exhibit very close scaling exponents. In Table 3 
the measured exponents are compiled, showing the best power law fits for 
structure functions of orders n = 2,4,6. We stress again the discreteness 
and monotonicity of the scaling exponents as assumed in Eq. (145); there 
is no saturation of the exponents as a function of j. Second, the measured 
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Fig. 18. Analysis on the real space. Log-log plot of Sqq (r) versus r (top curve), and 
of the three undecomposed longitudinal structure functions in the three directions 

(2) 

x,y,z (three bottom curves). The straight line gives the best fit slope Q = 0.7. 
Inset: logarithmic local slopes of the same curves in the main body of the figure 
(same symbols). Notice that only the projected curve shows a nice plateau 




10 100 
log(r) 



Fig. 19. Log-log plot of the second order structure function in all sectors with a 
strong signal. Symbols refer to sectors (j,m) as follows: (0, 0), (-I-); (2, 2), (x); (4, 0), 
(□); (4,2), (★); (6,0), (o); (6,2), (■). The statistical and numerical noise affecting 
the S0(3) projection is estimated as the threshold where the j = 6 sector starts to 
deviate from the monotonic decreasing behavior, i.e. 0(10"'^) 
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exponents in the sectors j = 4 and j = 6 are anomalous, i.e. they differ from 
the dimensional estimate given in 46. Unfortunately, from the RKF data it 



n 


i = o 


J =2 


j = 4 


j = 6 




Ci") ^ n/3 


c5"^-(n + 2)/3 


Cf^-(n + 4)/3 


Ci"^-(n + 6)/3 


2 
4 
6 


0.70 (2) — 0.66 
1.28 (4) — 1.33 
1.81 (6) — 2.00 


1.1 (1) — 1.33 
1.6 (1) — 2.00 
2.1 (1) — 2.33 


1.65 (5) — 2.00 
2.25 (10) — 2.66 
2.50 (10) — 3.33 


3.2 (2) — 2.66 
3.1 (2) — 3.33 

3.3 (2) — 4.00 


Tab] 


e 3 



Summary of the numerical and experimental findings for the scaling exponents in 
the isotropic and anisotropic sectors. The values for the anisotropic sector j = 2 are 
taken from the experiments [10,31]. For the values extracted from the numerical 
simulation (columns j = 0, 4, 6), error bars are estimated from the oscillation of the 
local slopes [36, 139]. For the experimental data the error is given as the mismatch 
between the two experiments. For all sectors we also give the dimensional estimate 
Cj") = (n + i)/3 [15]. 

was not possible to obtain clean results for the j = 2 sector. This is because 
of the presence of an annoying change of sign in the projections S2m{f) for 
any m (and any order n). Still, the overall consistency of the foliation and 
hierarchical organization of scaling exponents can be checked by collecting 
the scaling exponents in the j = 2 sector from the two sets of experiments 
[10, 31] previously discussed. In Fig. 20 we show both numerical data and the 
experimental values as extracted from [10,31]. The resulting picture is fully 
coherent : experimental data coming from the j — 2 sector fit well in the global 
trend. As one can see from table 3 all the anisotropic sectors show anomalous 
scaling laws. 

7.2.1 Universality of Anisotropic Fluctuations 

The third numerical experiment that we discuss here is devoted to study uni- 
versal properties of anisotropic scaling. We have already commented that there 
is a nice qualitative and quantitative agreement between the values extracted 
for the j — 2 sector (the only one available from experimental data) from 
different experiments. In order to check whether this universality holds also in 
higher anisotropic sectors one has to rely on DNS. In [131] a first direct com- 
parison between anisotropic scaling of longitudinal structure functions from 
two different homogeneous systems, the RKF and a homogeneous Rayleigh- 
Benard convective flow was reported. 

A Homogeneous Rayleigh-Benard system is a convective cell with fixed linear 
mean temperature profile along the vertical direction. The flow is obtained 
by decomposing the temperature field as the sum of a linear profile plus a 
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Fig. 20. Scaling exponents, Q > of structure functions of order n = 2,4,6 for 
isotropic and anisotropic sectors. Prom the DNS of RKF we have : isotropic sec- 
tor, J = (+); anisotropic sectors, j = 4 (x) and j = 6 (★). From the experimental 
data [10, 31], we have j = 2 (□). For an estimate of error bars see Table 3. 

fluctuating part, T(x,y, z;t) = T'{x,y, z;t) + {AT/2 — zAT/H), where H is 
the cell height and AT the background temperature difference. The evolution 
of the system can be described by a modified version [144] of the Boussinesq 
system [145]: 



dtu + (m • V) w = -Vp + i/V^w + agT'z 

AT 

dtT' + (w • V) T' = kV^T' - —V,. 

H 

where a is the thermal expansion constant, v and k the kinematic viscosity 
and the thermal difFusivity coefficients, and g is the acceleration due to gravity. 
In [131] fully periodic boundary conditions were used for the velocity field, u, 
and temperature, T', fields. 

Anisotropic effects in the Rayleigh-Benard system were analyzed in [131] start- 
ing from the stationary equation for the second order velocity structure func- 
tions; the extension of Karman-Howarth equation in the presence of a buoy- 
ancy term [146]. The result is, neglecting for simplicity tensorial symbols: 

{5u{vf)^ er + agzr ■ {5T{r) Su{r)) (149) 

j=0,l,... j=0 j=l ® j=l,2,.. 

where e denotes the energy dissipation, (^^(r)^) and (5T(r) ^^(r)), the general 
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third-order velocity correlation and temperature-velocity correlation, respec- 
tively. In Eq. (149) for each term the value of its total angular momentum, 
j, is indicated. Notice that the energy dissipation term in (149) has a non- 
vanishing limit, for high Re, only in the isotropic sector, j = 0. On the other 
hand, the buoyancy coupling, agz, brings only angular momentum j — 1. 
Due to the usual rule of composition of angular momenta we have that the 
buoyancy term, agz ■ (5T(r) 6u{r)), has a total angular momentum given by: 
jtot = 1®J = J) J + Using the angular momenta summation rule for j, 

one can decompose the previous equation obtaining the following dimensional 
matching, in the isotropic sector: 

{Su{r)^) j=o ~ e r -|- agzr{Su{r) ST{r))j=i + . . . 

and in the anisotropic sectors, j > 0: 

{5u{r)^)j ~ ag zr {5u{r) 5r(r))(j_i) + 



(150) 



where sub-dominant contributions coming from the j and j + 1 sectors of 
{6v{r) 5T{r)) are neglected. 

In the isotropic sector the buoyancy term is sub-dominant with respect to the 
dissipation term at scales smaller than the Bolgiano length. Lb = {e)^^'^N~^/'^{ag)~^^'^ 
where N is the rate of temperature dissipation. This is the case for the nu- 
merical simulation presented in [131], where velocity fluctuations are closer 
to the typical Kolmogorov scaling, Su{r) ~ r^^^, rather than to the Bolgiano- 
Obukhov scaling [1], du{r) ~ r^/^. 

Regarding the anisotropic sectors, Eq.(150) is the dimensional prediction for 
the system, consistent with the anisotropic properties of the buoyancy term, 

sector by sector. 

In [131] the S0(3) decomposition was applied in this system to velocity struc- 
ture functions (44) and to objects 

G'(^'^)(r) = ([(tx(r) -um ■ {T {r) - T (0))) ^Y.G%'\r) 

jm 



The dimensional matching of Eq. (150) can be extended to any order, giving: 

'5S(r)-rG(r,%)(r). 



Denoting with Xj^'^^ the anisotropic scaling exponents of the buoyancy terms, 
^jrra '* (^) ~ ''"^^ 'we get the dimensional estimate: 



— i + xf-i'^\ (dimensional prediction) (151) 
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In [131] it was shown that this dimensional prediction is not obeyed; the 
exponents cj"'' appear to be systematically smaller than the prediction (151). 
Interestingly, enough the log-log plots computed in the sectors j = 4, 6 show 
a good qualitative agreement with those calculated in the RKF of ref [36] 
as can be seen in fig. 21 where we compare the projection on the j — 4: 
sector of structure functions of different orders. Similar results are obtained 




10 r 100 



Fig. 21. Log-log plot of compensated anisotropic j = 4, m = projections 
S4'o{r)/r'^*' vs r, for both HRB and RKF flows. Top curves refer to p = 2: the 

' (2) 

best fit exponents which compensate HRB and RKF curves are Q ' = 1.7 and 
Q = 1.66, respectively. Curves in the middle refer to the same quantities but for 
p = 4: compensation has been obtained with ([^^ = 2.05 for HRB, and (^^^ = 2.2 
for RKF. Bottom curves refer to p = 6: here = 2.3 for HRB, and (^^^ = 2.5 for 
RKF. Notice that the curves of the two flows are compensated with very similar 
values of the exponents (within 10%). Inset: the same but for j = 4, m = 2, com- 
pensation has been done with the same values used for j = A,m = 0, to show the 
independence of the scaling exponents from the choice of the reference axis labeled 
with m . 

for j — 6 sector. These preliminary findings, if confirmed by other independent 
measurements, would support universality for anisotropic scaling exponents in 
three dimensional turbulence. 



7.3 Scaling of Longitudinal and Transversal structure functions 

As discussed in Sect. 2.6 there exists experimental and numerical data sug- 
gesting that longitudinal and transversal structure functions in supposedly 
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isotropic flows show different scaling exponents [30, 62-65]. One needs to clearly 
distinguish between experimental and numerical data. The former can never 
be considered fully isotropic; the best one can do is to try to perform a multi- 
fit procedure to clean out sub-leading anisotropic contributions as already 
explained in details in Sect. 6. This fitting procedure is, of course, affected by 
experimental errors which cannot be eliminated. Therefore it is quite danger- 
ous to make any firm conclusion about supposed different scaling exponents 
of longitudinal and transversal isotropic scaling on the basis of only exper- 
imental data . Numerical data are not much safer. Here anisotropy can be 
much better controlled. With isotropic forcing the only source of anisotropy 
is the 3-dimensional grid whose effect is usually too small to explain possi- 
ble discrepancies between longitudinal and transversal scalings. Indeed some 
state-of-the-art isotropic DNS indicate the possibility of different scaling ex- 
ponents both for inertial range structure functions [64] and for coarse grained 
energy and enstrophy measures [62, 147]. In Table 4 we summarize the best- fit 
values of the scaling exponents measured in [64]. The small scale fiuctuations 
were probed [147] by comparing the scaling of the coarse grained energy dissi- 
pation over a box of size r, e(r) (Eq. 17), and of the coarse grained enstrophy 
dissipation: a;(r, x) — ^ I\y\<r '^(^^ + v) where u!{x) is the local enstrophy 
dissipation. Different scaling exponents were measured for the averaged quan- 
tities, ((e)^(r)), (a;^(r)). Being scalar quantities one expects that in isotropic 
ensembles they would not have different exponents. From the theoretical point 
of view, different scaling exponents of longitudinal and transverse structure 
functions in isotropic ensembles are unlikely. In the language of the SO (3) 
decomposition it amounts to the scaling exponents depending on the g'-index 
which labels different basis functions with the same rotation properties. In 
the exactly solvable models examined before, this had never happened. In 
general one would need a different symmetry to lift the degeneracy of differ- 
ent q dependent basis functions. At this point this problem remains somehow 
unsettled. New numerical tests on larger grids and/or with a better resolved 
viscous behavior are needed before a firm statement can be made. 
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6 


8 


10 


An) 

So 


0.701 (14) 


1.29 (3) 


1.77 (4) 


2.17 (7) 


2.53 (9) 


An) 

So 


0.709 (13) 


1.27 (2) 


1.67 (4) 


1.93 (9) 


2.08 (18) 



Table 4 



Measured values of the longitudinal (first raw) and transverse (second raw) scaling 
exponents at Rex = 460 taken from [64]. One should note that the scaling range 
displayed by the scaling plots in [64] are relatively short, indicating that finite Re 
effects may still be rather important 
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l.Jf. Anisotropies in decaying turbulence 



Decaying turbulence has attracted the attention of various communities and 
is often considered in experimental, numerical and theoretical investigations 
[1,5, 148]. It is in fact quite common that even experiments aimed at stTidying 
stationary properties of turbulence involve processes of decay. Important ex- 
amples are provided by a turbulent flow behind a grid (see [149] and references 
therein) or the turbulent flow created at the sudden stop of a grid periodi- 
cally oscillating within a bounded box [150] . In the former case, turbulence is 
slowly decaying going farther and farther away from the grid and its charac- 
teristic scale becomes larger and larger (see [149] for a thorough experimental 
investigation). Whenever there is sufficient separation between the grid-size 
Lin and the scale of the tunnel or the tank Lq 3> a series of interest- 
ing phenomenological predictions can be derived. For example, the decay of 
the two-point velocity correlation function, for both isotropic and anisotropic 
flows, can be obtained under the so-called self-preservation hypothesis (see 
[1] chapter XVI). That posits the existence of rescaling functions allowing to 
relate correlation functions at different spatial and temporal scales. By in- 
serting the rescaling function into the equations of motion, asymptotic results 
can be obtained both for the flnal viscosity-dominated regime and for the 
intermediate asymptotic when nonlinear effects still play an important role. 

Here, we review some recent attempts to investigate the decay of three-dimensional 
homogeneous and anisotropic turbulence by direct numerical simulations of 
the Navier-Stokes equations in a periodic box [151] for both short and large 
times. The initial conditions are taken from the stationary ensemble of the 
Random Kolmogorov Flow discussed in the previous subsection. Here the cor- 
relation length-scale of the initial velocity field Lj„ is of the order of the size 
of the box Lo ^ L^n- 

On the one hand, one is interested in the long time decay regime where the typ- 
ical interesting questions are: (i) how do global quantities, such as single-point 
velocity and vorticity correlations, decay? (ii) What is the effect of the outer 
boundary on the decay laws ? (iii) Do those quantities keep track of the initial 
anisotropy ? (iv) As for the statistics of velocity differences within the inertial 
range of scales, is there a recovery of isotropy at large times ? (v) If so, do 
strong fluctuations get isotropic at a faster/slower rate with respect to those of 
average intensity ? (vi) Do anisotropic -and isotropic- fluctuations decay self- 
similar ly ? (vii)If not, do strong fluctuations decay slower or faster than typical 
ones ? On the other hand, the interest in the early stages of the decay is led by 
a hope of establishing a link between the small-scale velocity statistics in this 
phase and in forced turbulence. If such links existed, they would shed addi- 
tional light on the universality of forced turbulence. As turbulence decays, the 
effective Re decreases, while the viscous characteristic scale and time increase. 



116 



In [151] an offline analysis at fixed multiples {0, 1, 10, 10^, 10^ 10^ 10^ 10^} tq 
of the initial large-scale eddy turnover time Tq = Lq/u^^^^ was performed. 
A first hint on the restoration of isotropy at large times can be extracted from 
the analysis of single point quantities as 



Eii = Ui{x, t)ui{x, t), Qii = u;i{x, t)u;i{x, t). 



Here with we denote the average over spatial coordinates only, whereas 
(• • •) indicates the average over both initial conditions and space. The sym- 
metric matrices Eii{t) and ^lii{t) can be diagonalized at each time-step and 
the eigenvalues Ei(t), E2(t), E^it) and fii(t), 1^2 (t), ^s^t) can be extracted. The 
typical decay of Ei{t) and for i = 1, . . . , 3 is shown in Fig. 22. During the 




Fig. 22. Log- log plot of the eigenvalues of energy and vorticity matrices vs. time, 
expressed in tq unit. 

self-similar stage, t G [10, 10®], the energy eigenvalues fall off as £'{1,2.3} ~ t^'^, 
as expected for the decay in a bounded domain [149,152]. The enstrophy 
eigenvalues, ^^{1,2,3} decay as t~^^^^ as predicted from a simple dimensional 
argument [151]. To focus on the process of recovery of isotropy in terms of 
global quantities one may track the behavior of two sets of observables 

which vanish for isotropic statistics. Their rate of decay is therefore a direct 
measurement of the return to isotropy. The energy matrix En is particularly 
sensitive to the large scales while small-scale fluctuations are sampled by flu. 
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Fig. 23. Log-log plot of the anisotropy content at the large scales (Ai2-E'(t), Ai^Eit), 
top curves) and the small scales [Au^it], Aisfl{t), bottom curves) as a function of 
time, expressed in tq unit. The large-scale (small-scale) anisotropy content is defined 
as the mismatch between the eigenvalues of the single-point velocity (vorticity) 

correlation. 

As seen from Fig. 23, both large and small scales begin to isotropize after 
roughly one eddy turnover time and become fully isotropic (within statistical 
fluctuations) after 100 eddy turnover times. However, small scales show an 
overall degree of anisotropy much smaller than the large scales. 

Concerning small scales properties, in [151] a simple anisotropic generalization 
of the self-preservation hypothesis (see, e.g. Ref. [5]) was proposed: 



Here with Vj^ (t) wc take explicitly into account the fact that large-scale 
velocity properties may depend in a nontrivial way on both (j, m) and the 
order n. Furthermore, Lj^{t) accounts for the possibility that the characteris- 
tic length scale depend on the S0(3) sector. In analogy with the observations 
made in the stationary case [10, 15, 29, 30, 34-36, 55] a scaling law was postu- 
lated: 



The time behavior is encoded in both the decay of the overall intensity, ac- 



sSir,t)^vSit)fi:;i{r/LjUt)). 




(152) 
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counted by the prefactors a^-^(t), and the variation of the integral scales 
Lj„i{t). The representation Eq. (152) is the simplest one fitting the initial 
time statistics for t = and agreeing with the evolution given by the self 
preservation hypothesis in the isotropic case. The power law behavior for 
fjm {^1 Ljm{t)) can be expected only in a time-dependent inertial range of 
scales ri{t) <S r <S Li{t). As for the exponents appearing in (152), their val- 
ues are expectedly the same as in the stationary case. Concerning the time 
evolution, it seems difficult to disentangle the dependence due to the decay 
of o!jll{t) from the one due to the growth of the integral scale Ljm{t). The 
existence of a running reference scale, Ljm{t) introduces some non-trivial re- 
lations between the spatial anomalous scaling and the decaying time proper- 
ties, and those relations might be subject to experimental verification. In the 
case discussed in [151], the fact that the initial condition has a characteristic 
length-scale comparable with the box size, simplifies the matter. Indeed we 
expect that Ljm{t) ^ Lq, and the decay is due only to the fall-off of aj^{t). 
An obvious shortcoming is that the width of the inertial range LQ/rj{t) shrinks 
monotonically in time, thereby limiting the possibility of precise quantitative 
statements. 



7.4-1 Long time decay 

A quantitative way to define the temporal rate of recovery of isotropy at a 
fixed scale in the inertial range is given by the dimensionless ratio: 



nW(r,t)^5£M^H"\ (153) 



^jm V' > V — (n) 



In Fig. 24, we plot 11^^ (r, i) at r = 80 for structure functions of order 
n — 2,4,6 and for the most intense anisotropic sector, (j, m) = (4,0). All 

anisotropic sectors, for all measured structure functions, decay faster than the 
isotropic one. The measured slope in the decay is about S^"'' ~ 0.3 for all 
n, within the statistical errors. Note that these results agree with the simple 
picture that the time-dependence in (152) is entirely carried by the prefactors 
aj^(t) and the value of the integral scales Ljm{t) is saturated at the size of 
the box. Indeed, by assuming that large-scale fiuctuations are almost Gaus- 
sian we have that the leading time- dependence of a^-^^ is given by Oj^aoo" 
For the isotropic sector, a'^^'' ~ {cLm)"'i plugging that in (153), one get: 
n^m (^)^) ~ (^fmi't) / i't) ~ with S ~ 0.3(±0.1) independent of n. The 
quality of data is insufficient to detect possible residual cfi^ccts due to Ljm{t), 
which would make depend on n and j because of spatial intermittency. 
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Fig. 24. Hierarchical organization of anisotropic fluctuations at long times. Log- log 
plot of the anisotropic projections normalized by the corresponding isotropic pro- 
jection (see text), at two fixed scales r = 80 and r = 40 (inset) for n = 2,4, 6 in 
the anisotropic sector j = 4, m = 0. Symbols read as follows: n4Q^ (full box); n4Q^ 
(star); II^q^ (empty box). The straight line is with S ~ 0.3. Same symbols in 
the inset. 

Let us denote with V{A, r; t) the probability to observe a given longitudinal 
fluctuation, 6ui{r,t) = A in the direction r at a given time, t. For any given 
fixed value A and for any given time, t, we can project ^(A, r; t) on the S0(3) 
basis functions: 

ViA, r;t)^Yl E T^jm (r, A; t) Yj^r), (154) 

j=0 m=-j 

where now the projection, Vjm {f^, A; t) play the role of an effective PDF 
for each SO (3) sector. The projection of any longitudinal structure function, 
S'*^"^(r,t) on any sector, (j, m) can be reconstructed from the corresponding 
projection of the PDF on the same sector, Vjm (r. A; t), by averaging over all 
possible A: 

Sfl (r, t)^ j dAA^Vjm (r. A; t) 
which establish the link between decomposition (44) and (154). 

The interesting fact that the decay properties of the anisotropic sectors are 
almost independent of n indicates that a non-trivial time dependence in the 
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shape of the PDF's Vjm (f, A; t) for j' > must be expected. The most accurate 
way to probe the rescahng properties of Vjm {i^, A; t) in time is to compute the 
generahzed flatness: 



cr(n) 
jm 



ir,t) 



(n) 
J 



Were the PDF projection in the {j,m) sector self-similar for t 3> to. then 
^jm t) would tend to constant values. This is not the case for anisotropic 
fluctuations, as it is shown in Fig. (25). The curves Kj^{r,t) are collected for 
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Fig. 25. Log-log plot of the generalized flatness, xj"] (r, t) of order n = 4, 6 for both 
the isotropic (two bottom curves) , and the anisotropic sector (j = 4, m = 0) (two 
top curves) at r = 80, and as a function of time. In the inset we plot the same 
quantities, in the same order, at a different inertial range scale, r = 40. 

two fixed inertial range separations, r — SO and r = 40 (inset), for two dif- 
ferent orders, n — 4,6 and for both the isotropic and one of the most intense 
anisotropic sectors (j = 4, m = 0) . The isotropic flatness tends toward a con- 
stant value for large t. Conversely, its anisotropic counterparts are monotoni- 
cally increasing with t, indicating a tendency for the anisotropic fluctuations 
to become more and more intermittent as time elapses. Also the behavior 
in Fig. 25 is in qualitative agreement with the observation previously made 
that all the time dependence can be accounted for by the prefactors a^-^(t). 
Indeed, assuming that the length scales Ljm{t) have saturated and that the 
large scale PDF is close to Gaussian, it is easy to work out the prediction 
Kjmir,t) - r=(^-"/2)^ ^ 2(n/2 - 1). We conclude this section by a 
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brief summary of the results. It was found that isotropic fluctuations persist 
longer than anisotropic ones, i.e. there is a time-recovery, albeit slower than 
predicted by dimensional arguments, of isotropy during the decay process. It 
was also found that isotropic fluctuations decay in an almost self-similar way 
while the anisotropic ones become more and more intermittent. Qualitatively, 
velocity configurations get more isotropic but anisotropic fluctuations become, 
in relative terms, more "spiky" than the isotropic ones as time elapses. 

l.Ji..2 Short-time decay 

It is interesting to notice that it is possible to relate the small-scale universal 
properties of forced turbulent statistics to those of short-time decay for an en- 
semble of initial configurations [151]. As already remarked, one cannot expect 
an universal behavior for all statistical observables, as the very existence of 
anomalous scaling is the signature of the memory of the boundaries and/or 
the external forcing throughout all the scales. Indeed, the main message we 
want to convey here is that only the scaling exponents of both isotropic and 
anisotropic small-scale fiuctuations are universal, at least for forcing concen- 
trated at large scales. The prcfactors are not expected to be so. There is 
therefore no reason to expect that quantities such as the skewness, the kurto- 
sis and in fact the whole PDF of velocity increments or gradients be universal. 
This is the same situation that we discussed in great details in previous sec- 
tions for the passive transport of scalar and vector fields. However, carrying 
over the analytical knowledge developed for linear hydrodynamic problems 
involve some nontrivial, yet missing, steps. For the Navier-Stokcs dynamics, 
linear equations of motion appear when we consider the whole set of correla- 
tion functions as discussed in Sect. 4.3. These equations can be rewritten in a 
schematic form: 

a^C^") = r("+^)C("+^) + i/D(")C(") + (155) 

where r*^"'+^) is the integro-differential linear operator coming from the inertial 
and pressure terms, C^"-* is a shorthand notation for a generic (n)-th order cor- 
relator and is the linear operator describing dissipative effects. Finally, 
F^") is the correlator involving increments of the large-scale forcing / and 
of the velocity field. The balance between inertial and injection terms cannot 
lead to anomalous scaling. A natural possibility is that a mechanism similar to 
the one identified in linear transport problems be at work in the Navier-Stokes 
case as well. The anomalous contributions to the correlators would then be as- 
sociated with statistically stationary solutions of the unforced equations (155). 
The scaling exponents would a fortiori be independent of the forcing and thus 
universal. As for the prefactors, the anomalous scaling exponents are positive 
and thus the anomalous contributions grow at infinity. They should then be 
matched at the large scales with the contributions coming from the forcing to 
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ensure that the resulting combination vanish at infinity, as required for corre- 
lation functions. The aim here is not to prove the previous points but rather to 
test whether they fail: the Navicr- Stokes equations, being integro-differential 
and non-local, might directly couple inertial and injection scales and spoil the 
argument. This effect might be particularly relevant for anisotropic fiuctu- 
ations where infrared divergences may appear in the pressure integrals (see 
Sect. 5.3). In order to investigate the previous point, we performed two sets 
of numerical experiments in decay. 

The first set. A, is of the same kind as in the previous section, i.e. we in- 
tegrated the unforced Navier-Stokes equations with initial conditions picked 
from an ensemble obtained from a forced anisotropic stationary run. Statisti- 
cal observables are measured as an ensemble average over the different initial 
conditions. The ensemble at the initial time of the decay process therefore 
coincides with the stationary state in forced runs. If correlation functions are 
indeed dominated at small scales by statistically stationary solutions of the 
unforced equations then the field should not decay. Specifically, the field should 
not vary for times smaller than the large-scale eddy turnover time tq. Those 
are the times when the effects of the forcing terms start to be felt. Note that 
this should hold at all scales, including the small ones whose turnover times 
are much faster than tq. 

The second set of numerical simulations (set B) takes the same initial condi- 
tions but for the random scrambhng of the phases : Wj(k) — > Pji(k)w;(k) exp(i^;(k)), 
with ^i(k) i.i.d. random variables. In this way, the spectrum and its scaling 
exponent are preserved but the wrong organization of the phases is expected 
to spoil the statistical stationarity of the initial ensemble. As a consequence, 
two different decays are expected for the two sets of initial conditions. In par- 
ticular, contrary to set A, set B should vary at small scales on times of the 
order of the eddy turnover times Tr ~ r^^^. This is exactly what has been 
found in the numerical simulations for both isotopic and anisotropic statistics 
as can be seen for the anisotropic case in Fig. 26, where the temporal behav- 
ior of longitudinal structure functions of order 2 and 4 is shown. The scahng 
exponents of the contributions responsible for the observed behavior at small 
scales are thus forcing independent. 

To conclude, the data presented here support the conclusion that nonlocal 
effects peculiar to the Navier-Stokes dynamics do not spoil arguments on uni- 
versahty based on analogies with passive turbulent transport. The picture 
of the anomalous contributions to the correlation functions having universal 
scaling exponents and non-universal prefactors follows. 
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Fig. 26. Top: Temporal decay of the second-order anisotropic structure function 
{r,t), rescaled by its value at t = 0. Here r = 30, inside the inertial range. The 
two curves refer to the time evolution of the structure function starting from the 
forced-stationary velocity fields (squares, set A) and from the randomly dephased 
velocity fields (circles, set B). Time is normalized by the integral eddy turnover time. 
Notice that for set B we observe changes on a time scale faster than the integral eddy 
turnover time. That is to be contrasted with the case A, where structure functions 
are strictly constant in time up to an integral eddy turnover time. Bottom: The 
same curves but for the fourth-order structure function. 

8 Concluding Discussion, 

In this review we presented a mathematical framework in which anisotropy 
in turbulence can be studied, and we have tested its utility in the context 
of experimental analysis, numerical simulations and analytical models. The 
basic idea is to express the various statistical quantities of turbulence (e.g., 
structure functions, correlation functions) in terms of their projections on the 
different sectors of the S0(3) group. 

The utility of the SO (3) decomposition should be assessed in two main aspects. 
The first aspect is its functionality as a tool for characterizing anisotropy, 
whereas the second, and deeper aspect, is its physical relevance and the the- 
oretical and analytical advantages that are gained by using it. 

As a tool for describing anisotropy, the SO (3) decomposition is probably the 
most natural and general method. It is of high resolution - it subdivides the 
observed anisotropy into different sectors - the (j, m) sectors. The weights 
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of the various sectors give us a fine resolution of the anisotropy structure. 
Instead of having one measure for anisotropy (e.g., the overall percentage of 
anisotropy), we have an infinite set of numbers that compose a detailed profile 
of the anisotropy structure. 

The SO (3) decomposition is also very general. It is applicable to any physical 
observable that has a well-defined transformation under rotations. These can 
be, for example, correlation functions, structure functions or Green's functions 
(response functions). The observables themselves may depend on any number 
of space coordinates or even be space independent. They may also be scalars, 
vectors or tensors. Any such quantity can be presented as a sum of parts that 
belong to the different (j, m) sectors of the rotation group. Additionally, since 
the SO (3) decomposition is invariant to isotropic operations, it is invariant 
to the most common operations that we use. For example, to obtain the nth 
order longitudinal structure function we can take the full nth order structure 
function (which is a tensor) and contract it with n unit vectors in the direction 
of the separation distance. Since this operation is linear and isotropic, it will 
preserve the (j, m) sectors of the full structure function. That is, the (j, m) 
sector of the full structure function will be transformed into the same (j, m) 
sector of the longitudinal structure functions. The same thing happens for 
operations such as differentiations (for example when we look at moments of 
the gradient fields), space averaging, time averaging, coordinate fusion etc... 
From the pure theoretical point of view the first and most obvious advantage 
of using the SO (3) decomposition is its elegance and the overwhelming simpli- 
fication that it offers in analytical calculations. The S0(3) decomposition may 
also have a deeper physical justification if it produces universal quantities such 
as distinct scaling behavior in the anisotropic sectors. There are several differ- 
ent evidences that suggest that this is indeed the case. Experimental results 
clearly show that a better scaling is achieved if we take higher j components 
into account. Additionally, different experimental setups seem to lead to the 
same numerical values of the anisotropic exponent. This is a strong support 
for the hypothesis that the anisotropic sectors of the structure functions have 
universal exponents. Another support for the idea that the SO (3) decomposi- 
tion exposes universal quantities, comes from numerical simulations. In DNS, 
the SO (3) decomposition can be performed directly (since the velocity field 
is accessible in every point in space and time) which makes the results much 
less ambiguous. We clearly see that even in the very moderate Reynolds num- 
bers of the simulations, a scaling behavior is detected once we use the SO (3) 
decomposition. In some cases, without the SO (3) decomposition no scaling 
behavior is seen at all. Furthermore, the resulting exponents in the isotropic 
sector are remarkably similar to the experimental values which are measured 
at very high Reynolds numbers. This is a strong indication that at least the 
isotropic sector has a universal profile, and therefore by disentangling it from 
the anisotropic sectors we get universal results. In other sectors of the rota- 
tion group, the scaling behavior is not as good, and in some sectors there is 
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no scaling at all. However, in those sectors where scaling was detected, the 
scaling exponents seem to agree with the theoretical and experimental predic- 
tions. Sectors with same j and different ms had the same scaling exponents 
(when scale- invariance was observed). And finally, all exponents increased as 
a function of n (order of the structure-function) as well as j'. It is still not 
clear whether the "bad", non scale-invariant behavior that was detected in 
some sectors, is a result of the poor Reynolds numbers of the simulations, 
or is a genuine effect that tells us that the foliation picture is incomplete. A 
further research with higher resolution is probably needed to settle this issue. 
In the Navier-Stokes case one can prove a "weak foliation" . Weak fohation is 
an approximate foliation that happens in the case of weak anisotropy, when 
we linearize the anisotropic part of the theory around its isotropic part. In 
such case, the linearized anisotropic part of the theory is subject to a linear 
and isotropic equation whose kernel contains the isotropic solution, and hence 
foliation occurs. This is a very robust approximation since it holds for virtually 
any nonlinear and isotropic theory in the case of weak anisotropy. Additionally, 
we know that as the Reynolds number increases, the statistics becomes more 
and more isotropic and therefore the linear approximation becomes better and 
better. 

The S0(3) decomposition has a physical relevance also in the presence of 
strong anisotropy. If we merely consider the direct application of the Navier- 
Stokes and continuity equations on the various structure functions, we get 
a family of linear and isotropic constraints. These constraints are valid in- 
dependently of the amount of anisotropy in the system. Their linearity and 
isotropy lead to foliation and hence one can discuss them in every sector inde- 
pendently. In some sectors, they are sufficient to determine the full solution, 
whereas in others they can reveal some general properties of the solution. 
For example, the isotropic sector of the third-order structure function is com- 
pletely determined by the different constraints, and is given by the well-known 
4/5 law of Kolmogorov. Notice that because of foliation, this is true also in 
the presence of anisotropy, which means that the 4/5 law holds also in an 
anisotropic turbulence. In the j = 2 sectors, on the other hand, the third-order 
structure-function is given by two undetermined scalar functions, whereas in 
the J = 4, 6, ... it is given by three. Another example is the j — 1 sectors in 
the second-order correlation function, which must all vanish. 

To conclude, the framework of the SO (3) decomposition provides an elegant 
and efficient way to describe anisotropy in turbulence. It also greatly simpli- 
fies many analytical calculations that involve anisotropic quantities, mainly 
through the mechanism of foliation. This mechanism is present in simplified 
models of turbulence, and may also be valid approximately in Navier-Stokes 
turbulence. It predicts that the anisotropic sectors of the statistics have uni- 
versal properties such as scaling exponents. Further research is needed to ac- 
curately measure these anisotropic exponents in experiments as well as in 
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numerical simulations. 

A quantitative computation of the anisotropic exponents in Navier-Stokes tur- 
bulence from first principles may very well be an illusive goal. Nevertheless, 
we believe that the general principles of the physics behind these exponents 
and anisotropy in general can be understood. 
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A The General Form of the 2nd Rank Tensor 

In this appendix we discuss the general structure of the second rank correlation 
functions 

F"^(r) = (m"(x + r)M'^(x)) , (A.l) 

In (A.l) homogeneity of the flow is assumed, but not isotropy. Note that this 
object is more general than the structure function S°'^ in being nonsymmetric 
in the indices, and having no definite parity. We wish to find the basis functions 
Bgj^{r), with which we can represent F"^(r) in the form: 

^°^(r) = E S.>n(rX^(f) (A.2) 

and derive some constraints among the functions aqjm{i^) that result from 
incompressibility. We shall see, that due to the isotropy of the incompressibility 
conditions, the constraints are among aqj^(r) with the same j,m only. 
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We begin by analyzing the incomprcssibility condition: An incompressible flow 
with constant density is characterized by the relation: 9q,-u"(x, t) = as a 
result, one immediately gets the following constraints on F"'^{r): 

9«F"^(r) = 0, 9^F"^(r) = 0. 

Plugging the trial tensor (A. 2) into the last two equations we obtain 2 equa- 
tions connecting the different aqjrn' 

da E a,,jm{r)Bfa^) =0 d^Yl ^<i,3m{r)B^,U^) = (A.3) 

q,jm q,jm 



We first notice that the differentiation action is isotropic. As a result, if T"'^ (r) 
is some arbitrary tensor with a definite j, m transformation properties, then 
the tensor daT°'^{r) will have the same j^m transformation properties. Com- 
ponents with different j, m are linearly independent. Therefore equations (A.3) 
should hold for each j', m separately. 

Next, we observe that (A.3) are invariant under the transformation F"'' — > 
F^"^. As a result, the symmetric and anti-symmetric parts of F°^^ should satisfy 
(A.3) independently. To see that, let us write F"^ as a sum of a symmetric 
term and an anti-symmetric term: F"^ = Fg^ + Pa^ i we then get: 



daF'^f = d^Ff + daFf = d^Fl'^ - daF^"" = 
dpF'^f^ = df^Ff + dpFf = 

from which we conclude: (9o,F^^ = daF^^ = 0. Finally, (A.3) is invariant under 
the transformation F°'^{r) — )• F"^(— r) and as a result the odd parity and 
the even parity parts of F"^ should fulfill (A.3) independently. We conclude 
that a necessary and sufficient condition for (A.3) to hold is that it holds 
separately for parts with definite j,m, definite symmetry in the a, (3 indices 
and a definite parity in r: 



a summation is over 



a5jm(r)5j„(f) = 

^q,3m '^i*'^ definite S3mimetries 



where the summation is over q such that -B^j^ has a definite indices symmetry 
and a definite parity. 

According to (37) we can write these -B^^^ as: 

(1) {—y parity, symmetric tensors: 
. <^(r)^r-^-5-/'$,„,(r). 
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• b: 



,a/3 
7,jm 



(r) 



(2) (—y parity, anti-symmetric tensors: 



(3) 



parity, symmetric tensors 



-j-i 
-i+i 



(4) (— )-''''^ parity, anti-symmetric tensors: 

In order to differentiate tliese expressions we can use tlie following identities: 



which give rise to: 

r"aa$j^(r)=j$j„,(r) . 

From this point, it is a matter of simple (though somewhat lengthy) algebra to 
derive the differential constraints among agjrni^)- The results arc as follows: 

(1) qe {1,7,9,5} 

«'ijm(^) - jr'^CLijm + ja'Tjra " fr~^a7,jm + Ogj™ + 2r"^a9j„ = , (A.4) 



r aijrn + a 



h,jm + 3r ^ayjm + (j - 1) a'^,jm - {f - 3j + 2) r ^a^jm = . 



(2) q e {3} 

«3jm - Jr'^asjm = , (A.5) 
«3,im + r'^aajm = . 

These equations have no solutions other than: asjrn{T^) — 0. 

(3) ge{8,6} 

O'sjm + ^r'^asjm + {j - l)a6jm - {f - 2i + 1) r"^a6jm = . (A.6) 

(4) q e {4, 2} 



(A.7) 
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There are obviously more unknowns than equations, since we merely exploited 
the incompressibihty conditions. Nevertheless, we believe that the missing 
equations that arise from the dynamical hierarchy of equations will preserve 
the distinction between aqjrn of different j, m (again, due to the isotropy of 
these equations). Note also, that the above analysis holds also for the second- 
order structure function 

5"'^(r) = (K(x + r) - m"(x)] [^'^(x + r) - m^(x) 



Only that in this case we should only consider the representations q = 1, 7, 9, 5 
for even j and the representations g = 8, 6 for odd j. This follows from the fact 
that S^^{r) is symmetric with respect to its indices and it has an even parity 
in r. Also, in that case, it is possible to go one step further by assuming a 
specific functional form for the aqj^(r). We know that the S"'^{y) is expected 
scale in the inertial range, and we therefore may assume: 

where Cq^jm are just numerical constants. If we now substitute this definition 
into the equations (A. 4, A. 6), we get a set of linear equations among the Cq^jm- 
These relations can be easily solved and give us two possible tensors for even 
j {q = 1,7,9,5) and one tensor form for odd j (from q = 8,6). This kind 
of approach was taken in the two-probes experiment which is described in 
Sect. 6. 



B Anisotropy in d-dimensions 



To deal with anisotropy in d-dimensions we need classify the irreducible rep- 
resentations of the group of all d-dimensional rotations, SO(d) [153], and then 
to find a proper basis for these representations. The main linear space that we 
work in (the carrier space) is the space of constant tensors with n indices. This 
space possesses a natural representation of SO(d), given by the well known 
transformation of tensors under d-dimensional rotation. 

The traditional method to find a basis for the irreducible representations of 
SO(d) in this space, is using the Young tableaux machinery on the subspace of 
traceless tensors [153, 154]. It turns out that in the context of the present work, 
we do not need the explicit structure of these tensors. Instead, all that matters 
are some relations among them. A convenient way to derive these relations is 
to construct the basis tensors from functions on the unit (i-dimensional sphere 
which belong to a specific irreducible representation. Here also, the explicit 
form of these functions in unimportant. All that matters for the calculations 
is the action of the Laplacian operator on these functions. 
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Let us therefore consider first the space Sd of functions over the unit d- 
dimensional sphere. The representation of SO(d) over this space is naturally 
defined by: 

On'ifiu) = "^{TZ-^u) , (B.l) 

where ^(■u) is any function on the d-dimensional sphere, and TZisa d-dimensional 
rotation. 

Sd can be spanned by polynomials of the unit vector il. Obviously (B.l) does 
not change the degree of a polynomial, and therefore each irreducible rep- 
resentation in this space can be characterized by an integer j = 0, 1,2, . . ., 
specifying the degree of the polynomials that span this representation. At this 
point, we cannot rule out the possibility that some other integers are needed 
to fully specify all irreducible representations in Sd and therefore we will need 
below another set of indices to complete the specification. 

We can now choose a basis of polynomials {Yj^„{u)} that span all the irre- 
ducible representations of SO(d) over Sd- The index a counts all integers other 
than j needed to fully specify all irreducible representations, and in addition, 
it labels the different functions within each irreducible representation. 

Let us demonstrate this construction in two and three dimensions. In two 
dimensions a is unneeded since all the irreducible representation are one- 
dimensional and are spanned by Yj{u) = e^^'^ with being the angle between 
u and the the vector ei = (1,0). Any rotation of the coordinates in an angle 
00 results in a multiplicative factor e^'^° . It is clear that Yj{u) is a polynomial 
in u since Yj{u) = [u ■ p]^ where p = (1, i). In three dimensions a = m where 
m takes on 2j + 1 values m — —j, + Here Yj^rn e*"*'^Pj"(cos ^) 

where and 9 are the usual spherical coordinates, and Pp is the associated 
Legendre polynomial of degree j — m. Obviously we again have a polynomial 
in u of degree j. 

We now wish to calculate the action of the Laplacian operator with respect 
to u on the Yj^a-{u). We prove the following identity: 

u^d^da^^iu) = + d- 2)Yj,„{u) . (B.2) 

One can easily check that for d = 3 (B.2) gives the factor j(j + 1), well known 
from the theory of angular-momentum in Quantum Mechanics. To prove this 
identity for any d, note that 

\u\'^-^d^\u\%,a{u)) = . (B.3) 
This follows from the fact that the Laplacian is an isotropic operator, and 
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therefore is diagonal in the Yj^cr- The same is true for the operator 
But this operator results in a polynomial in u of degree j — 2, which is spanned 
by 5^',(t' such that j' < j — 2. Therefore the RHS of (B.3) must vanish. Ac- 
cordingly we write 

d'^\u^\Yj^^{u) + 2d°'\u^\d°'Yj^^ + \u^\d%^^{u) = . (B.4) 

The second term vanishes since it contains a radial derivative u'^da operating 
on Yj^„{u) which depends on u only. The first and third terms, upon elementary 
manipulations, lead to (B.2). 

Having the 5^>('u) we can now construct the irreducible representations in 
the space of constant tensors. The method is based on acting on the Yj^„{u) 
with the isotropic operators u", 9" and (5"'^ Due to the isotropy of the above 
operators, the behavior of the resulting expressions under rotations is similar 
to the behavior of the scalar function we started with. For example, the tensor 
fields S°'^Yj^a-{u): d°'d^Yj^„{u) transform under rotations according to the (j, a) 
sector of S6(d). 

Next, we wish to find the basis for the irreducible representations of the space 
of constant and fully symmetric tensors with n indices. We form the basis 

^ai..,.„ ^ _ _ _ 3<n. (B.5) 

Note that when j and n are even 5"^'^^'"" no longer depends on u, and is 
indeed fully symmetric by construction. Simple arguments can also prove that 
this basis is indeed complete, and spans all fully symmetric tensors with n 
indices. Other examples of this procedure for the other spaces are presented 
directly in the text. 

Finally let us introduce two identities involving the Bnj^a- The first one is 

Zn^j = [n{n + d-2)- j{j + d-2)]. (B.7) 
It is straightforward to derive this identity using (B.2). The second identity is 
Y: = ^^j:;-'''" , 3<n-2. (B.8) 

This identity is proven by writing -u" in (B.5) as u'^u"'^^, and operating with 
the derivative on u^. The term obtained as u^d""^ . . . d'^"u^~'^Yj^„{u) vanishes 
because we have n derivatives on a polynomial of degree n — 2. It is worthwhile 
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noticing that these identities connect tensors from two different spaces. The 
space of tensors with n indices and the space of tensors with n — 2 indices. 
Nevertheless, in both spaces, the tensors belong to the same (j, a) sector of 
the SO(d) group. This is due to the isotropy of the contraction with 5"i"2 in 
the first identity, and the contraction with 6°'^°'^ in the second identity. 



C Full Form for the j — 2 Contribution for the Homogeneous Case 



In this appendix we focus on the decomposition of second order tensorial 
structure functions up to j = 2. For this purpose we define: 

The J = is well-known and given explicitly by: 



CqT 



''0 



(2) 



a B 



(C.l) 



where Q ~ 0.68 is the known universal scaling exponent for the isotropic 
contribution and Cq is an unknown coefficient that depends on the boundary 

conditions of the flow. For the j = 2 sector which is the lowest contribution to 
anisotropy to the homogeneous structure function, the m = (axisymmetric) 
terms were derived from constraints of symmetry, even parity (because of 
homogeneity) and incompressibility on the second order structure function 
[29] 



S j=2,m=o{'''^) 



A2) 



(Cf ^ - 2)5"/^ - Cf + 6)5-/^1^ + 2Cf (Cf - 2) 



^(2)/ . (2) _ .J. 



(C.2) 



(2) 

where Q is the universal scaling exponent for the j = 2 anisotropic sector 
and a and b are independent unknown coefficients to be determined by the 
boundary conditions. We would now like to derive the remaining m — ±1, 
and m = ±2 components 



'^2m — Z^"'q,2,mr ^ ^q,2,m{'^) 
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As usual the q label denotes the different possible ways of arriving at the same 
j and runs over all such terms with the same parity and symmetry (a conse- 
quence of homogeneity and hence the constraint of incompressibility) . In all 
that follows, we work closely with the procedure outhned in [23] . Following the 
convention in [23] the ^'s to sum over are q — {1, 7, 9, 5}. The incompressibility 
condition dau" = coupled with homogeneity can be used to give relations 
between the agjm for a given {j, m). That is, for j = 2, m = — 2 . . . 2 



(Cf ^ - 2)ai,2,m + 2(Cf - 2)a7,2m + (Cf + 2)a9,2,m = (C.3) 

Ol,2,m + (C2 + 3)a7,2,m + C2 Ct5,2,m = 0. 

We solve the above equations in order to obtain a^^2,m and a7,2m in terms of 
linear combinations of ai,2,m and aQ^2m- 



«i,2..([Cf '1^ - Cf ' - 2) + a9.2.„([Cf ]^ + 5d^' + 6) 
2Cf'(cf ' -2) 

_ ai,2m(2 — C-P) — a9,2,m(2 + C2^'') 
Q7,2,m — -—py\ ; ■ 

2(d') - 2) 

Using the above constraints on the coefficients, we are now left with a linear 
combination of just two linearly independent tensor forms for each m 



-^2^ = «9,2,.r^^ (2 + CnB?iJr) + 2Cf (Cf - 2)BSl^Jr) 

+ i[C?^?-CP-2)Bli^{r)]. (C.5) 



The task remains to find the explicit form of the basis tensor functions -B^ f m(^) > 

ge{l,7,9,5},me{±l,±2} 

•5it™(^)=r-25"Vy2„(f) 

• Brimir) = r-2[r"9'3 + rf^d''yY2m{f) 

•Bf^„Xr)=r-^r^r^-^Y2^{r) 

.Bti^{f)^d'^d^T%m{r) 

We obtain the m = {±1, ±2} basis functions in the following derivation. We 
first note that it is more convenient to form a real basis from the r'^Y2m{r) since 
we ultimately wish to fit to real quantities and extract real best-fit parameters. 
We therefore form the r'^Y2k{r) [k = —1, 0, 1) as follows: 
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rVa o(f) = rVs o(f) = cos' 9 ^ 



— r 



r'l^2+i(r) 



2 

2 (^cos — i sin 0) cos 6* sin + (cos + i sin 0) cos sin ^ 

2 

: r' cos ^ sin ^ cos = rari 
2 (cos — i sin 0) cos 6* sin 6' — (cos + i sin 0) cos 6* sin ^ 



-2i 

cos sin 6 sin = r3r2 



r 



2 Y2 2{r)-Y2.2{r) 
2i 

2 (cos 2(f) + i sin 20) sin' 9 — (cos 2(f) — i sin 20) sin' ^ 



2i 

r' sin 20 sin' 9 = 2rir2 



2 (cos 20 + i sin 20) sin' 9 + (cos 20 — i sin 20) sin' 9 
~^ 2 

= r' cos 20 sin' 9^rl-rl (C.6) 

This new basis of r'^Y2k{r) is equivalent to using the r'ljm(r) themselves as 
they form a complete, orthogonal (in the new /c's) set. We omit the normal- 
ization constants for the spherical harmonics for notational convenience. The 
subscripts on r denote its components along the 1 (m), 2 {p) and 3 {n) di- 
rections, m denotes the shear direction, p the horizontal direction parallel to 
the boundary and orthogonal to the mean wind direction and n the direction 
of the mean wind. This notation makes it simple to take the derivatives when 
we form the different basis tensors and the only thing to remember is that 



a°ri = 9°^(r . m) = m'^ 
= (9"(r -p) 

= a"(r . n) (C.7) 



We use the above identities to proceed to derive the basis tensor functions 
Bii-i{r) ^r-H'^^ir • n){r - m) 

Bf2,-iir) =r-'[(r'*m'' + r%")(r • n) + {r^'n^ + r^n°)(r • m)] 
Bgi-iir) ^r-\'^r^(r - n)(r - m) 
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Bf^^,{r)^r-h'^rP{r.n){r.p) 
Bii-2{'^) = 2r-'5"^(r • m)(r • p) 

B7i-2{r) = 2r-^[(rV + r'^p")ir - m) + {r"m'^ + r^m")(r • p)] 
Bg2 -2(*^) = 2r~^r°r^(r • m)(r • p) 
Bti-2{r) = 2(mV + ^^V) 

57,2,2(»') =2r"^[(r"m^ + rV°)(r • m) - (rV + rV)(T- ' P)] 
B9i2ir) =r-^r''r^[(r • m)^ - (r • p)^] 

^5i2(r=) = 2(m"m^ - pV) (C.8) 



Note that for each dimension k the tensor is bihnear in some combination of 
two basis vectors from the set m, p and n. Substituting these tensors forms 
into Eq. C.5 we obtain the full tensor forms for the j — 2 non-axisymmetric 
terms, with two independent coefficients for each k. 



.(2) 



'S'j=2,fe=l(^, 



■■ag,2,-ir'^ [ - C2^H2 + C?V"^[K^'^ + r^m")(r • n) 
+ (r^n'^ + r^n")(r • m)] + 2Ci^\cP - 2)r-^r"r'^{r • n)(r • m) 
+ ([Cf ]' + 5Cf ^ + 6)(n"m^ + n^m")" 
+ ai,2,-ir^^'^ [2C?^(C2'^ - 2)r-25"'^(r • n)(r • m) 
-d^^ld^^ - 2)r-2[(r"m'^ + r^m")(r • n) + (r^n'^ + r'^n")(r • m)] 
+ ([d'^]' - d'^ - 2)(n"m^ + n/^m")" 
= a9,2,ir^2'' [ - (P{2 + d'V"'[(^V + rV)(r' • n) 
+ (r"n^ + r^n")(r • p)] + 2d'\d'^ - 2)r-Vr^(r • n)(r • p) 
+ ([dy + Sd'^ + 6)(nV + riV)" 



+ ai,2,i/2'' 2Cr(Cr - 2)r-^5«'^^(r • n)(r ■ p 



.(2).^(2) 



-2 rati/ 



- d'Vf^ - 2)r^"[(rV + rV)(r- • n) + (r"n^ + r''n")(r • p; 



(2)l2 



a(2) 
S2 



2)(nV + ^V) 



'5g2,.=-2(r-) =a9,2,-2r^^" f - 2d'^(2 + d'V-'[(rV + rV)(r- • m) 



.(2) 



(2)^ 



+ (r"m^ + r'^m")(r • p)] + 2Cr(Cr - 2)r-Vr^(r • p)(r • m) 



(2)/^(2) 
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+ ([Cff + 5Cf + 6)(mV + mV) 
+ ai,2,-2/^'' [2Cf (Cf - 2)r-2<5"'^(r • m)(r • p) 

- 2Cf ^(C^^^ - 2)r-2[(rV + rV)(^ • ^^i) + (r"m^ + r^m")(r • p)] 
+ 2([Cf - Cf ^ - 2)(mV + ^V)" 

^jf2,fe=2(r-) =a9,2,2r^^'' [ - 2Cf ^2 + C?V-^[(r«m^ + r^m'^)(r • m) 

- (rV + rV)(r- • p)] + 2Cf ^ - 2)r-V°r^[(r • m)' - (r • p)'] 
+ 2([d'^]' + Sd'^ + 6)(m"m^ - pV) 

+ a,2,2r^^^72Cf(Cf -2)r-V^[(r. 



m)^ — (r • p)^] 



- 2Cf ^ - 2)r-"[(r°m'' + r''m")(r • m) - (rV + rV)(r" • P)] 



+ 2([C2'^]' - - 2){m''m^ - p^p") 



(C.9) 



Now we want to use this form to fit for the scahng exponent Q in the structure 
function S^^{r) from data set I where a = j3 = Z and the azimuthal angle of 
r in the geometry is (j) = tt/2. 



5 5^2,^=1 0, 71/2) = a,,2,ir^^"\-2CP{C? + 2) sin^cos^ + 2Cf (Cf - 2) cos^ ^sin^] 
S%,^,^_,{r,9,n/2)^0 

S%,^,^,{r, 9, n/2) = a9,2,2r^^'V2Cf (Cf^ - 2) cos^ ^sin^ 9] + ai,2,2r^^'V2Cf^(Cf - 2) sin^ ^] 



We sec that choosing a particular geometry eliminates certain tensor contri- 
butions. In the case of set I we are left with 3 independent coefficients for 
m 7^ 0, the 2 coefficients from the m — contribution (Eq. C.2), and the sin- 
gle coefficient from the isotropic sector C.l, giving a total of 6 fit parameters. 
The general forms in C.9 can be used along with the k = (axisymmetric) 
contribution C.l to fit to any second order tensor object. For convenience, the 
table shows the number of independent coefficients that a few different exper- 
imental geometries we have will allow in the j — 2 sector. It must be kept in 
mind that these forms are to be used only when there is known to be homo- 
geneity. If there is inhomogeneity, then we cannot apply the incompressibility 
condition to provide constraints in the various parity and symmetry sectors 
and we must in general mix different parity objects, using only the geometry 
of the experiment itself to eliminate any terms. 
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K 


^ = 7r/2,Q! = /3 = 3 


^ = 0,a = /3 = 3 


(/, = 0,a = /3 = 1 


= O,Q! = 3,/3 = 1 


(7 /= u 


(7 — U 




(7 — U 


(7 /= u 


u — u 


(7 /= u 


(7 — U 


n 
u 


9 


9 


9 


9 


9 


9 
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_i 

J. 








1 





1 





2 


2 


1 


1 
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2 


2 





2 





2 


2 


2 





Total 


5 


2 


5 


2 


5 


4 


6 


2 


Table C 


.1 



The number of free coefficients in the j = 2 sector for homogeneous turbulence and 
for different geometries 



D The j=l Component in the Inhomogeneous Case 



D.l Antisymmetric Contribution 



We consider the tensor 

T"^(r) = {u'^{x + r) - u'^{x)){u^{x + r) + uf^ix))). 

This object is trivially zero for a = j3.\n the experimental setup, we measure 
at points separated in the shear direction and therefore have inhomogencity 
which makes the object of mixed parity and symmetry. We cannot apply 
the incompressibility condition in same parity /symmetry sectors as before to 
provide constraints. Wc must in general use all 7 irreducible tensor forms. This 
would mean fitting for 7 x 3 = 21 independent coefficients plus 1 exponent (,{ 
in the anisotropic sector, together with 2 coefficients in the isotropic sector. 
In order to pare down the number of parameter we are fitting for, we look at 
the antisymmetric part of T°'^{r) 

2^a/3(^) ^ — ^ {u^(^x)u^{x + r)) - {u^{x)u''{x + r)) 

which will only have contributions from the antisymmetric j — 1 basis tensors. 

These are 

• Antisymmetric, odd parity 

Bf,,m = r-^[r"9^ - r^9'^]ryi,^(r) (D.l) 



• Antisymmetric, even parity 
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B^i^^r-h^P^d,rY,,m{r) (D.2) 

As with the j = 2 case we form a real basis rYi^k{r) from the (in general) 
complex rYi^jn{r) in order to obtain real coefficients in the fits. 



rYi^k=o{r) ^ rYi,o(r) = r cos6' = rg 

^ ... yi,i(r) + Fi,i(r) . , . , 
rYi^k=i[r) — r — ■ — — ■ = r sm ^ sm (p — r2 

^ n,-i(r)-yi,i(r) . , , 

And the final forms are 



-03,1,0 


(f) 


= r" 




-04,1,0 


(f) 


= r 


-2ga/3M^^(^ • n) 


2,1,0 


(f) 


= r" 




^3,1,1 


(f) 


= r 




-04,1,1 


(r) 


= r" 


-V^%(r.p) 


-02,1,1 


(r) 


= r" 




-03,1,-1 


(r) 


= r" 




nap 
-04,1,-1 


(r) 


= r 


-2e°^'^r^(r • m) 


Tjal3 
-02,1,-1 


(r) 


= r' 





(D.3) 

Note: For a given /c the representations is symmetric about a particular axis 
in the coordinate system chosen (l=m (shear), 2=p (horizontal), 3=n (mean- 
wind)) We now have 9 independent terms and we cannot apply incompress- 
ibility in order to reduce the number of independent coefficients in the fitting 
procedure. We use the geometric constraints of the experiment to do this. 
• = (vertical separation), a = 3, /? = 3 



4\,o(^,^,0 = O) = -sin^ 
4\i(r,^^,0 = O) = l 
Bl\^_^{r,9,(t> = Q)^cose (D.4) 

There are no contributions from the reflection-symmetric terms in the j = 
isotropic sector since these are symmetric in the indices. The helicity term in 
J = also doesn't contribute because of the geometry. So, to lowest order 

r"'^(r) =r/ii(r) =a3,i,o(r)(-sin^) + a2,i,i(r) + a3,i,-i(r)cose 
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We have 3 unknown independent coefficients and 1 unknown exponent to fit 
for in tlie data. 



D.2 Symmetric Contribution 
We consider tlie structure function 



in tfie case wfiere we fiave fiomogeneous flow. Tfiis object is symmetric in the 
indices by construction, and it is easily seen that homogeneity implies even 
parity in r: S°'^{r) = S^"{r) and S'"'^(— r) = S°'^{r). We reason that this 
object cannot exhibit a j = 1 contribution from the 5*0(3) representation in 
the following manner. Homogeneity allows us to use the incompressibility con- 
dition: daS°'^ — and dpS"'^ — 0, separately on the basis tensors of a given 
parity and symmetry in order to give relationships between their coefficients. 
For the even parity, symmetric case we have for general J > 2 just two basis 
tensors and they must occur in some linear combination with incompressibil- 
ity providing a constraint between the two coefficients. However, for j' = 1 we 
only have one such tensor in the even parity, symmetric group. Therefore, by 
incompressibility, its coefficient must vanish. Consequently, we cannot have a 
j = 1 contribution for the even parity (homogeneous), symmetric structure 
function. Now, we consider the case as available in experiment when r has 
some component in the inhomogeneous direction. Now, it is no longer true 
that S°''^{r) is of even parity and moreover it is also not possible to use in- 
compressibility as above to exclude the existence of a j = 1 contribution. We 
must look at all j = 1 basis tensors that are symmetric, but not confined to 
even parity. These are 
• Odd parity, symmetric 



S^^{r) = ((«"(cc + r) - u''{x)){u^{x + r) - u^{x))) 





rd°'d^rYik(r) = 



(D.5) 



• Even parity, symmetric 
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We use the real basis of r~^Yik{r) which are formed from the r^^Yim{r) . Both 
-^5ifc(^) -^6ifc(^) vanish because of the taking of the double derivative 
of an object of single power in r. We thus have 4 different contributions to 
symmetric J = 1 and each of these is of 3 dimensions {k = —1, 0, 1) giving in 
general 12 terms in all. 



<,o(r) 


= r~ 


-^5"^(r • n) 








— r~ 










— r' 


-3j.a^P(^r • n) 






Btioir) 


= r~ 




•P) 


-(rV + ^V)(T"'^«) 


BtUr) 


— r' 


-^S'^^ir ' p) 






B?Ur) 


— r' 










= r' 








Bf^,^ir) 


= r~ 




• n) 


- (rV + r^n")(r -m 


Bfi,-iir) 


= r' 


-^r'^(r • m) 






Brl-iir) 


= r' 










= r' 


-'^r"r'^{r • m) 








= r" 


"2[(rV + rV)(r" • 


n)- 


(r"n'^ + r'^n")(r .p)] 



These are all the possible j — 1 contributions to the symmetric, mixed parity 
(inhomogeneous) structure function. 

For the experimental setup II, we want to analyze the inhomogeneous structure 
function in the case a — (3 = 3, and azimuthal angle = (which corresponds 
to vertical separation) and we obtain the basis tensors 



BfioiO) 


= cos 9 


Bjxoi^) 


= 2cos^ 


^9lo(^) 


= cos^ 6 


i?8,ll(^) 


= —2 cos 9 sin 9 


Bfi-AO) 


— sin 9 


BlUi^) 


— cos^ 9 sin 9 



Table D.l gives the number of free coefficients in the symmetric j = 1 sec- 
tor in the fit to the inhomogeneous structure function for various geometric 
configurations. 
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k 


^ = 0,a = /3 = 3 


(f) = 0^a = P = l 


(l) = 0,a = 3,P = 1 


(9 7^ 


6' = 


9 y^O 


9 = 


9^0 


9 = 





3 


3 


2 


1 


2 





1 


1 





1 











-1 


2 





3 





2 


1 


T<)(al 


G 


O 

o 


G 


1 


1 


1 



Table D.l 

The number of free coefficients in the symmetric j = 1 sector for inhomogeneous 
turbulence and for different geometries. 

E The Matrix Form of the Operator of the Linecir Pressure Model 

Using the basic identities of the ^jm{r) functions (see [23]), 

a short calculation yields: 



;CC"(r) = K>'''{r)d^d,C"{r) = Dx' f2c'/ + 2(2 + e)^ - (2 + e)(j + l)(j + 2) 
+Dx 

Therefore, in matrix notation, the Kraichnan operator can be written as: 



r r 
24' + 2(2 + 6)^ + 2(2 + - 2(2 + e)j{j - l)%]B%^{r) . 



2 



I— 




2Dr' 



1 



1 






14' 



1 0' + i)(i + 2) 




Letting 



(E.l) 



(E.2) 
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and applying a Laplacian to PT'*, we get, 



j4' + j^+j(2j-1)^-j(j + 1)^ 



+ 



r r 



Hence in matrix notation, 



"(;;)-(r;)(:;)^H-."r;')(::;) 



Now that the matrix forms of the Kraichnan operator and of the Laplacian of 
the projection operator have been found, we can combine these two results to 
find the matrix form of the LHS of Eq. (100). To this aim let us define. 




and from Eq. (E.1,E.3) we get 




where the number in parenthesis denotes the order of the derivative. The 
matrices Mj are given by: 



M4 = P2IK2 , (E.4) 

M3 = 2eP2K2 + P2l^i + P1K2 , 

M2 = e(e - 1)^2^2 + 2(e - l)P2Ki + P2K0 

+ePiK2 + PiKi + P0K2 , 
Ml = (e - l)(e - 2)P2Ki + 2(e - 2)P2Ko 

+(e - l)PiKi + PiKo + PoKi , 
Mo = (e - 2)(e - 3)P2Ko + (e - 2)PiKo + PqKq . 
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To find the RHS of Eq. (100) we expand the "forcing" A°'{r) in terms of the 
spherical vectors Bij^, B2jm, 



(E.5) 



and applying a Laplacian we find the matrix form of d A°'(r): 



Or+f/i-o-+i)(j+2)M ^ 




If 



lfl-j{j - l);^/2 




(E.6) 



At this point it is worthwhile to remember that the forcing term A°'{r/L) is 
assumed to be analytic. As a result for r/L <^ 1 its leading contribution in 
the (j, m) sector is proportional to d°'r^Yjjn{r) ~ r^~^. However d'^A°'{r/L) 
is also analytic, and must therefore also scale like r^~^ for small r, instead of 
like r^~^ which could be the naive dimensional guess. 

To proceed we restrict ourselves to finding the solution in the inertial range 
and beyond. In these ranges the dissipative term Kd'^d^C"{r) is negligible and 
can be omitted, thus reaching Eq. (101) for Ci(r) and C2(r). 
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